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Preface 


Students commonly find difficulty with problems in fluid mechanics. They 
may misunderstand what is required or misapply the solutions. This book is 
intended to help. It is a collection of problems in elementary fluid mechanics 
with accompanying solutions, and intended principally as a study aid for under- 
graduate students of chemical engineering — although students of all engineering 
disciplines will find it useful. It helps in preparation for examinations, when 
tackling coursework and assignments, and later in more advanced studies of the 
subject. In preparing this book I have not tried to replace other, fuller texts on 
the subject. Instead I have aimed at supporting undergraduate courses and 
academic tutors involved in the supervision of design projects. 

In the text, worked examples enable the reader to become familiar with, and 
to grasp firmly, important concepts and principles in fluid mechanics such as 
mass, energy and momentum. The mathematical approach is simple for anyone 
with prior knowledge of basic engineering concepts. I have limited the prob- 
lems to those involving incompressible, Newtonian fluids and single-phase 
flow through pipes. There is no attempt to include the effects of compressible 
and non-Newtonian fluids, or of heat and mass transfer. I also held back from 
more advanced mathematical tools such as vectorial and tensorial mathematics. 

Many of the problems featured have been provided by university lecturers 
who are directly involved in teaching fluid mechanics, and by professional 
engineers in industry. I have selected each problem specifically for the light it 
throws on the fundamentals applied to chemical engineering, and for the confi- 
dence its solution engenders. 

The curricula of university chemical engineering degree courses cover the 
fundamentals of fluid mechanics with reasonable consistency although, in 
certain areas, there are some differences in both procedures and nomenclature. 
This book adopts a consistent approach throughout which should be recogniz- 
able to all students and lecturers. 

Thave tailored the problems kindly contributed by industri: 
commercial secrets and to ensure that the nature of each problem is clear. There 























to safeguard 


is no information or detail which might allow a particular process or company 
to be recognized. All the problems use SI units. As traditional systems of units 
are still very much in use in industry, there is a table of useful conversions. 
Fluid mechanics has a jargon of its own, so I have included a list of definitions. 

There are nine chapters. They cover a range from stationary fluids through 
fluids in motion. Each chapter contains a selected number of problems with 
solutions that lead the reader step by step. Where appropriate, there are prob- 
lems with additional points to facilitate a fuller understanding. Historical refer- 
ences to prominent pioneers in fluid mechanics are also included. At the end of 
each chapter a number of additional problems appear; the aim is to extend the 
reader's experience in problem-solving and to help develop a deeper under- 
standing of the subject. 

I would like to express my sincere appreciation to Dr Robert Edge ( formerly 
of Strathclyde University), Mr Brendon Harty (Roche Products Limited), Dr 
Vahid Nassehi (Loughborough University), Professor Christopher Rielly 
(Loughborough University), Professor Laurence Weatherley (University of 
Canterbury), Dr Graeme White (Heriot Watt University), Mr Martin Tims 
(Esso UK plc) and Miss Audra Morgan (IChemE) for their assistance in 
preparing this book. I am also grateful for the many discussions with profes- 
sional engineers from ICI, Esso and Kvaerner Process Technol 

The text has been carefully checked. In the event, however, that readers 
uncover any error, misprint or obscurity, I would be grateful to hear about it. 
Suggestions for improvement a 















so welcome. 


Carl Schaschke 
April 2000 


List of symbols 


The symbols used in the worked examples are defined below. Where possible, 
they conform to consistent usage elsewhere and to international standards. SI 
units are used although derived SI units or specialist terms are used where 
appropriate. Specific subscripts are defined separately. 


Roman Term Sl or preferred unit 
area of pipe or orifice m? 
area of channel or tank m? 
breadth of rectangular weir m 
constant — 


velocity of sound 
Chézy coefficient 





coefficient — 
concentration gr! 
diameter m 
impeller diameter m 
fraction — 
friction factor — 
depth of body below free surface m 
force N 
gravitational acceleration ms? 
head m 
slope of channel = 
constant E 
fundamental dimension for length — — 
length 


mass loading 

mass 

mass flowrate 

mean hydraulic depth 
fundamental dimension for mass — 
channel roughness ms 





SEPP POOH TW DA ABWORC AAA wes 


OYZ 


t^ mmc 
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m 
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number of pipe diameters 
rotational speed 
specific speed 
pressure 

power 

wetted perimeter 
volumetric flowrate 
radius 

frictional resistance 
radius 

depth 

suction specific speed 
thickness of oil film 
time 

fundamental dimension for time 
torque 

velocity 

volume 

width 

work 

principal co-ordinate 
distance 

principal co-ordinate 
distance 

principal co-ordinate 
static head 


ratio of pipe to throat diameter 
film thickness 
finite difference 
absolute roughness 
pump efficiency 
angle 

friction factor 
dynamic viscosity 
kinematic viscosity 
3.14159 

density 

surface tension 
shear stress 

friction factor 
angular velocity 


rps 
m3/45-3/2 


Nm? 


m 
ms! 





radians s~! 





Fluid mechanics and 
problem-solving 


Fluid mechanics forms an integral part of the education of a chemical engineer. 
The science deals with the behaviour of fluids when subjected to changes of 
pressure, frictional resistance, flow through pipes, ducts, restrictions and 
production of power. It also includes the development and testing of theories 
devised to explain various phenomena. To the chemical engineer, a knowledge 
of the behaviour of fluids is of crucial importance in cost-effective design and 
efficient operation of process plant. 

Fluid mechanics is well known for the large number of concepts needed to 
solve even the apparently simplest of problems. It is important for the engineer 
to have a full and lucid grasp of these concepts in order to attempt to solve prob- 
lems in fluid mechanics. There is, of course, a considerable difference between 
studying the principles of the subject for examination purposes, and their appli- 
cation by the practising chemical engineer. Both the student and the profes- 
sional chemical engineer, however, require a sound grounding. It is essential 
that the basics are thoroughly understood and can be correctly applied. 

Many students have difficulty in identifying relevant information and funda- 
mentals, particularly close to examination time. Equally, students may be hesi- 
tant in applying theories covered in their studies, resulting from either an 
incomplete understanding of the principles or a lack of confidence caused by 
unfamiliarity. For those new to the subject, finding a clear path to solving a 
problem may not always be straightforward. For the unwary and inexperienced, 
the opportunity to deviate, to apply incorrect or inappropriate formulae or to 
reach a mathematical impasse in the face of complex equations, is all too real. 
The danger is that the student will dwell on a mathematical quirk which may be 
specific purely to the manner in which the problem has been (incorrectly) 
approached. A disproportionate amount of effort will therefore be expended on 
something irrelevant to the subject of fluid mechanic: 

Students develop and use methods for study which are dependent on their 
own personal needs, circumstances and available resources. In general, 
however, a quicker and deeper understanding of principles is achieved when a 
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problem is provided with an accompanying solution. The worked example is a 
recognized and widely-used approach to self-study, providing a clear and 
logical approach from a distinct starting point through defined steps, together 
with the relevant mathematical formulae and manipulation. This method bene- 
fits the student by appreciation of both the depth and complexity involved in 
reaching a solution. 

While some problems in fluid mechanics are straightforward, unexpected 
difficulties can be encountered when seemingly similar or related simple prob- 
lems require the evaluation of a different but associated variable. Although the 
solution may require the same starting point, the route through to the final 
answer may be quite different. For example, determining the rate of uniform 
flow along an inclined channel given the dimensions of the channel is straight- 
forward. But determining the depth of flow along the channel for given parame- 
ters in the flow presents a problem. Whereas the former is readily solved analyti- 
cally, the latter is complicated by the fact that the fluid velocity, flow area and a 
flow coefficient all involve the depth of flow. An analytical solution is no longer 
possible, thus requiring the use of graphical or trial and error approaches. 

There are many similarities between the governing equations in heat, mass 
and momentum transport and it is often beneficial to bring together different 
branches of the subject. Other analogies between different disciplines are also 
useful, although they must be applied with care. In fluid mechanics, analogies 
between electrical current and resistance are often used, particularly in dealing 
with pipe networks where the splitting and combining of lines can be likened to 
resistors in parallel and in series. 

Some applications of fluid mechanics require involved procedures. 
Selecting a pump, for example, follows a fairly straightforward set of well- 
defined steps although the lengthy procedure needed can become confusing. It 
is important to establish the relationship between the flowrate and pressure, or 
head, losses in the pipework connecting process vessels together. With fric- 
tional losses due to pipe bends, elbows and other fittings represented by either 
equivalent length pipe or velocity heads, pumping problems therefore require 
careful delineation. Any pump calculation is best reduced to the evaluation of 
the suction pressure or head and then of the discharge head; the difference is the 
delivery head required from the pump. For a sizing calculation, all that is really 
needed is to determine the delivery head for the required volumetric flowrate. 
As in many process engineering calculations dealing with equipment sizing, the 
physical layout plays an important part, not only in standardizing the method 
for easy checking but also in simplifing the calculations. Obviously there will 
s requiring more detail but, with a bit of attention, such deviations from 
practice can easily be incorporated. 
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Finally, the application of fluid mechanics in chemical engineering today 
relies on the fundamental principles largely founded in the seventeenth and 
eighteenth centuries by scientists including Bernoulli, Newton and Euler. 
Many of today’s engineering problems are complex, non-linear, three- 
dimensional and transient, requiring interdisciplinary approaches to solution. 
High-speed and powerful computers are increasingly used to solve complex 
problems, particularly in computational fluid dynamics (CFD). It is worth 
remembering, however, that the solutions are only as valid as the mathematical 
models and experimental data used to describe fluid flow phenomena. There is, 
for example, no analytical model that describes precisely the random behaviour 
of fluids in turbulent motion. There is still no substitute for an all-round under- 
standing and appreciation of the underlying concepts and the ability to solve or 
check problems from first principles. 





Contents 


Preface 
List of symbols. 
Fluid mechanics and problem-solving 


1 Fluid statics 

Introduction 

1.1 Pressure at a point 

12 Pressure within a closed vessel. 

13 Forces within a hydraulic ram 

l4 Liquid-liquid interface position in a solvent separator 

15 LiquidHiquid interface measurement by differential pressure 
1.6 Measurement of crystal concentration by differential pressure 
17 Pressure within a gas bubble 

18 Pressure measurement by differential manometer 

19 Pressure measurement by inverted manometer 

1.10 Pressure measurement by single leg manometer 

1.11 Pressure measurement by inclined leg manometer 

1.12. Archimedes’ principle 

1.13. Specific gravity measurement by hydrometer 

1.1 Transfer of process liquid to a ship. 

Further problems 


2 Continuity, momentum and energy 
Introduction 

2. Flow in branched pipes 

22 Forces on a U-bend 

23 Pressure rise by valve closure 

24' The Bernoulli equation 

25 Pressure drop due to enlargements 


35 
35 
36 
38 


2 


xi 






























































26 Pipe entrance head loss 
27 Force on a pipe reducer 
28 Vortex motion 

29 Forced and free vortices 
Further problems. 


3  Laminar flow and lubrication 

Introduction 

3. Reynolds number equations 

32 Laminar boundary layer 

33 Velocity profile in a pipe 

34 Hagen-Poiseuille equation for laminar flow in a pipe 

35 Pipe diameter for laminar flow 

3,6 Laminar flow through a tapered tube 

37 Relationship between average and maximum velocity in a pipe 
3,8 Relationship between local and maximum velocity in a pipe 
39 Maximum pipe diameter for laminar flow 

3.10 Vertical pipe flow 

3.11 Film thickness in a channel 

3.12. Flow down an inclined plate 

3.13. Flow down a vertical wire 

3.14 Flow and local velocity through a gap 

3.15 Relationship between local and average velocity through a gap. 
3.16 Relationship between average and maximum velocity through a gap 
3.17 Shear stress for flow through a gap 

3.18 Flat disc viscometer 

3.19 Torque on a lubricated shaft 

3.20 Lubricated collar bearing 

Further problems 


4 Dimensional analysis 

Introduction 

41 Flow through an orifice 

42 Flow over notches 

43 Scale-up of centrifugal pumps 

44 Frictional pressure drop for turbulent flow in pipes 
45 Scale model for predicting pressure drop in a pipeline 
Further problems 


5 Flow measurement by differential head 
Introduction 





4 
49 
5l 
53 
57 


61 
61 
63 


8888 


n 
D 
75 
7 
78 
79 
81 
82 
84 
86 
87 
88 
89 
9 
95 
95 


To 


13 
13 





5.1 Pitottube 

52 Pitot traverse 

53 Horizontal venturi meter 

54 Orifice and venturi meters in parallel 

55 Venturi meter calibration by tracer dilution 

56 Differential pressure across a vertical venturi meter 
57 Flow measurement by orifice meter in a vertical pipe 
5,8 Variable area flowmeter 

59 Rotameter calibration by venturi meter 

Further problems 


6 Tank drainage and variable head flow 
Introduction 

6.1 Orifice flow under constant head 

62 Coefficient of velocity 

63 Drainage from tank with uniform cross-section 
64 Tank drainage with hemispherical cross-section 
65 Tank drainage with cylindrical cross-section 

6.6 Drainage between two reservoirs 

6.7 Tank inflow with simultaneous outflow 

6.8 Instantaneous tank discharge 

69 Instantaneous tank inflow with outflow 

6.10 Tank drainage through a horizontal pipe with laminar flow 
Further problems 


7 Open channels, notches and weirs 
Introduction 

7.1 Chézy formula for open channel flow 

7.2. Flow ina rectangular open channel 

73 Depth of flow in a rectangular channel 

7.4 Economical depth of flow in rectangular channels 
75 Circular channel flow 

76 Maximum flow in circular channels 

73 Weirs and rectangular notches 

78 Depth of a rectangular weir 

79 Instantaneous flow through a rectangular weir 
7.10 Flow through a triangular notch 

7.11 Tank drainage through a V-notch 

7.12. Flow through a trapezoidal notch 

Further problems 


15 
n 
19 
22 
124 
126 
128 
130 
134 
136 


41 
Ml 
143 
145 
147 
149 
151 
153 
155 
157 
159 
161 
163 


167 
167 
169 
WW 
173 
175 
177 
179 
180 
182 
184 
186 
188 
189 
190 


xiii 





8 Pipe friction and turbulent flow 

Introduction. 

8.1 Economic pipe diameter 

82 Head oss due to friction. 

83 General frictional pressure drop equation applied to laminar flow 
84 Blasius equation for smooth-walled pipes 

8.5 Prandtl's universal resistance equation for smooth-walled pipes 
8.6 Pressure drop through a rough-walled horizontal pipe 

87 Discharge through a siphon 

88 Flow through parallel pipes 

89 Pipes in series: flow by velocity head method 

8.10 Pipes in series: pressure drop by equivalent length method 

8.11 Relationship between equivalent length and velocity head methods 
8.12. Flow and pressure drop around a ring main 

8.135 Tank drainage through a pipe with turbulent flow 

8.14 Turbulent flow in non-circular ducts 

8.15 Head loss through a tapered section 

8,16 Acceleration of a liquid in a pipe 

Further problems 


9 Pumps 

Introduction 

9,1 Centrifugal pumps 

92 Centrifugal pump matching 

93 Centrifugal pumps in series and parallel 

94 Cavitation in centrifugal pumps 

9.5 Net positive suction head: definition. 

9,6 Net positive suction head: calculation 1 

97 Specific speed 

9.8 Net positive suction head: calculation 2 

99 Effect of reduced speed on pump characteristic 

9,10 Duty point and reduced speed of a centrifugal pump 
9.11 Power, impeller diameter, speed and delivered head 
9,12 Suction specific speed 

9.13 Reciprocating pumps 

9,14 Single-acting reciprocating piston 

9,15 Discharge from reciprocating pumps. 

9.16 Rotary pumps 

Further problems 


xiv 


193 
193 
196 
197 
199 
200 
202 
204 


209 
a2" 
214 
217 
219 
221 


226 
228 
230 


Glossary of terms 

Selected recommended texts 

Nomenclature and preferred units 

Useful conversion factors 

Physical properties of water (at atmospheric pressure) 
Losses for turbulent flow through fittings and valves 
Equivalent sand roughness of pipes 

‘Manning coefficient for various open channel surfaces 
Moody plot 

Index 


275 
282 
284 
285 
288 
289 
291 
292 
293 
294 


xv 








‘The scientist describes what is: 
the engineer creates what never was." 
Theodore von Kármán (1881-1963) 


‘I hear, and I forget 

I see, and I remember 
Ido, and I understand.’ 
Anonymous 








Fluid statics 


Introduction 


Fluids, whether moving or stationary, exert forces over a given area or surface. 





Fluids which are stationary, and therefore have no velocity gradient, exert 
normal or pressure forces whereas moving fluids exert shearing forces on the 
surfaces with which they are in contact. It was the Greek thinker Archimedes 
(c287BC-c212BC) who first published a treatise on floating bodies and 
provided a significant understanding of fluid statics and buoyancy. It was not 
for another 18 centuries that the Flemish engineer Simon Stevin (1548-1620) 
correctly provided an explanation of the basic principles of fluid statics, Blaise 
Pascal (1623-1662), the French mathematician, physicist and theologian, 
performed many experiments on fluids and was able to illustrate the funda- 
mental relationships involved. 

In the internationally accepted SI system (Systeme International d'Unités), 
the preferred derived units of pressure are Newtons per square metre (Nm?) 
with base units of kgm-!s-!. These units, also known as the Pascal (Pa), are 





relatively small. The term is therefore frequently used to represent one 
hundred thousand Newtons per square metre (105 Nm"? or 0.1 MPa). Many 
pressure gauges encountered in the process industries are still to be found cali- 
brated in traditional systems of units including the Metric System, the Absolute 
English System and the Engineers’ English System. This can lead to confusion 
in conversion although many gauges are manufactured with several scales. 
Further complication arises since the Pascal is a relatively small term and ST 
recommends that any numerical prefix should appear in the numerator of an 
expression. Although numerically the same, Nmm-? is often wrongly used 
instead of MNm-. 

It is important to note that the pressure of a fluid is expressed in one of two 
ways. Absolute pressure refers to the pressure above total vacuum whereas 











gauge pressure refers to the pressure above atmospheric, which itself is a var 
able quantity and depends on the local meteorological conditions. The atmo- 





spheric pressure used as standard corresponds to 101.3 kNm~? and is equivalent 
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to approximately 14.7 pounds force per square inch, or a barometric reading of 
760 mmHg. The pressure in a vacuum, known as absolute zero, therefore corre- 
sponds to a gauge pressure of 101.3 kNm? assuming standard atmospheric 
pressure. A negative gauge pressure thus refers to a pressure below atmospheric. 

The barometer is a simple instrument for accurately measuring the atmo- 
spheric pressure. In its simplest form it consists of a sealed glass tube filled with 
a liquid (usually mercury) and inverted in a reservoir of the same liquid. The 
atmospheric pressure is therefore exerted downwards on the reservoir of liquid 
such that the liquid in the tube reaches an equilibrium elevation. Above the 
liquid meniscus exists a vacuum, although in actual fact it corresponds to the 
vapour pressure of the liquid, In the case of mercury this is a pressure of 10 
kNm~ at 20°C. 

In addition to gauges that measure absolute pressure, there are many devices 
and instruments that measure the difference in pressure between two parts in a 
system. Differential pressure is of particular use for determining indirectly the 
rate of flow of a process fluid in a pipe or duct, or to assess the status of a 
particular piece of process equipment during operation — for example, identi- 
fying the accumulation of deposits restricting flow, which is important in the 
case of heat exchangers and process ventilation filters. 

Although there are many sophisticated pressure-measuring devices avail 
able, manometers are still commonly used for measuring the pressure in vessels 
or in pipelines, Various forms of manometer have been designed and generally 
are either open (piezometer) or closed (differential). For manometer tubes with 
a bore of less than 12 mm, capillary action is significant and may appreciably 
raise or depress the meniscus, depending on the manometric fluid. 

Finally, while fluids may be described as substances which offer no resis- 
tance to shear and include both gases and liquids, gases differ from liquids in 
that they are compressible and may be described by s mple gas laws. Liquids 
are effectively incompressible and for most practical purposes their density 
remains constant and does not vary with depth (hydrostatic pressure). At ultra 
high pressures this is not strictly true. Water, for example, has a 3.3% 
compressibility at pressures of 69 MNm-? which is equivalent to a depth of 7 
km. It was Archimedes who first performed experiments on the density of 
solids by immersing objects in fluids, The famous story is told of Archimedes 
being asked by King Hiero to determine whether a crown was pure gold 
throughout or contained a cheap alloy, without damaging the crown. Suppos- 
edly, while in a public bath, Archimedes is said to have had a sudden thought of 
immersing the crown in water and checking its density. He was so excited that 
he ran home through the streets naked shouting ‘Eureka! Eureka! — I have 
found it! I have found itl’. 
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1.1 Pressure at a point S 
Determine the total force on a wall of an open tank 2 m wide containing fuel 
oil of density 924 kgm™ at a depth of 2 m. 





Solution 

To determine the pressure at a point in the static liquid below the free surface, 
consider the equilibrium forces on a wedge-shaped element of the liquid. 
Resolving in the x-direction 


pAyLsin 0 - p,AyAz =0 


where 


sino - dt 
L 


Then 
Fem 
Resolving in the z-direction 
F * pAyLcos 0 — p; Axhy =0 


where 
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coso = À* 
L 


and the weight (force due to gravity) of the element is 
F = mg 


2 Ax: 
P735 





Ayg 


If the element is reduced to zero size, in the limit this term disappears because it 
represents an infinitesimal higher order than the other terms and may be 
ignored. Thus 


P=P2 


Note that the angle of the wedge-shaped element is arbitrary. The pressure p is 
therefore independent of 0. Thus, the pressure at a point in the liquid is the same 
in all directions (Pascal's law). To determine the pressure at a depth H, the equi- 
librium (upward and downward) forces are 


Pam Axdy + pAxAyHg. — pAxAy =0 
which reduces to 
P =Pam + PgH 
The pressure (above atmospheric) at the base of the tank is therefore 
p =pgH 
=924xg x2 
=18.129 x 10? Nm? 
The total force exerted over the wall is therefore 


p.P 
2 


_18129x10° x2x2 
2 


= 36258 x 10° N 


The total force is found to be 36.26 kKNm™?. 
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1.2 Pressure within a closed vessel 

A cylindrical vessel with hemispherical ends is vertically mounted on its axis. 
The vessel contains water of density 1000 kgm- and the head space is pressur- 
ized to a gauge pressure of 50 kNm~. The vertical wall section of the vessel has 
a height of 3 m and the hemispherical ends have radii of 1 m. If the vessel is 
filled to half capacity, determine the total force tending to lift the top dome and 
the absolute pressure at the bottom of the vessel. 


Hz25m 








Solution 
The total vertical force, F, tending to lift the dome is the pressure applied over 


the horizontal projected area 
F =p, nr? 

where py is the gauge pressure within the vessel. That is 
F 250x10) x nx P 


=156 x 10° N 
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Note that above the liquid surface the pressure in the head space is exerted 
uniformly on the inner surface of the vessel. Below the liquid, however, the 
pressure on the vessel surface varies with depth. The absolute pressure (pres- 
sure above a vacuum) at the bottom of the vessel is therefore 


P= Pam + Py + PH 
=101.3x 10? +50x 10° +1000x g x 25 


= 175.3 x 10° Nm 
The force tending to lift the dome is 1.56 MN and the pressure at the bottom of 
the vessel is 175.3 kNm7?. 

Note that, unlike the gas pressure which is exerted uniformly in the head 
space, the analysis to determine the hydrostatic forces acting on the submerged 
curved surface (lower domed section) requires resolving forces in both the 
vertical and horizontal directions. The magnitude of the horizontal reaction on 
the curved surface is equal to the hydrostatic force which acts on a vertical 
projection of the curved surface, while the magnitude of the vertical reaction is 
equal to the sum of the vertical forces above the curved surface and includes the 
weight of the liquid. In this case, however, the vessel is symmetrical such that 
the hydrostatic force is in the downward direction. The downward force 
imposed by the gas and liquid is thus 


2n? 





F=(py + pgh nr? +pg 


2 3 
= (50,000 + 1000 x g x 15)x ex 12 +1000% g x EXE 


= 223,854 N 
= 223.8 KN 








FLUID STATICS 


1.3 Forces within a hydraulic ram 

A hydraulic ram consists of a weightless plunger of cross-sectional area 
0.003 m? and a piston of mass 1000 kg and cross-sectional area 0.3 m?. The 
system is filled with oil of density 750 kgm™. Determine the force on the 
plunger required for equilibrium if the plunger isat an elevation of 2 m above 
the piston. 





Solution 
For the piston, the pressure at the datum elevation xx is 


eti 
a 


Px 


where F} is the force of the piston and 4; is the area of the piston. This pressure 
is equal to the pressure applied by the plunger at the same datum elevation. That 
is 


F, 
Pis = + Po 
a 


where F3 is the force on the plunger, az is the area of the plunger and H is the 
elevation of the plunger above the datum. Therefore 
in 


F; 
= + Pog 
a a 
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Rearranging 
F| 
Fy E E 
a 


«ems, (eror 


-750x x2] 


=54N 


The force required for equilibrium is found to be 54 N. Note that if no down- 
ward force is applied to the weightless plunger, the plunger would rise to an 
elevation of 4.44 m. 

The hydraulic ram illustrated is an example of a closed system in which the 
pressure applied by the piston is transmitted throughout the hydraulic fluid 
(oil). The principle of pressure transmission is known as Pascal's law after 
Pascal who first stated it in 1653. Hydraulic systems such as rams, lifts and 
jacks are based on this principle and are useful for lifting and moving purposes. 
It is usual in such hydraulic systems to replace the piston with compressed air. 
The force applied is then controlled by the applied air pressure. High pressures 
can therefore be achieved, as in the case of hydraulic presses, in which the force 
exerted against à piston in turn exerts the force over a smaller area. For 
example, the plunger shown corresponds to a diameter of 62 mm over which an 
equilibrium pressure of 18 kNm-? is applied. If it were to be connected to a 
shaft 18 mm in diameter, then the force exerted over the area of the shaft would 
correspond to 222 kNm-? — a factor of 12 times greater. 
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14 Liquid-liquid interface position in a 

solvent separator 

Mixtures which contain two mutually insoluble organic and aqueous liquids 
are to be separated in a separator which consists of a vertical chamber with 
overflow and underflow. The mixture is, fed slowly to the separator in which the 
aqueous phase, of constant density 1100 kgm, is discharged from the under- 
flow at the base of the chamber to a discharge point 50 cm below the overflow 
level in the chamber. The organic phase can vary in density from 600-800 
kgm. Determine the minimum height of the chamber, H, which can be used if 
the organic phase is not to leave with the aqueous phase. If the height H is made 
equal to 3 m, determine the lowest possible position of the interface in the 
chamber below the overflow. 









Vent to 
atmosphere. 


E — [50m 
Overflow Organic EN I 


Solution : 
The separator is assumed to operate at atmospheric pressure. Equating the pres- 
sure in the chamber and discharge point for the maximum possible depth (in 
metres) for the organic phase in the chamber gives 
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Po8H =Paqg(H —05) 
where p, and P ag are the densities of organic and aqueous solutions, respec- 
tively. Rearranging 
OSP aq 


Hz——. 


Pag 7 Po 


To ensure no loss of organic phase with the aqueous phase, the height of the 
chamber is greatest for the highest possible organic density (800 kgm-3). 
Therefore 
dn 05 x 1100 
1100 — 800 


=22m 


For a fixed length of chamber of 3 m, the interface between the two pha: 
determined from the pressure in the chamber and discharge point. That is 





Po8H, * pau gH» — pau g(H —05) 
where 

H =H; +H, 
Therefore 

PogHi +Pagg(H -Hy)-pqug(H -05) 


Rearranging, the interface position is at its lowest position when the organic 
phase has a density of 600 kgm-?. That is 


05p ag 


H, =— 
Paq “Po 

= 05x1100 
1100 -600 


=llm 


The maximum depth is found to be 2.2 m and the interface below the overflow 

is found to be 1.1 m. Ideally, the feed point to the chamber should be located at 

the liquid-liquid interface to ensure quick and undisturbed separation. Where 

the density of the organic phase is expected to vary, either an average position 

i the position corresponding to the most frequently encountered density may 
used. 
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1.5 Liquid-liquid interface measurement by 

differential pressure 

Aqueous nitric acid is separated from an insoluble oil in a vessel. Dip legs 
extend into both phases through which air is gently discharged sufficient to 
overcome the hydrostatic pressure. Determine the position of the interface 
berween the legs if the legs are separated a distance of 1 m for which the differ- 
ential pressure between the legs is 10 kNm"7. The densities of oil and nitric acid 
are 900 kem? and 1070 kgmr^2, respectively. 





Nitric acid 


Solution 

The use of dip legs is an effective way of measuring liquid densities, liquid- 
liquid interface positions and detecting the presence of solid material in liquids. 
As it has no moving or mechanical parts it is essentially maintenance free and it 
has therefore found application in the nuclear industry amongst others. In this 
application, the dip legs are used to determine the position of the liquid-liquid 
interface in which the densities of the two phases are assumed to be constant. 
The differential pressure between the legs is 


Ap zpogH, + Pn SH 
where p, andp,, are the densities of the oil and nitric acid and where the fixed 
distance between the ends of the dip legs is 


H=H,+H, 


=Im 
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Eliminating H3 and rearranging 


_ AP —PngH 
(po —Pnd8 


Ay 





The depth is found to be 30 cm below the upper dip leg. 

Note that a single dip leg can be used to determine the depth of liquids of 
constant density in vessels in which the gas pressure applied is used to over- 
come the hydrostatic pressure. For cases in which the density of the liquid is 
likely to vary, due to changes in concentration or the presence of suspended 
solids, the density can be determined using two dip legs of different length, the 
ends of which are a fixed distance apart. In the more complicated case of two 
immiscible liquids in which the densities of both phases may vary appreciably, 
it is possible to determine the density of both phases and the location of the 
interface using four dip legs with two in each phase. 

In practice, itis necessary to adjust carefully the gas pressure until the hydro- 
static pressure is just overcome and gas flows freely from the end of the dip 
leg(s). Sensitive pressure sensing devices are therefore required for the low 
gauge pressures involved. Fluctuating pressure readings are usually experi- 
enced, however, as the gas bubbles form and break off the end of the leg. 
Conversion charts may then be used to convert a mean pressure reading to 
concentration, interface position or liquid volume, as appropriate. 
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1.6 Measurement of crystal concentration by 
differential pressure 

The concentration of sodium sulphate crystals in a liquid feed to a heat 
exchanger is determined by a differential pressure measurement of the satu- 
rated liquid in the vertical leg feeding the heat exchanger. If the pressure 
measurements are separated by a vertical distance of 1.5 m, determine the 
density of the solution with crystals and the fraction of crystals for a differential 
pressure measurement of 22 kNm~. The density of saturated sodium sulphate 
is 1270 kem? and density of anhydrous solution sulphate is 2698 kgm3. 











Heat 
exchanger 





Condensate 
He 














Solution 
Assuming no differential pressure loss due to friction in the leg, the differential 
pressure is due to the static pressure between the pressure measurement points. 
That is 

Ap =pgH 


where p is the density of the solution. 
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Rearranging 
Ap 


poH 


_ 22x10 
g x15 
=1495 kgm? 


The density of the solution with crystals is 1495 kgm™. This density is greater 
than that of the saturated sodium sulphate solution alone and therefore indicates 
the presence of crystals for which the fractional content is found from 


P=SfiPs +frPc 


wherep , is the density of saturated solution, p , is the density of crystals, and f, 
and f» are the respective fractions where 


Sith =i 
Eliminating fı 
hs P-P; 
Pe -Ps 
1495 -1270 
2698 — 1270 
=0157 


That is, the crystal content is found to be 15.7%. This is, however, an overesti- 
mate since frictional effects of the flowing liquid in the leg are ignored. Where 
they can not be ignored the differential pressure is modified to 


Ap -pg(H - Hj) 


where H ;, is the head loss due to friction. 
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1.7 Pressure within a gas bubble 

A small gas bubble rising in an open batch fermenter has a radius of 0.05 cm 
when it is 3 m below the surface. Determine the radius of the bubble when it is 1 
m below the surface. It may be assumed that the pressure inside the bubble is 
2 or above the pressure outside the bubble, where r is the radius of the bubble 
and Gis the surface tension of the gas-fermentation broth and has a value of 
0.073 Nm. The pressure and volume of the gas in the bubble are related by the 
expression pV =c where c is a constant. 


Fd d deese 





Solution 


At a depth of 3 m, the pressure within the bubble, p}, is dependent on the pres- 
sure at the free surface, the hydrostatic pressure and surface tension effect. Thus 


2 
Pi =Pam +p8H; +E 
n 


2x0073 


-1013x10? 41000x g x 3 
5x107* 


7131022 x 10? Nm? 
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At a depth of 1 m, the pressure inside the bubble, p>. is 
26 
P2 = Parm *pHa 


2x 0073 
pea 
n 


—101.3 x 10? 1000 x g x1 


0146 
n 


=11111x10° + 


Since pV is a constant, then 


PAV) =P2V2 
where for a spherical bubble 

nite? pte! 
That is 

pit =p} 
Therefore 


0146 
131022 x 10? x (5x 1075 [ns 105 + rd. 5 


"n 
The cubic equation can be solved analytically, by trial and error or by assuming 
that the second term in the brackets is substantially small, reducing the effort 
required for solution to yield a bubble radius of approximately 0.053 mm. 
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1.8 Pressure measurement by differential manometer 
Determine the pressure difference between two tapping points on a pipe 
carrying water for a differential manometer reading of 20 cm of mercury. The 
specific gravity of mercury is 13.6. 





Solution 

The differential or U-tube manometer is a device used to measure the difference 
in pressures between two points and consists of a transparent U-tube, usually 
made of glass, and contains a manometric fluid such as mercury. It is typically 
used to measure the pressure drop of moving fluids due to friction along pipes 
or due to obstacles in pipelines such as flow measuring devices, fittings and 
changes in geometry. The pressure difference of the process fluid is indicated 
by the difference in levels of the manometric fluid between the two vertical legs 
of the U-tube which, at the datum elevation xx, are 


py * pg(H, * H)- py * psHy * py, 8H 


where P yg is the density of mercury and p the density of water. 
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Rearranging, the differential pressure Ap between the legs is 
Ap = Pi —P2 
= pg t pu, gH —pg(H, * H) 
=(PHg -P)8H 
= (13,600 — 1000) x g x 02 


224721 x 10? Nm? 


The differential pressure is 24.7 kNm"?. 

Note that the location of the manometer below the pipe, H;, is not required in 
the calculation. In practice it is important to allow sufficient length in the legs to 
prevent the manometric fluid reaching the tapping point on the pipe for high 
differential pressures. Filled with mercury, differential manometers can typi- 
cally be used to measure differentials up to about 200 kNm^? or with water to 
about 20 kNm"?. Where a temperature variation in the process fluid is 
expected, it is important to allow for density-temperature variation of the mano- 
metric fluid, which can affect readings. 

In general, the U-tube differential manometer as a pressure-measuring 
device is largely obsolete. There are many sophisticated methods and pressure- 
measuring devices now used by industry. But the differential manometer 
continues to be a useful tool in the laboratory and for testing purposes. 
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1.9 Pressure measurement by inverted manometer 

A laboratory rig is used to examine the frictional losses in small pipes. Deter- 
mine the pressure drop in a pipe carrying water if a differential head of 40 cm is 
recorded using an inverted manometer. 








Solution 

The inverted manometer avoids the use of a manometric fluid and instead uses 
the process fluid (water in this case) to measure its own pressure. It consists of 
an inverted U-tube with a valve into which air or an inert gas can be added or 


vented. Here, the pressure at the datum elevation xx, in left and right hand legs 
is 


Dj -pg(H Hy) py —pgH, —paj, gH 


where p is the density of water and p ,;, is the density of air. Rearranging, the 
differential pressure Ap is therefore 


Ap = Pi ~P2 
= pg (H Hi) - psH, - pai, gH 


=(P -Pair )8H 


19 














FLUID MECHANICS 


Since the density of air is in the order of 1000 times less than that of water, it 
may therefore be reasonably assumed that the differential pressure is approxi- 
mated to 


Ap = pg 
=1000x g x 04 


-394Nm^? 


The differential pressure is found to be 3.9 kNm~. As with the differential 
manometer, the elevation of the manometer, Hj, is not required in the calcula- 
tion. In practice, however, it is important to ensure a reasonable position of 
liquid levels in the legs. This is best achieved by pressurizing the manometer 
with air or inert gas using the valve, where for high pressures the density may 
become appreciable and should be taken into consideration. In the case of air, 
the error in the calculation is unlikely to be greater than 0.5%. In the case illus- 
trated, the density of water corresponds to a temperature of 10°C for which the 
density of air at atmospheric pressure is 1.2 kgm”, If this had been taken into 
account, it would have yielded a differential pressure of 3919 Nm~ or an error 
of 0.12%. A more significant error is likely to be due to the effects of tempera- 
ture on density and may affect the result by as much as 1%, Other errors are 
likely to be caused by defining the top level of the manometric fluid in the 
vertical leg due to its meniscus. A column-height accuracy of 0.025 mm is, 
however, generally achievable with the keenest eye reading. 
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1.10 Pressure measurement by single leg manometer 

A mercury-filled single leg manometer is used to measure the pressure drop 
across a section of plant containing a process fluid of density 700 kgm~. The 
pressure drop is maintained by an electrical device which works on an on/off 
principle using a contact arrangement in a narrow vertical tube of diameter 2 
mm while the sump has a diameter of 2 cm. If the pressure drop across the plant 
is to be increased by 20 kNnr?, determine the quantity of mercury to be 
removed from the sump if the position of the electrical contact cannot be 
altered. 


Plant 


roce € 


Signal to 
pressure 
control 
mechanism. 







Sensing 
device 


| Electrical 
contact 


Solution 

The single leg manometer uses a sump or reservoir of large cross-section in 
place of one leg. When a differential pressure is applied, the level in the leg or 
tube rises due to a displacement from the sump. The ratio of leg to sump area is 
generally needed for particularly accurate work but is ignored for most 
purposes since the area of the sump is comparatively larger than that of the leg. 
The device in this case operates when the level of mercury in the tube falls, 
breaking the electrical circuit. The pressure control mechanism therefore 
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receives a signal to increase the pressure difference. When the mercury level 
rises, the opposite occurs. An increase in pressure drop of 20 kNm^ therefore 
corresponds to an increase in difference in level of mercury of 


E- 
(pu, -P)s 


| 20x10? 
(13,600 — 700) x g 


=0.158 m 


The volume of mercury to be removed to ensure the contact is still just made is 
therefore 


2 
yay 
4 
2 
= mx 0.02" , 158 
= 496 x 107 m> 


That is, the volume to be removed is approximately 50 ml. Note that if the 
displacement of mercury from the sump into the tube is taken into account then 
this would correspond to a drop in level in the sump, H ,, of 


a 
Hy =H 


= (20002), o1ss 
002 
2158x107 m 


This is very small and ignoring it is justified. 
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1.11 Pressure measurement by inclined leg manometer 
An oil-filled inclined leg manometer is used to measure small pressure changes 
across an air filter in a process vent pipe. If the oil travels a distance of 12 cm 
along the leg which is inclined at an angle of 20? to the horizontal, determine 
the gauge pressure across the filter. The density of oil is 800 kgmr. 


Pi 





Original level -- - 


Solution 

This instrument is useful for measuring small differential pressures and 
consists of a sump of manometric fluid (oil) with a leg extended down into it 
and inclined at some small angle. Applying a differential pressure across the 
sump and the leg results in a displacement of the manometric fluid into the leg, 
the distance the manometric liquid travels up along the leg being a measure of 
differential pressure and is 


Ap =p1 -P2 
=pg(H; +H) 

If the oil is displaced from the sump up along the leg by a distance L, the corre- 
sponding drop in level in the sump, H ;, is therefore 

a, = 

A 

Also, the vertical rise of the oil is related to length by the sine of the angle of the 
inclined leg. That is 

H3 - Lsin8 
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The differential pressure is therefore 


ap =pa( SE» tsino) 
a ` 
=pgL| — + sin 0 

wpa 


As no details are provided regarding the dimensions of the manometer, the 
cross-sectional area of the oil sump, A, is therefore assumed to be very much 
larger than the area of the leg, a. The equation therefore reduces to 


Ap = pgLsin @ 
=800x g x 012x sin 20° 


-32 Nm? 


The differential pressure is found to be 322 Nm. 

The device is particularly useful for measuring small differential pressures 
since if the terms inside the brackets are kept small it allows the length along the 
inclined leg, L, to be appreciable. If, for a given differential pressure, the equiv- 
alent movement of manometric liquid up a vertical leg would be A, say, then the 
ratio of movements L to h 

LUN V 
ho +sin@ 

A 


can therefore be considered as a magnification ratio. 
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1.12 Archimedes’ principle 

A vessel containing a process material with a combined total mass of 100 kg is 
immersed in water for cooling purposes. Determine the tension in the cable of 
an overhead crane used to manoeuvre the fully immersed container into its 
storage position if the bulk density of the vessel is 7930 kgm. 





Container 


Solution 
Consider a body of mass m, immersed in the liquid such that the net downward 
force is the difference between the downward and upward forces, That is 


F=m,8 -mg 


where m is the mass of water displaced. This is known as Archimedes’ prin- 
ciple and states that when a body is partially or totally immersed, there is an 
upthrust equal to the weight of fluid displaced. For the immersed object, the net 
downward force is taken by the tension in the cable and can be determined 
where the mass of the container and water displaced is related to volume by 
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where p, is the bulk density of the container and p is the density of water. 
Rearranging, the mass of water displaced by the container is therefore 





The tension in the cable is therefore 


F = meg] pP 
Pe 
i] 


-100x g x|1- 25 
7930 


=857N 


That is, the tension in the cable is 857 N. Note that the tension in the cable when 
the vessel is lifted out of the water is 


f-2meg 
=100x g 
-981N 


The buoyancy effect therefore reduces the tension in the cable by 124 N. 
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1.13 Specific gravity measurement by hydrometer 

A hydrometer floats in water with 6 cm of its graduated stem unimmersed, and 
in oil of SG 0.8 with 4 cm of the stem unimmersed. Determine the length of stem 
unimmersed when the hydrometer is placed in a liquid of SG 0.9. 





Solution 

Hydrometers are simple devices for measuring the density or specific gravity of 
liquids routinely used in the brewing industry to determine quickly the conver- 
sion of sugar to alcohol in fermentation. They consist of a glass tube which have 
a weighted glass bulb and graduated stem of uniform diameter and float in the 
liquid being tested. The density or specific gravity (SG) is usually read directly 
from the graduated stem at the depth to which it sinks. For no net downward 
force. the vertical downward forces acting on the body are equal to the upthrust. 
Thus 


mg -m,yg 
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where m and m, are the mass of liquid and hydrometer, respectively. Thus, 
Archimedes’ principle for a floating body states that when a body floats, it 
displaces a weight of fluid equal to its own weight. The displacement by the 
hydrometer is therefore 


mpg =pg(L-x)a 
=Po8(L -xo)a 


where a is the cross-sectional area of the stem, p and p,, are the densities of 
water and oil, and x and x,, are the lengths of stem unimmersed in the respective 
liquids. Therefore 


pe(L — 5)a - pua — x) 


Rearranging, the length of hydrometer is therefore 


L= PYT Poto 
Po =P 
.. 1000 x 006 — 800 x 004 
E 800 — 1000 
=014m 


For the hydrometer immersed in a liquid of SG 0.9 (900 kgm), let the length 
of stem remaining unimmersed be x; . Therefore 


1000 x g x (014 —006)x a 2900 x g x (014 —x;)x a 


Solving, x, is found to be 0.0511 m. That is, the length of stem above the liquid 
of SG 0.9 is 5.11 cm. 
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1.14 Transfer of process liquid to a ship 

A liquid hydrocarbon mixture of density 950 kgm is transferred by pipeline to 
a ship at a loading terminal. Prior to transfer, the ship has an unloaded 
displacement of 5000 tonnes and draft of 3 m. Transfer of the hydrocarbon is at 
a steady rate of 125 nh. If the sea bed is at a depth of 5.5 m, determine the 
quantity delivered and time taken if the ship requires at least 1 m of clearance 
between the sea bed and hull to manoeuvre away safely from the loading 
terminal. 

















Solution 
Applying Archimedes' principle, the ship prior to transfer displaces its own 
weight of sea water. That is 


msg = mg 
= pATig 
where m, and m are the mass of ship and sea water displaced. p is the density of 


sea water, A is the water plane area and T} is the depth of the ship below the 
waterline. After transfer 


msg + Mpeg -pAT)g 


where m),. is the mass of hydrocarbon mixture. After transfer, the ship is clear 
from the sea bed by 1 m. Combining these two equations, the mass of hydro- 
carbon transferred is 
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n 
my —m,—-l 


«sae x(5-1) 


-225x10f kg 
The transfer time is therefore 
Mhe 
PhO 
_ 25x10° 
950 x 125 
= 2105 h 


t= 


That is, a transfer of 2500 tonnes of hydrocarbon mixture is completed in 21.05 
hours. 

It should be noted that the approach illustrated is rather simplistic. No 
account is made for the dimensions of the ship in terms of its length and beam 
nor the variation of the water plane area with depth. The beam is an important 
dimension in terms of stability where the stability is dependent on the relative 
position of the ship’s centre of gravity and centroid of the displaced volume 
called the centre of buoyancy. A ship is unstable and will capsize when, for a 
heel of up to 10°, a line drawn vertically up from the centre of buoyancy is 
below the centre of gravity — a point known as the metacentre. 

The safe transfer of liquids to and from tanks within ships requires a careful 
sequence of operation. Tidal effects on moored ships and the effects of the 
liquid free surface in the tanks must also be taken into consideration. It was the 
British politician Samuel Plimsoll (1824-1898) who was responsible for 
getting legislation passed to prohibit ‘coffinships’ — unseaworthy and over- 
loaded ships — being sent to sea. The Merchant Sea Act of 1874 included, 
amongst other things, enforcement of the painting of lines, originally called 
Plimsoll marks and now known as load line marks, to indicate the maximum 
load line which allows for the different densities of the world's seas in summer 
and winter. 
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Further problems 


(1) Explain what is meant by gauge pressure and absolute pressure. 


(2) A hydraulic press has a ram of 10 cm diameter and a plunger of 1 cm 
diameter. Determine the force required on the plunger to raise a mass of 500 kg 
on the ram. 


Answer: 49.05 N 
(3) The reading of a barometer is 75.5 cm of mercury. If the specific gravity 
of mercury is 13.6, convert this pressure to Newtons per square metre. 
Answer: 100,792 Nm? 
(4) A rectangular tank 5 m long by 2 m wide contains water to a depth of 2 m. 


Determine the intensity of pressure on the base of the tank and the total pressure 
on the end, 


Answer: 19.6 kNm-?, 39 kNm-? 
(5) Determine the total pressure on a vertical square sluice, of 1 m square, 
positioned with its top edge 3 m below the level of water, 

Answer: 34.3 kNm? 
(6) A tube is filled with water to a depth of 600 mm and then 450 mm of oil of 


SG 0.75 is added and allowed to come to rest. Determine the gauge pressure at 
the common liquid surface and at the base of the tube. 


Answer: 3.3 kNm"?, 9.2 kNm^? 


(7) Show that when a body is partially or totally immersed in a liquid, there is 
an upthrust on the body equal to the weight of the liquid displaced. 


(8) Show that a floating body displaces a weight of the liquid equal to its own 
weight. 


(9) A U-tube has a left-hand leg with a diameter of 5 cm and a right-hand leg 
with a diameter of | cm and inclined at an angle of 24°. If the manometer fluid is 
oil with a density of 920 kgm-? and a pressure of 400 Nm is applied to the 
left-hand leg, determine the length by which the oil will have moved along the 
right-hand leg. 


Answer: 9.9 cm 
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(10) Determine the absolute pressure in an open tank containing crude oil of 
density 900 kgm^? at a depth of 5 m. 


Answer: 145.4 kNm-? 


(11) An open storage tank 3 m high contains acetic acid, of density 1060 
kgm, and is filled to half capacity. Determine the absolute pressure at the 
bottom of the tank if the vapour space above the acid is maintained at atmo- 
spheric pressure, 


Answer: 117 kNm? 


(12). A differential manometer containing mercury of SG 13.6 and water indi- 
cates a head difference of 30 cm. Determine the pressure difference across the 
legs. 


Answer: 37.1 kNm? 


(13) A U-tube contains water and oil. The oil, of density 800 kgm™?, rests on 
the surface of the water in the right-hand leg to a depth of 5 cm. If the level of 
water in the left-hand leg is 10 cm above the level of water in the right-hand leg, 
determine the pressure difference between the two legs. The density of water is 
1000 kgm-3. 


Answer: 589 Nm? 


(14). A separator receives continuously an immiscible mixture of solvent and 
aqueous liquids which is allowed to settle into separate layers, The separator 
Operates with a constant depth of 2.15 m by way of an overflow and underflow 
arrangement from both layers. The position of the liquid-liquid interface is 
monitored using a dip leg through which air is gently bubbled. Determine the 
position of the interface below the surface for a gauge pressure in the dip leg of 
20 kNm. The densities of the solvent and aqueous phases are 865 kgm™ and 
1050 kgm-^, respectively, and the dip leg protrudes to within 5 cm of the 
bottom of the separator. 


Answer: 90 cm 
(15) A hydrometer with a mass of 27 g has a bulb of diameter 2 cm and length 


8 cm, and a stem of diameter 0.5 cm and length 15 cm. Determine the specific 
gravity of a liquid if the hydrometer floats with 5 cm of the stem immersed. 


Answer: 1.034 
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(16) Two pressure tapping points, separated by a vertical distance of 12.7 m, 
are used to measure the crystal content of a solution of sodium sulphate in an 
evaporator. Determine the density of the solution containing 25% crystals by 
volume and the differential pressure if the density of the anhydrous sodium 
sulphate is 2698 kgm and the density of saturated sodium sulphate solution is 
1270 kgm^. 

Answer: 1627 kem-?, 203 kNm? 


(17) A vacuum gauge consists of a U-tube containing mercury open to atmos- 
phere. Determine the absolute pressure in the apparatus to which it is attached 
when the difference in levels of mercury is 60 cm. 


Answer: 21.3 kNm? 


(18) Determine the height through which water is elevated by capillarity in a 
glass tube of internal diameter 3 mm if the hydrostatic pressure is equal to 40/d 
where o is the surface tension (0.073 Nm~!) and d is the diameter of the tube. 


Answer: 9.9 mm 


(19) Explain the effect of surface tension on the readings of gauges of small 
bore such as piezometer tubes. 


(20) A ship has a displacement of 3000 tonnes in sea water, Determine the 
volume of the ship below the water line if the density of sea water is 1021 
kgm-3, 


Answer: 2938 m? 


Q1) A closed cylindrical steel drum of side length 2 m, outer diameter 1.5 m 
and wall thickness 8 mm is immersed in a jacket containing water at 20°C 
(density 998 kgm). Determine the net downward and upward forces when the 
drum is both full of water at 20°C and empty. The density of steel is 7980 
kgm", 

Answer: 5.17 kN, -29.4 kN 


(22) Anoil/water separator contains water of density 998 kgm:^? to a depth of 
75 cm above which is oil of density 875 kgm™ to a depth of 75 cm. Determine 
the total force on the vertical side of the separator if it has a square section 1.5 m 
broad. If the separator is pressurized by air above the oil, explain how this will 
affect the answer. 


Answer: 16 kN 
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Continuity, 
momentum 
and energy 


Introduction 

With regard to fluids in motion, it is convenient to consider initially an ideal- 
ized form of fluid flow. In assuming the fluid has no viscosity, it is also deemed 
to have no frictional resistance either within the fluid or between the fluid and 
pipe walls. Inviscid fluids in motion therefore do not support shear stresses 
although normal pressure forces still apply. 

There are three basic conservation concepts evoked in solving problems 
involving fluids in motion. The conservation of mass was first considered by 
Leonardo da Vinci (1452-1519) in 1502 with respect to the flow within a river. 
Applied to the flow through a pipe the basic premise is that mass is conserved. 
Assuming no loss from or accumulation within the pipe, the flow into the pipe 
is equal to the flow out and can be proved mathematically by applying a m: 
balance over the pipe section. The flow of incompressible fluids at a steady r: 
is therefore the simplest form of the continuity equation and may be readily 
applied to liquids. 

The conservation of momentum is Newton’s second law applied to fluids in 
motion, and was first considered by the Swiss mathematician Leonhard Euler 
(1707-1783) in 1750. Again, by considering inviscid fluid flow under steady 
flow conditions, calculations are greatly simplified. This approach is often 
adequate for most engineering purposes. 

The conservation of energy was first considered by the S 
Daniel Bernoulli (1700-1782) in 1738 to describe the conservation of 
mechanical energy of a moving fluid in a system. The basic premise is that the 
total energy of the fluid flowing in a pipe must be conserved. An energy balance 
on the moving fluid across the pipe takes into account the reversible pressure- 
volume, kinetic and potential energy forms, and is greatly simplified by consid- 
ering steady, inviscid and incompressible fluid flow. 
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2.1 Flow in branched pipes 

Water flows through a pipe section with an inside diameter of 150 mm at a rate 
of 0.02 m?s-!. The pipe branches into two smaller diameter pipes, one with an 
inside diameter of 50 mm and the other with an inside diameter of 100 mm. If 
the average velocity in the 50 mm pipe is 3 ms^!, determine the velocities and 
flows in all three pipe sections. 


Pipe2  d,- 50mm 





Solution 

The continuity equation is effectively a mathematical statement describing the 
conservation of mass of a flowing fluid where the mass flow into a pipe section 
is equal to the mass flow out. That is 


Pia), 7 p202v2 
For an incompressible fluid in which the density does not change, the volu- 
metric flow is therefore 

avi =a 


For the branched pipe system in which there is no loss or accumulation of the 
incompressible process fluid (water), the flow through the 150 mm diameter pipe 
(Pipe 1) is equal to the sum of flows in the 50 mm (Pipe 2) and 100 mm diameter 
pipes (Pipe 3). That is 


Qi =Q +03 
nd? nd? 
EB aie iius 
4 * -4 


Rearranging, the velocity in the 100 mm diameter pipe is therefore 
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This corresponds to a flow of 


2 
3 

=—3y 

Q3 a 





_ «x0? 
4 


x18 


=0014m's"! 


Similarly, the velocity and flow can be found for the other two pipes and are 
given below. 


Pipe | Pipe 2 Pipe 3 
Diameter, mm 150 50 100 
Velocity, ms“! 1.13 3.00 1.80 
Flowrate, m3s~! 0.020 0.006 0.014 
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2.2 Forces on a U-bend 

A horizontal pipe has a 180* U-bend with a uniform inside diameter of 200 mm 
and carries a liquid petroleum fraction of density 900 kgm~ at a rate of 150 
mtr, Determine the force exerted by the liquid on the bend if the gauge pres- 
sure upstream and downstream of the bend are 100 kNm-? and 80 kNm-?, 
respectively. 


Solution 

The thrust exerted by the flowing liquid on the horizontal bend is resolved in 
both the x- and y-directions. Assuming that the gauge pressures of the liquid are 
distributed uniformly in the U-bend, then resolving the force in the x-direction 
gives 


F, — py cos 0, — pyay cos 8 - pQ(v; cos» —v, cos0j) 
and in the y-direction 

Fy 2 pyaj sin 8j paa» sin65 —pQ(vs sin 65 — v, sin0j) 
The respective upstream and downstream pressure forces are 


2 
pid 21x 105 x $027 


=3141N 
and 


2 
paa; - 8x10* x T502. 


=2513 N 
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For the uniform cross-section, the average velocity remains constant. That is 
x= wy 
EA 
nd? 
150 


A 
x 


—onx02) 
=1.33 ms! 
The momentum fluxes are therefore 
pQv, = pQv2 


= 900 x 159 x 1.33 
3600 


=499N 


For the liquid entering the 180° bend the angle 8, is 0° and for the liquid leaving 
5 is 180°, The resolved force in the x-direction is therefore 


F, — 3141 x cos0°-2513 x cos180°+ 499 x (cos 180°-cos0°) 
=5554.N 
Since sin 0° and sin 180° are equal to zero, the force in the y-direction is 
Fy =0 


Although not taken into consideration here, the reaction in the vertical direction 
F. can also be included where the downward forces are due to the weight of the 
bend and the fluid contained within it. 
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2.3 Pressure rise by valve closure 

A valve at the end of a water pipeline of 50 mm inside diameter and length 500 m 
is closed in 1 second giving rise to a uniform reduction in flow. Determine the 
average pressure rise at the valve if the average velocity of the water in the 
pipeline before valve closure had been 1.7 ms~!. 


Solution 

When a liquid flowing along a pipeline is suddenly brought to rest by the 
closure of a valve or any other obstruction, there will be a large rise in pressure 
due to the loss of momentum causing a pressure wave to be transmitted along 
the pipe. The corresponding force on the valve is therefore 


F2m* 
t 


where v/t is the deceleration of the liquid and the mass of water in the pipeline 
is 


m=paL 
Thus 
v 
F 2 paL— 
B t 


2 
= 1000x 005 5003 T 


=1669 N 


corresponding to a pressure on the valve of 


a|% 


p 


è 


F 
nd? 

_ 4x 1669 
mx 0.05? 





= 850.015 x 10° kNm? 


The average pressure on the valve on closure is found to be 850 kKNm™?. Serious 
and damaging effects due to sudden valve closure can occur, however, when 
the flow is retarded at such a rate that a pressure wave is transmitted back 
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along the pipeline. The maximum (or critical) time in which the water can be 
brought to rest producing a maximum or peak pressure is 


where c is the velocity of sound transmission through the water. With no resis- 
tance at the entrance to the pipeline, the excess pressure is relieved. The pres- 
sure wave then travels back along the pipeline reaching the closed valve at a 
time 21/c later. (The period of 2L/c is known as the pipe period.) In practice, 
closures below values of 2L/c are classed as instantaneous. In this problem, 
the critical time corresponds to 0.67 seconds for a transmission velocity of 1480 
ms-! and is below the 1.0 second given. The peak pressure can be significantly 
greater than the average pressure on valve closure with the pressure wave being 
transmitted up and down the pipeline until its energy is eventually dissipated. It 
is therefore important to design piping systems within acceptable design limits. 
Accumulators (air chambers or surge tanks) or pressure relief valves located 
near the valves can prevent potential problems. 

The peak pressure resulting from valve closures faster than the pipe period 
can be calculated (in head form) from 


Hz 


8 
This basic equation, developed by the Russian scientist N. Joukowsky in 1898, 
implies that a change in flow directly causes a change in pressure, and vice 
versa. The velocity of sound transmission, c, is however variable and is 
dependent upon the physical properties of the pipe and the liquid being 
conveyed. The presence of entrained gas bubbles markedly decreases the effec- 
tive velocity of sound in the liquid. In this case, the peak head is 


_ 17x 1480 
8 
= 2565 m 





H 


which corresponds to a peak pressure of 2516 kNm~?. However, the 
Joukowsky equation neglects to consider the possible rise due to the reduction 
in frictional pressure losses that occur as the fluid is brought to rest. It also does 
not consider the pressure in the liquid that may exist prior to valve closure — all 
of which may well be in excess of that which can be physically withstood by the 


pipe. 
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2.4 The Bernoulli equation 

An open tank of water has a pipeline of uniform diameter leading from it as 
shown below. Neglecting all frictional effects, determine the velocity of water 
in the pipe and the pressure at points A, B and C. 


Ls 
Free surface 





Solution 
The Bernoulli equation (named after Daniel Bernoulli) is 
2 
v v. 
Pi Suus s Pra as: 
P8 2g Pg 2 


The first, second and third terms of the equation are known as the pressure head, 
velocity head and static head terms respectively, each of which has the funda- 
mental dimensions of length. This is an important equation for the analysis of 
fluid flow in which thermodynamic occurrences are not important. It is derived 
for an incompressible fluid without viscosity. These assumptions give results of 
acceptable accuracy for liquids of low viscosity and for gases flowing at 
subsonic speeds when changes in pressure are small. 

To determine the velocity in the pipe, the Bernoulli equation is applied 
between the free surface (point 1) and the end of the pipe (point 2) which are 
both exposed to atmospheric pressure. That is 


Pi =P2 =Patm 


The tank is presumed to be of sufficient capacity that the velocity of the water at 
the free surface is negligible. That is 


v1 50 


Therefore 
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The average velocity is the same at all points along the pipeline. That is 
Va V4 vg yp 


The pressure at A is therefore 


2 
v» 
Pa =P8| zi -za -7-A 
A =P8| z1 -ZA 2g 
2 
=1000g x ieee 
2g 
= 17,658 Nm 
The pressure at B is 
2 
v 
Pa =P8| zi -žp -2 
B =P8| zı -zB 2g 
2 
-1000g x |o — 198^ 
28 


--1962 Nn? 


Finally, the pressure at C is 


"e 
Pc 7P8| zi —zc 
1 c 2g 
2 
-1000g x | 15 198" 
2g 


--16677 Nn? 


The average velocity in the Pipeline is 1.98 ms-! and the pressures at points A, 
B and C are 17.658 kNm~, -1.962 kNm= and -16.677 kNm-2, respectively. 
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2.5 Pressure drop due to enlargements 

Water flows through a pipe with an inside diameter of 5 cm at a rate of 10 m?h-! 
and expands into a pipe of inside diameter 10 cm. Determine the pressure drop 
across the pipe enlargement. 


a 
Flow p uw D Pr 
w a uiua 





Solution 
If a pipe suddenly enlarges, eddies form at the corners and there is a permanent 
and irreversible energy loss. A momentum balance across the enlargement 
gives 

Paz + PQv, = pra, + PQV2 
Rearranging, the pressure drop is therefore 
j 2201 =¥2) 


a 


p2-P 


where the average velocity in the smaller pipe is 








and in the larger pipe is 
(4 
No (3 Jn 
2[995 | i4 
01 


=0.352 ms! 
The pressure drop is therefore 
1000 x p x (141 - 0352) 
paps 
axo? 
4 
234 Nm? 
Applying the Bernoulli equation over the section, the head loss is 
2. 2 
Es ¥1 ~2 , Pi-P2 
2g ps 


| 
Y1 Y2 _ PQ) -v2) 


2g a»pg 
d aid 
LYLY? _ vaw v3) 
2g 2g 
which reduces to 
2 
H, = (vy =v) 
2g 


2 


$! 
ecl 1255. 
2g vi 


From continuity for an incompressible fluid 


aw =a 
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Then 


2 
v 

H, i2] 
8 2 


or in terms of diameter for the circular pipe 


$ 
v 

H, =|1- d 
2g dy 


Therefore 


= 0.057 m 


The pressure drop is therefore 
Apy = psi, 
=1000 x g x 0057 


=559 Nm” 


The pressure drop is 559 Nm"?. Note that for a considerable enlargement where 
dy >>a, the head loss tends to 


jus 
LU 
That is, the head loss due to an enlargement is equal to one velocity head based 
on the velocity in the smaller pipe. This is often referred to as the pipe exit head 
loss. Note that although there is a loss of energy (or head) there may not neces- 
sarily be a drop in fluid pressure because the increase in cross-section causes a 
reduction in velocity and an increase in pressure. 
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2.6 Pipe entrance head loss 
Derive an expression for the entrance loss in head form for a fluid flowing 
through a pipe abruptly entering a pipe of smaller diameter. 








The permanent and irreversible loss of head due to a sudden contraction is not 
due to the sudden contraction itself, but due to the sudden enlargement 
following the contraction. Consider, therefore, a pipe of area a, which reduces 
to areaa). The fluid flowing into the narrow pipe is further contracted forming a 
vena contracta. At this point the areaa,,. is related to the smaller pipe area by a 
coefficient of contraction as 


[Lr 


Beyond the vena contracta, the fluid expands and fills the pipe. The head loss 
due to this expansion is 


fue ve =¥2P 
2g 


and is known as the Carnot-Borda equation after the French mathematicians 
Lazare Nicolas Marguerite Carnot (1753-1823) and Jean Charles Borda 
(1733-1799). From continuity 


4592 = cV yc 
= Cea ye 
Therefore 


Ti 
Tae ES 
G 
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The constant k is found by experiment. For a sudden contraction, the head loss 
is close to 


and is usually referred to as the entrance head loss to a pipe. Experimental 
values are 


ala, 0 02. 04 06 08 10 
k 0.5 0.45 036 021 007 0 
0.5 
k 
0 
0 10 
aya, 
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2.7 Force on a pipe reducer 

Water flows through a pipe of inside diameter 200 mm ata rate of 100 m3h-!. If 
the flow abruptly enters a section reducing the pipe diameter to 150 mm, for 
which the head loss is 0.2 yelocity heads based on the smaller pipe, determine 
the force required to hold the section in position. Upstream of the reducer, the 
gauge pressure is 80 kNm-. 











aF. 
Flow 1 
= 
pep, ^ 
Solution 
The velocities in the larger and smaller pipe are 
4Q 
n=: 
2 
nd) 
4x 10 
= 3600 
nx 02? 
= 0884 ms“! 
and 
na < 
2 
ndi 
45,100. 
- 3600 
nx 015? 
=157 ms! 


The head loss at the reducer is based on the velocity in the smaller pipe as 
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2 

v 
Hy 2022- 

28 


2 
202x157. 
28 
70025 m 
The pressure in the 200 mm diameter pipe is 80 KNm~. The pressure in the 150 
mm diameter pipe is found by applying the Bernoulli equation 
2 2 
v v 
PY gi PD. 2 L 
ps 28 pg 2g 


Rearranging 
pn «Pot -v3) -pgH , 


- 80x 10? + OO v ossa? -1572) - 1000 x g x 0025 


= 78,913 Nm? 
The upstream and downstream pressure forces are therefore 
Pia, = 80x 10° x x02 
=2513 N 
and 
pray = 78,913 x 7% 015 
=1394 N 


The force in the x-direction is therefore 
Fy =pPQUV2 -v1)-p14, + p20 


= 1000x O x (1.57 -0884) -2513 + 1394 
3600 


= —1100 N 


A force of 1.1 kNm-? in the opposite direction to flow is required to hold the 
reducing section in position. 
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2.8 Vortex motion 
Derive an expression for the variation of total head across the streamlines of a 
rotating liquid. 


p+Ap 
E v+ dv 
dr Streamline 
p Í V7 streamline 


r 


Solution 

Consider an element of liquid of length L and width dr between two horizontal 
streamlines of radii r and r +dr and which have corresponding velocities v and 
v +dv. The difference in radial force is equal to the centrifugal force. That is 


2, 
pLv^dr 
yes SELS 
r 


from which the pressure head is therefore deduced to be 
Ap v?dr 


Pg er 


The radial rate of change of pressure head is therefore 


Ap 
di 
er) v? 
dr gr 


while the radial change of velocity head is 





+dv y 
2gdr 
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The rate of change of total head with radius is therefore 


dr gr gdr 


gv dp 

gir dr 
This is an important result based on a horizontal moving fluid and can be used 
to determine the variation of head (or pressure) with radius for both forced and 
free vortex motion. In a free vortex, the fluid is allowed to rotate freely such as 


inthe case of a whirlwind, flow round a sharp bend or drainage from a plughole. 
There is a constant total head across the streamline. Thus 


2 
aH vo aL uv 


dH =0 


such that vr is a constant. For free vortex flow it can be shown that the velocity 
increases and pressure decreases towards the centre. 

Ina forced vortex in which a fluid is rotated or stirred by mechanical means, 
the tangential velocity is directly proportional to the streamline radius as 


v =or 


where ( is the angular velocity. For forced vortex flow it can be shown that the 
free surface is a paraboloid. 


52 





(CONT ry, MOMENTUM AND ENERGY 


2.9 Forced and free vortices 


An impeller of diameter 50 cm rotating at 60 rpm about the vertical axis inside 
a large vessel produces a circular vortex motion in the liquid. Inside the 
impeller region the motion produces a forced vortex and a free vortex outside 
the impeller with the velocity of the forced and free vortices being assumed 
equal at the impeller edge. Determine the level of the free surface at a radius 
equal to the impeller and a considerable distance from the impeller shaft above 
the liquid surface depression. 


Free vortex Forced vortex Free vortex 


— F 


Free surface 














Datum 


Solution 


For two-dimensional flow, the rate of change of total head H with radius r is 
given by 


(see Problem 2.8, page 51) 
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For the forced vortex, the tangential velocity is related to angular velocity by 
v-2or 
and thus 
dv =a@dr 
Therefore 
an 20r 
dr g 
Integrating, the difference in total head between two streamlines of radii r and 


nis 


H. 2945 
Tuin 


n 
to give 


Hiie 
m-n -CD 


From the Bernoulli equation applied at the free surface (p; =p2), the total head 
is 
v 2 = v? 
H -H => +z -z 
2 1 2g 2-2 








oho 
B tza- 
28 
Combining the equations for total head gives 
woe - r?) 
z-z = — —L 
2g 
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Since at the centre of the paraboloid r; =0, the elevation of the free surface z at 
the radius of the impeller r; is therefore 





_ (on 
£2—z1 2g 
—QnNn Y 
2g 
2x nx 6 x025 
S 60. 
2g 
=0126 m 
For the free vortex 
dH =0 
Therefore 
Br dk g 
y r 
Integrating 


y 
log, “2 + log, 2 =0 
vi 7i 


then for the free surface 
c-vr 
where c is a constant. At the edge of the impeller (7; ) this is 
c -v3n 
- ar? 
= 2nNrp 


=2x rx É 025? 
60 


= 0.393 
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The tangential velocity at the edge of the impeller is 


c 
v => 
n 


and at a distant point from the impeller is 


d 
v == 
25 
Applying the Bernoulli equation at the free surface (p3 =p3) then 


2 
v3 - 


28 


2 2 





237723 = 





Since the radius rs is large then 
c 





23 7227 


2 
2gr; 























0.3937 
“Ox g x025? 
=0126 m 
The total depression in the rotating liquid is therefore 
z3 -Z1 = (z3 — 22) + (22 — 21) 
= 0.126 m+ 0126 m 
-0252m 


The total depression is found to be 25.2 cm. 
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Further problems 

(1) Water flows upwards through a pipe which tapers from a diameter of 200 
mm to 150 mm overa distance of 1 m. Neglecting friction, determine the rate of 
flow if the gauge pressure at the 200 mm section is 200 kNm^? and at the 150 
mm section is 150 kNm~2, Water has a density of 1000 kgm^?. 


Answer: 0.192 m3s-! 


(2) Show that for a liquid freely discharging vertically downwards from the 
end of a pipe, the cross-sectional area of the vena contracta, az, to pipe area, aj, 
separated by a distance z is 


(3) A pipe of inside diameter 100 mm is suddenly enlarged to a diameter of 
200 mm. Determine the loss of head due to this enlargement for a rate of flow of 
0,05 m3s~!. 


Answer: 1.16 m 


(4) Water is discharged from a tank through an external cylindrical mouth- 
piece with an area of 100 cm? under a pressure of 30 kNm-2. Determine the rate 
of discharge if the coefficient of contraction is 0.64. 


Answer: 0.0645 m?s-! 


(5)  Wateris added to a process vessel in the form of a jet and directed perpen- 
dicularly against a flat plate. If the diameter of the jet is 25 mm and the jet 
velocity is 10 ms-!, determine the power of the jet and the magnitude of the 
force acting on the plate. 


Answer: 491 W, 49,1 N 
(6) A process liquid of density 1039 kgm~ is fed continuously into a vessel 


as a jetat a rate of 12 m?h-. If the jet, which has a diameter of 25 mm, impinges 
on a flat surface at an angle of 60° to the jet, determine the force on the plate. 


Answer: 20.4 N 





(7)  Ajetof water 50 mm in diameter with a velocity of 10 ms! strikes a series 
of flat plates normally. If the plates are moving in the same direction as the jet with 
a velocity of 7 ms-!, determine the pressure on the plates and the work done. 


Answer: 30 kNm2, 412 W 
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(8) Show that the efficiency n of a simple water wheel consisting of flat 
plates attached radially around the circumference in which a fluid impinges 
tangentially is 
21 ~va)v2 
i 


where v, is the velocity of the jet and v; is the velocity of the plates. 


(9) Determine the efficiency of the plates in Further Problem (7). 


Answer: 42% 


(10) Show that the maximum efficiency of the water wheel described in 
Further Problem (8) in which a jet impinges normally on its flat vanes is 50%. 


(11) Water is discharged through a horizontal nozzle at a rate of 25 litres per 
second, If the nozzle converges from a diameter of 50 mm to 25 mm and the 
water is discharged to atmosphere, determine the pressure at the inlet to the 
nozzle and the force required to hold the nozzle in position. 


Answer: 150 kNm-?, 133 N 


(12) Two horizontal pipes with inside diameters of 4 cm are connected by a 
smooth horizontal 90° elbow. Determine the magnitude and direction of the 
horizontal component of the force which is required to hold the elbow at a 
flowrate of 25 litres per second and at a gauge pressure of 3 atmospheres. Air at 
atmospheric pressure surrounds the pipes. 


Answer: 1063 N, 45° 
(13) A pipeline of inside diameter 30 cm carries crude oil of density 920 


kgm- at a rate of 500 m?h-!. Determine the force on a 45° horizontal elbow if 
the pressure in the elbow is constant at 80 KNm-?. 


Answer: 4122 N, 22.5° 
(14) A 200 mm inside diameter pipe carries a process liquid of density 1017 


kgm^? at a rate of 200 m?h7!. Determine the magnitude and direction of the 
force acting on a 90? elbow due to momentum change only. 


Answer: 141 N, 45° 
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(15) A pipeline with a diameter of 90 cm carries water with an average 
velocity of 3 ms~!. Determine the magnitude and direction of the force acting 
on a 90* bend due to momentum change. 


Answer: 8097 N, 45° 
(16) Acetic acid with a density of 1070 kgm™ flows along a pipeline at a rate 
of 54 m?h-!, The pipeline has an inside diameter of 100 mm and rises to an 
elevation of 5 m. Determine the kinetic energy per unit volume of the acid and 


the pressure head at the elevated point if the gauge pressure at the lower eleva- 
tion is 125 kNm?. 


Answer: 1951 Jm, 16.9 m 
(17) Distinguish between a free and forced vortex. 


(18) Show that the surface of a liquid stirred within a cylindrical vessel 
forming a forced vortex is a paraboloid. 


(19) Determine the difference in pressure between radii of 12 em and 6 cm of a 
forced vortex rotated at 1000 rpm. 


Answer: 59.2 kNm? 


(20) In a free cylindrical vortex of water, the pressure is found to be 200 
kNm~ at a radius of 6 cm and tangential velocity of 6 ms~!. Determine the 
pressure at a radius of 12 cm. 


Answer: 213.5 kNm^? 


(21). A cylindrical tank of radius 1 m contains a liquid of density 1100 kgm? 
to a depth of 1 m. The liquid is stirred by a long paddle of diameter 60 cm, the 
axis of which lies along the axis of the tank. Determine the kinetic energy of the 
liquid per unit depth when the speed of rotation of the paddle is 45 rpm. 


Answer: 904 Jm7! 


(22) Comment on the consequence of assuming inviscid fluid flow in terms of 
flow through pipes and vortex motion in tanks. 
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Laminar flow 
and lubrication 


Introduction 

Flowing fluids may exhibit one of two types of flow behaviour that can be 
readily distinguished. In streamline or laminar flow, fluid particles move along 
smooth parallel paths or layers (laminae) in the direction of flow with only 
minor movement across the streamlines caused by diffusion. 

Over the centuries, the existence of laminar and turbulent flow has been 
studied extensively by many prominent scientists. In 1839, it was first noted by 
the German hydraulics engineer Gotthilf Heinrich Ludwig Hagen that laminar 
flow ceased when the velocity of a flowing fluid increased beyond a certain 
limit. With much work on the subject over the following three decades, Hagen 
finally concluded in 1869 that the transition from laminar to turbulent flow was 
dependent on velocity, viscosity and pipe diameter. Around the same time, 
French physician and physicist, Jean Louis Marie Poiseuille, whilst researching 
the effects of blood flow in veins, reported similar independent work on the 
viscosity and pressure drop of water in capillaries and reached the same mathe- 
matical conclusions as Hagen. 

It was not until 1883 that British scientist Osborne Reynolds showed by both 
dimensional analysis and experiment that the transition depends on a single 
dimensionless parameter which bears his name. This is given by 





where Reis the Reynolds number, p is the density, v the average velocity and 
u the viscosity of the fluid, and d is the diameter of pipe. 

Reynolds” experiments involved injecting a trace of coloured liquid into the 
flow of water in a horizontal glass tube. At low flowrates the coloured liquid 
was observed to remain as discrete filaments along the tube axis, indicating 
flow in parallel streams. At increased flow, oscillations were observed in the 
filaments which eventually broke up and dispersed across the tube. The critical 
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value of Reynolds number for the break-up of laminar flow was found to be 
about 2000 while turbulent flow was not found to occur until above 4000. The 
so-called transition zone, which lies between Reynolds numbers of 2000 and 
4000, was found to be a region of fluid streamline instability and 
unpredictability. 

The majority of mathematical problems involving laminar flow tend to be 
straightforward. The basis of most calculations involves applying a simple 
equilibrium force balance to the flowing fluid with fully developed laminar 
flow. Without considering inertial forces, the acceleration of particles in the 
fluid caused by the pressure gradient is retarded by the viscous shear stresses set 
up by the velocity gradient perpendicular to the direction of flow. This 
approach forms the basis of the independent work of both Hagen and Poiseuille, 
from whence it is possible to determine important flow details such as velocity 
profiles and rates of flow through pipes, gaps and channels as well as the 
viscous and lubricating effects of fluids in bearings. 
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3.1 Reynolds number equations 

Establish equations for Reynolds number for the flow of a fluid with average 
velocity v in a pipe of inside diameter d, in terms of kinematic viscosity, volu- 
metric flowrate, mass flowrate and mass loading. 


Solution 
The Reynolds number is an important dimensionless parameter and is valuable 
for identifying whether the flow of a fluid is either laminar or turbulent. The 
Reynolds number is traditionally given as 
"d 
p P 


u 
It can be alternatively expressed in terms of kinematic viscosity which is 


T 


P 
to give a Reynolds number expression of 


fan s 
v 


Since the average velocity v is related to volumetric flowrate by 





Likewise, since volumetric flowrate is related to mass flowrate by density 
m 
Q=" 

p 


then Reynolds number is 


oa. 
und 
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Finally, since mass loading is related to mass flowrate by flow area 
4m 
nd? 

The Reynolds number is therefore 


Ree 


u 


It should be noted that the Reynolds number is of primary importance in consid- 
ering the flow of fluids through pipes and open channels, and around objects. 
As a dimensionless parameter the Reynolds number (Re = pvd /u ) represents 
the ratio of inertial to viscous forces in a flowing fluid and is directly ass ed 
with the boundary surface over which the fluid passes. For laminar flow in 
circular pipes (Re« 2000), the viscous forces dominate while the inertial forces 
are of little significance. In turbulent flow the reverse is true with the laminar 
sublayer being destroyed and the viscous forces being of little significance. 

The Reynolds number is also of importance in considering the viscous drag 
on submerged objects. For a Reynolds number below a value of 0.5 the precise 
shape of the object is of less importance than the viscosity of the fluid and the 
velocity of flow. For such conditions it is possible to derive a relationship for 
the drag on objects in which the inertial effects are not included. For Reynolds 
numbers above 0.5, however, the reverse is again true — the influence of 
viscosity diminishes while the influence of inertia progressively dominates 
with increasing values of Re. 
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3.2 Laminar boundary layer 
Determine the furthest distance into a pipe with an inside diameter of 8 mm and 
a well-rounded entrance at which fully developed laminar flow begins. 


Developing boundary layer 


Fully developed laminar flow 






_ 


Solution 

When a fluid with uniform velocity flow enters a well-rounded pipe, entrance 
particles close to the pipe wall are brought to rest due to the viscous properties 
of the fluid. This is known as the no-slip condition. A steep velocity gradient 
then exists in the fluid varying from zero at the pipe wall to the uniform velocity 
in the bulk of the fluid. As the fluid moves further into the pipe, the viscous 
retardation of the particles in adjacent layers gradually increases in thickness, 
For the steady flow through the pipe, the fluid near the centre of the pipe accel- 
erates until an equilibrium condition is reached for fully developed laminar 
flow with a parabolic variation of velocity with pipe radius. The region where 
the velocity is changing and one where it is uniform is known as the laminar 
boundary layer. Since the uniform velocity is approached asymptotically, the 
edge of the boundary layer is defined as a point where the fluid velocity reaches 
995; of its theoretical maximum value, The region prior to the section at which 
laminar flow is fully developed is known as the entrance transition length L’, for 
which a theoretical formula proposed by Henry L. Langhaar in 1942 and 
supported by good experimental agreement is given by 


L’ =0058Red 


Since the maximum entrance transition length occurs for a Reynolds number of 
2000 then 


L’ = 0058 x 2000 x 0.008 
=0928 m 


The furthest distance the fluid can flow into the 8 mm inside diameter pipe 
before fully developed laminar flow can exist is 0.928 m. 


65 









































— —— M — 
FLUID MECHANICS LAMINAR FLOW AND LUBRICATION 





3.3 Velocity profile in a pipe and applying the no-slip condition at the wall (v , =O at r=R) then 
Derive an expression for the local velocity of a fluid flowing with fully devel- v. -1 Ap & 
oped laminar flow through a horizontal pipe of radius R and sketch the velocity favs "ODE J rdr 
profile. 0 R 
where Ap is the pressure drop over the length of element (p; —p2). 
Integration gives 
B 
i ET [=] 
ae 
mL 2 Ip 
= 1 APR? - 2) 
4 L 











The velocity of the fluid therefore has a parabolic variation with pipe radius. 





Solution 


Consider a cylindrical element of the fluid to be well away from the wall on the 
centre line of the pipe where the difference of pressure force across the ends 
provides the driving force to overcome the friction on its outer surface. For 
uniform flow, an equilibrium force balance on the element is therefore 





(pı -pa)nr? =QrrdL 


That is Note that the maximum velocity of the fluid, Vmax» Occurs at the furthest 
Š distance from the wall (r =0). That is 
pnr? (n rat yn -qQndL d Ap p2 
tee = 
ALL 
The viscous shear stress is given by 
dv, 
t=p— 
E dr 
Let 
PAP 
6L L 
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3.4 Hagen-Poiseuille equation for laminar flow 

in a pipe 

A process vessel is to be supplied with glycerol of SG 1.26 and viscosity 1.2 
Nsm~ from a storage tank 60 m away. A pipe with an inside diameter of 12.6 
cm is available together with a pump capable of developing a delivery pressure 
of 90 kNm- over a wide range of flows. Determine the glycerol delivery rate 
using this equipment. Confirm that the flow is laminar. 








Solution 
To determine the rate of flow of the fluid through the pipe with laminar flow 
consider a thin annular element between radius r and r+dr for which the 
velocity can be taken as constant. The elemental volumetric flowrate is 
therefore 


dQ =v, 2mrdr 
The total rate of flow can therefore be found by integrating across the entire 
pipe radius 
R 
feo = fv, 2nrdr 
0 0 


Since the local velocity is related to radius by 


14,2. 2 
z——É (Rr see Problem 3.3, page 66) 
iar T ( pagi 


the integration is therefore 


v 


N R 
fao =Z fer? -r?)rar 
0 2n L 0 
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which gives 
o-ze n] 
Q L| 2 Acl 
= TAP pA 
8u L 


This equation was derived independently by both Hagen in 1839 and Poiseuille 
in 1840 and is known as the Hagen-Poiseuille equation. The rate of flow is 
therefore 





x 90x10). (0126 Y 
LR 90x10 [0126 
8x12 ^ 60 2 
= 0.00773 ms! 


As a check for laminar flow, the Reynolds number is 


Re= #2 
und 


4 x 1260 x 000773 
- ""IAGERGOIDG. (see Problem 3.1, page 63) 


=82 


and is below 2000, therefore confirming laminar flow. 
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3.5 Pipe diameter for laminar flow 

In the design of a small-scale bioprocess plant, a shear sensitive Newtonian 
liquid of density 1100 kgm-? and viscosity 0.015 Nsm~ is to be transported 
along a length of pipe at a rate of 4 litres per minute. Determine the diameter of 
the pipe required if the pressure drop along the pipe is not to exceed 100 Nm-2 
per metre length. Comment on the value of the Reynolds number. 


Solution 
Assuming laminar flow, the Hagen-Poiseuille equation is 
Q- e * Rt (see Problem 3.4, page 68) 


Rearranging, the pipe radius is 
Bu 
R=| 22 
S 
L 


8x 0015 x 2004 V 
—— 501. 
zx100 


—00126 m 
The diameter of the pipe is therefore 25.2 mm. To confirm laminar flow, the 
Reynolds number is 
Re = apo 
nud 
0004 


4x1100x —— 
60 


-A————ÉM— Prol -1, 
xx 0015 x 00252 (see Problem 3.1, page 63) 


=247 


The flow is therefore laminar with a Reynolds number of 247. 
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3.6 Laminar flow through a tapered tube 

A lubricating oil of viscosity 0.03 Nsm~ is delivered to a machine at a rate of 
1077 mis-! through a convergent, tapered tube of length 50 cm with an 
upstream diameter of 10 mm and downstream diameter of 5 mm. Determine the 
pressure differential which will maintain the flow. Entrance and exit losses, 
and inertia effects due to the change in velocity in the tube, may be neglected. 


L1 Not to scale 


5mm 








50cm. 


Solution 
For laminar flow, the rate of flow is given by the Hagen-Poiseuille equation 
Q-— Ap gà 


(see Problem 3.4, page 68) 
UL pag 


Over the length of the tapered tube, the differential pressure Ap increases with 
decreasing tube radius. From geometry 

R - 0005 (1 — L) 
where R and Lare both measured in metres. Substituting for R, the total pressure 
differential is therefore obtained by integrating over the length of tube. That is 

05 

= MQ J dL 

20005* 5 (1— L)* 








This is made easier using the substitution 





zB . 
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Thus 
du = -dL 
The integration then becomes 
-8uQ [du 
^ n005* 5 f 


Integrating with respect to u gives 


Q u> 
x0005^ 3 





That is 


0.5 
EESO. 1 
3n0005* | (1 —1)° 


8 «003 x 1077 ( 1 
= -1 





3xnx000^ (0-05? 


2285 Nm? 


The pressure differential is found to be 28.5 Nm-?. 
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3.7 Relationship between average and maximum 
velocity in a pipe 


Deduce the relationship between the average and maximum velocity for a fluid 
flowing with fully developed laminar flow through a horizontal pipe. 


Solution 
The average velocity of the fluid Q can be determined from the total rate of 
flow through the total cross-sectional area a available for flow 


Q 


v=% 


a 
Since the flow Q is given by the Hagen-Poiseuille equation (see Problem 3.4, 
page 68), then 
n Ap rz 
n 8u L 
=LA p 
8u L 
Also, since the maximum velocity occurs at the furthest point from the wall 
(r =0) (see Problem 3.3, page 66) 





eee Se 
mæ =a L 
then the ratio of average velocity to maximum velocity is 
1 Ap R 
v uL 
Vma 1AP p2 
ALL 
=t 
2 


That is, the average velocity of a flowing fluid with laminar flow is half the 
maximum velocity. 
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RADIUS 








0 VELOCITY 


Note that the average velocity can be alternatively expressed in terms of pipe 
diameter as 

i ms Ap d? 

324 L 
and maximum velocity as 
1 Ap 2 
v z——— d 
"* 60 L 
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3.8 Relationship between local and maximum 

velocity in a pipe 

Fuel oil of density 900 kgm? and viscosity 0.1 Nsm? flows through a pipe of 
inside diameter 100 mm. The velocity at the centre line is 2 ms-!. Derive an 
expression relating the local to maximum velocity of a fluid flowing with fully 
developed laminar flow through a horizontal pipe, and determine the shear 
stress in the fuel oil at the pipe wall and the pressure drop per unit length of 


pipe. 


Solution 
The local velocity is given by 


Vy re al =r?) (see Problem 3.3, page 66) 


and the maximum velocity occurs at the furthest point from the pipe wall (r =0). 


1 
ae Eid 
Therefore 
1 AP (pe spa 
Di Ml 0o. 
V max 1 App 
4u L 
2 


pes 
R 
Since the viscous shear stress is 
dv, 
t=p 
E dr 


the velocity gradient is therefore obtained by differentiation to give 





dv, r 
ar T mur 
Thus 


r 
T= -WV max Eus 
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At the pipe wall (r =R), the shear stress is therefore 


ty = me 


-2x01x2 
0.05 


(98 Nn? 
The negative sign indicates that the shear stress is in the opposite direction to 
flow. To verify laminar flow, the Reynolds number is. 

Re- pud 


* 
Mo 
u 
900x 2 x 01 
01 
- 900 


which is less than 2000. Finally, the pressure drop along the pipe is obtained 
from 





Ap Ava 
L R? 
_4x01x2 
005? 
-320 Nm? n^! 


The pipe wall shear stress is found to be 8 Nm? and pressure drop to be 320 
Nm per metre length of pipe. 
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3.9 Maximum pipe diameter for laminar flow 

A light oil of viscosity 0.032 Nsm~ and SG 0.854 is to be transferred from one 
vessel to another through a horizontal pipe. If the pressure drop along the pipe 
is not to exceed 150 Nm? per metre length to ensure a satisfactory flow and 
flow is always to be laminar, determine the maximum possible internal diam- 
eter of the pipe. 


Solution 
The maximum value of Reynolds number at which laminar flow can exist is 
2000. That is 


where the average velocity is given by 


Ap gà 
Su L 


Substituting, the Reynolds number is therefore 
1 Ap(dy 
su L2 


cm 
324 


or rearranging in terms of pipe diameter 


v= (see Problem 3.7, page 73) 


Re= 


Eo 


d? 


t! 


L 


1 
3 


| 328? fe 


Ap 
Pp 


1 
_ { 32x 0.0327 x 2000 )* 
854 x 150 


=008 m 


The maximum inside diameter is found to be 8 cm. 
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3.10 Vertical pipe flow 

A viscous liquid, of density 1000 kgm and viscosity 0.1 Nsm 2, flows between 
two storage tanks through 10 m of vertical pipe which has an inside diameter of 
5 em. The rate of flow is controlled by a valve at a rate of 14 m?hr!. Determine 
the pressure drop across the valve when the difference in levels between the two 
storage tanks is 12 m. Both ranks are open to atmosphere. Neglect losses in the 


pipe fittings. 











Solution 

Gravity has no effect for laminar flow in a horizontal pipe. For vertical flow, 
however, gravitational effects are in the form of hydrostatic pressure where, in 
this case, the pressure drop across the valve is 

128QuL 


APyaive 7 pg(L t H)— zi 





gx 14 x 01x10 


= 1000 g x (10 + 2)— 3 
xx005 
— 92,368 Nm? 


The pressure drop across the valve is found to be 92.4 kNm-?. 
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3.11 Film thickness in a channel 

Derive expressions for the velocity distribution, total flow and film thickness of 
a Newtonian liquid flowing under gravity down a surface inclined at an angle ® 
to the horizontal. Assume that the flow is laminar and fully developed. 





Solution 


Assuming steady flow, a negligible shear stress at the gas-liquid interface and a 
uniform pressure in the direction of flow (x-direction), an equilibrium force 
balance on the element gives 


WL(8 — y)pg sin 9 - TWL 
where the shear stress is given by 


tz da 
iS 


The velocity profile through the film can be found using the no-slip boundary 
condition (y =0 at v , =0). Therefore 


vy na 
jov. pg sing 
0 


j (6-06 
0 





Integrating 


y, zBESHO S. 
u 
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The volume flowrate of the liquid is given by 


è 
Jao [way 
0 0 


Therefore substituting for v y 


" ô 
à enel - 





E b 
Integrating 
3 2 ô 
Wpg sin@) dy" y 
277 [2 7€ 
B 0 
.. Wps sin 8 è 
3y 


Rearranging, the film thickness is therefore related to flowrate by 


1 
3 
5=|— #0 _ 
Wpg sin 8 
The film thickness can also be related to the average velocity by 


= 2 
"ow 


_ Wpg sin 0 ri 
3uóW 


= P8 sin 8 8 
3u 


Rearranging, the thickness of the film of liquid is therefore 


1 
"ere 
pg sin 8 


This equation has been found by experiment to hold for a Reynolds number 
(Re = pv6/1) less than 500. 
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3.12 Flow down an inclined plate 

In the first of two experiments water, of density 1000 kgmr? and viscosity 0.001 
Nsm, was pumped to the top of a flat plate and evenly distributed along one 
side of it. The plate was angled at approximately 60? to the horizontal and the 
flowrate was 3x1076 m?s-! per metre length of plate. It was found that the film 
thickness was 1.016x107* m. In the second experiment, a different liquid of 
density 1500 kgnr? was pumped to the top of the plate at 31079 m?s-! per 
metre length of plate and the film thickness was found to be 1.25x10- m. 
Determine the viscosity of the second liquid and the exact angle of the plate. 


Solution 
The rate of flow of liquid down the plate is given by 
Woe si 
gpa a 8 (see Problem 3.11, page 79) 


For the first experiment involving water, the exact angle of inclination is found 
by rearranging the equation to 


sing- —u. 
Wped? 
Lo 3x3x10 x1x 10% 
1x 1000 x g x (1016 x 1075) 
-0875 


The exact angle of inclination was therefore 
6-61* 


For the second experiment, the viscosity is found by rearranging the equation 
for flow to 


p = Weg sin® 5 
3 
AOO ADES aae a 
3x3x10% 


= 2.795 x 10-3 Nem 


The exact angle was 61° and the viscosity of the second liquid was found to be 
2.795 x 10-3 Nsm~?. 
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3.13 Flow down a vertical wire 

A film of a Newtonian fluid flows down the outside of a stationary vertical wire 
which has a radius R. Deduce an expression for the local velocity of the fluid, 
Vz, down the wire. 


= 








Film Wire 


Solution 
Consider a force balance on an element of the fluid between radii r and R +8 
assuming no shear effects at the free surface. Then 


PgLT((R + 3° =r?) = mL 
where the viscous shear stress is given by 


t= dv. 
dr 
Substituting and noting that the fluid is stationary (no-slip condition) at the wire 
surface (v , =O at r =R), the local velocity at radius r measured from the axis of 
the wire can be found from 
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“a Pei R+ 
je Pd 2 


to th 
-PE (gg? [9^ — prar 
9i R'OR 


Integrating with respect to r gives 


a7 
Pg 2 r ir 

z =| (R+ lo; =|= 

vz a 8) log, [r], FI, 


= far +8? io, (5 a = mj 


While the velocity of the fluid is zero at the wire surface, itis at its greatest at the 
free surface. The average velocity, rate of flow and shear stress within the fluid 
can similarly be determined and are useful results when considering the process 
of wire coating. This involves the drawing of a wire vertically up through a bath 
of viscous liquid taking with it a film of the liquid. In practice, however, the 
coating liquid — such as a polymer melt — is not usually likely to exhibit 
Newtonian behaviour and, on drying, density change and shrinkage must also 
be taken into consideration. 


83 


















































FLUID MECHANICS 


3.14 Flow and local velocity through a gap 

Derive an expression for the local velocity and total flow of a Newtonian fluid 
with fully developed laminar flow between two flat horizontal plates separated 
by a distance 2H. 

















Solution 
Assuming that the plates are sufficiently wide that edge effects can be neglected 
(that is, W>>2H ), an equilibrium force balance on the element is 


2yW(p, — p3) 21WdL 
That is 


öp 
2yW = <—dL ||=2tWdL 
y Z (n + SL J T 


where the shear stress is given by 
dv 


tape 
u dy 

Let 
3 oM 
6L L 


The local velocity v, at a distance y can therefore be found applying the 
no-slip boundary condition (v , =O at y =H). That is 


Integration with respect to y gives 
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The velocity of the fluid therefore has a parabolic variation with distance 
between the plates, with the maximum velocity at y =0. 


DISTANCE 


i 
VELOCITY = 





The total flow of fluid can be determined by considering the sum of the 
elemental flows through a narrow slit of thickness dy. From symmetry about 
the mid-point between the plates there are two such narrow slits. Thus 


H 
feo -2w [v.dy 
0 0 


Substituting for v y 





In addition to assuming that edge effects can be neglected, for incompressible 
viscous flow of Newtonian fluids through gaps or between plates it is also 
assumed that there is no pressure gradient across the fluid layer and that gravita- 


tional forces are negligible. 
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3.15 Relationship between local and average velocity 
through a gap 

Derive an expression relating the local velocity of a Newtonian fluid flowing 
with fully developed laminar flow between two horizontal parallel plates sepa- 
rated by a distance 2H in terms of the average velocity. 


Solution 
The total flow is related to average velocity by 
Q -2WHv 


The average velocity is therefore obtained by substituting an expression for 
laminar flow between the plates (see Problem 3.14, page 84) 


Q -2W Ap 5 
3u L 
to give 
aw AP 
pre a 
2WH 


s18 H? 
3p L 
The local velocity v , between the plates is given by 
1 2 
Vx =i 2 -y?) 


(see Problem 3.14, page 84) 
24 


Therefore 


1 Ap 

vy UL 
v 1 
3 


The local velocity is therefore related to average velocity by 


7) 


where y is the distance measured away from the middle of the gap. 
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3.16 Relationship between average and maximum 
velocity through a gap 

Deduce the ratio of maximum velocity to average velocity for a Newtonian fluid 
flowing with laminar flow through a gap between two flat horizontal plates. 


Solution 
The maximum velocity occurs at the furthest point from either surface (y —0) 
(see Problem 3.14, page 84). That is 


1 MP 2 
mL 


Since the average velocity is given by 


V max 


(see Problem 3.15, page 86) 


The maximum velocity is one and a half times the average velocity through the 
gap. 


DISTANCE 
H 


H 


VELOCITY 
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3.17 Shear stress for flow through a gap 

Oil, with a viscosity of 1.5 Nsm-?, flows with fully developed laminar flow 
through a gap formed by two horizontal parallel plates set 20 mm apart. Deter- 
mine the magnitude and direction of the shear stresses that act on the plates 
when the average velocity is 0.5 ms-!. 


Solution 
The velocity of the fluid is given by 
sies ni | (see Problem 3.15, page 86) 
E^ H 


The velocity gradient is therefore 


st) 











The shear stress at the wall is 225 Nm? and acts in the opposite direction to 
flow. 
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3.18 Flat disc viscometer 

A device used to measure the viscosity of a viscous liquid consists of a flat disc 
which rotates on a flat surface between which is sandwiched the liquid under 
investigation. The disc has a diameter of 5 cm and produces a shear stress of 
400 Nm^ for a rotational speed of 600 rpm. Determine the viscosity of the 
liquid and the torque on the rotating disc if the clearance between the disc and 
surface is 2 mm. 





Solution 
The viscosity of liquids can be determined experimentally using a device 
known as a viscometer. There are various types of viscometer available and all 
are based on the principles of laminar flow. The flat disc viscometer is perhaps 
one of the simplest forms, and can be used for quick estimates of viscosity. For 
the liquid sandwiched between the surface and the rotating disc, the shear stress 
is given by 

dv, 

dy 





T= 


or in terms of angular velocity and clearance between disc and plate 


or 
t=p— 
I 


where 
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0 -2nN 


-2x nx 800 
60 


= 628 rads! 


Rearranging in terms of viscosity 


_ tt 
y or 


_ 400 x 2x 109 
62.8 x 0025 


— 0509 Nsm? 


The torque is given by 
R 
T= nj wrdr 
0 
Substituting for shear stress 


R 
7 = HOT Bap 
o 


Integrating 


_ Tx 0509 x 628 


4 
2x2x10> ^ DS 


=98 x 1079 Nm 


The viscosity and torque are found to be 0.509 Nsm~ and 9.8x10? Nm. 
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3.19 Torque on a lubricated shaft 

A shaft 100 mm in diameter rotates at 30 rps in a bearing of length 200 mm, the 
surfaces being separated by a film of oil 0.02 mm thick. Determine the viscous 
torque on the bearing and power required to overcome the frictional resistance 
if the viscosity of the oil is 0.153 Nsm-?. 








Solution 
It is assumed that there are negligible end effects and that the shaft is lightly 
loaded so that the shaft runs concentrically. The resistance per unit area or shear 
stress is 


The oil film is very thin such that the viscous shear stress may be given in terms 
of angular velocity at radius r 





The tangential resistance (viscous drag) on the bearing is 


F=QnrL 
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The resisting torque is therefore 
T=Fr 
- Qnrir 


gun 2nrLr 





=u 
Ax r^ NL 
=e 


4x r? x 005? x 30x 02 


-0153x = 
2x10 


-226 Nm 
The power of the shaft is 
P=To 
=T2nN 
=226x2x mx 30 
= 42,600 W 


The torque on the bearing is found to be 226 Nm and the power required to 
overcome the frictional resistance is 42.6 kW. 

Note that if the radial clearance between the shaft and bearing is not small, 
the adjustment to torque is made 


Ax rer NL 


no-n 


where r, and rs are the radii of the shaft and bearing, respectively. 
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3.20 Lubricated collar bearing 

The axial thrust on a shaft is taken by a collar bearing, the outside and inside 
diameters of which are 200 mm and 125 mm, respectively. The bearing is 
provided with a forced lubrication system which maintains a film of oil of thick- 
ness 0.5 mm between the bearing surfaces. At a shaft speed of 300 rpm the. fric- 
tion loss of power amounts to 22 W. Determine the viscosity of the oil. 






Not to scale 





Bearings are mechanical devices used to reduce the friction between two parts 
in contact with one another and which move with respect to one another. Collar 
bearings are to be found on vertical rotating shafts and take a load on a 
stationary support. To reduce friction, a film of lubricating fluid separates the 
rotating collar and stationary support. The viscous resistance of the collar 
bearing can be obtained by assuming the face of the collar to be separated from 
the bearing surface by a thin film of oil of uniform thickness where the viscous 
shear stress is given by 


dv 
tap 
Bay 


For the thin film of oil, the shear stress may therefore be expressed in terms of 
angular velocity as 


rap" 
t 
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The tangential viscous force changes with radius. For a thin ring at a radius r 
- Qnrdr 
=u onrdr 
t 
The moment of tangential force on the ring is therefore 
=u OY onrdr r 
t 
For the whole bearing to consist of similar concentric rings, the total torque 
required to overcome the viscous resistance of the bearing is found by 
integrating 


= 2 ag 
Ry 





= - R$) 


The power loss is related to torque by 
P=To 
Therefore rearranging in terms of viscosity 


—— 
nor (R} - RD) 


2 x 0.0005 x 22 


wx(2x ax aj x (014 —00625*) 


= 00837 Nsm? 


The viscosity of the oil is found to be 0.0837 Nsm~?. 
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Further problems 


(1) Explain what is meant by a Newtonian fluid. 


(2) A cubic block of volume 1000 cm? is supported by a film of oil 0.1 mm 
thick on a horizontal surface. Determine the viscous drag if the block is moved 
across the surface with a velocity of 1 ms-!. The viscosity of the oil is 0.05 
Nsm^. 


Answer: 5N 


(3) A block is supported by a film of oil 0.2 mm thick on a horizontal surface. 
The area of the oil film is 100 cm2, the viscosity of the oil is 0.05 Nsm"?. Deter- 
mine the viscous drag at a velocity of 1 ms-!. 


Answer: 2.5 N 


(4) A fluid flows at a rate Q in a pipe of diameter d, and has a Reynolds 
number Re. If the fluid passes to a pipe of diameter d» , deduce the relationship 
between Re; and Rez. 


(5) Show that the local velocity of a liquid with laminar flow along a hori- 
zontal channel of width W in terms of flowrate Q and depth Scan be given by 


ai 


2 
where y is the depth measured from the base of the channel. 


(6) Show that the maximum velocity of a liquid with laminar flow along a 
horizontal channel of width W and depth can be given by 


V max 


2w 


(7) Sketch the velocity profile and write down the boundary conditions for 
the laminar movement of a Newtonian liquid between two large flat plates 
which are horizontal, parallel and 1 mm apart, if the top plate moves at +0.01 
ms-! and the lower plate moves at —0.01 ms-! and there is no net flow of fluid 
between the plates. 


(8) Outline the assumptions on which the flow of fluid with fully developed 
laminar flow through a pipe is based. 
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(9) Glycerol, of SG 1.26 and viscosity 1.2 Nsm~2, flows between two flat 
plates 2 cm apart. If the pressure drop along the plates is 500 Nm per metre 
length, determine the maximum velocity of the glycerol. 


Answer: 0.021 ms! 


(10) Show that the local velocity of a Newtonian fluid with laminar flow 
flowing through an annular space between two tubes can be given by 


we [| 
y= (Rr? = 2) LBL 


ve | 


where R is the inner radius of the outer tube, r is the outer radius of the inner 
tube, y is a radius between R and r, 4 is the dynamic viscosity of the fluid and 
Ap/L is the pressure drop per unit length through the annulus. 


1 Ap 
;ı =——]| (R 
Vy am 








(11) Sketch the velocity profile of a Newtonian fluid through an annular gap 
with laminar flow. 


(12) Sketch the velocity profiles for the flow of two immiscible liquids A and 
B along a horizontal channel, The viscosities and thickness of layers may be 
assumed to be identical and where the density of liquid A is less than the density 
of liquid B. 


(13) Show that the flowrate of the Newtonian liquid flowing down the outside 
of a vertical stationary wire of radius R may be given by 


5 4 
eG ae Atte) 








4 


where dis the thickness of the liquid film. 
(14) Show that the local velocity v, of Newtonian liquid flowing with fully 


developed laminar flow down a plate inclined at an angle @ to the horizontal 
may be given by the expression 


pg sin®}(  3Qu 
y, PETS OMBRE [VL 
H Pg sin 8 2 


where y is the depth measured from the surface of the plate. 
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(15) Derive an expression for the velocity distribution for a Newtonian fluid 
flowing with fully developed laminar flow through a narrow gap formed by two 
vertical parallel plates and sketch qualitatively the velocity and shear stress 
profiles. 


(16) A liquid flows in laminar flow along a horizontal rectangular channel 
which has a width W and a depth à, where W 26. Show that well away from the 
entrance, exit and side walls of the channel, the liquid velocity along the 
channel is given by 


60 


vy = y =y?) 


where Q is the volumetric rate of flow and y is the distance from the bottom of 
the channel. 


(17). A liquid flows in laminar flow at a rate Q along a horizontal tube which 
has an inside diameter d. Show that under steady flow conditions and far away 
from the entrance to or exit from the tube, the liquid velocity v , at a radius r is 
given by the expression 


where Q is the volumetric rate of flow and d is the diameter of the tube. 


(18) A cylindrical bob with a diameter of 100 mm and length 150 mm revolves 
within a cup with a clearance of 1 mm with the annular space filled with glyc- 
erine. If the bob rotates at 120 rpm for which a torque of 2.3 Nm is recorded, 
determine the viscosity of the glycerine. 


Answer: 0.079 Nsm? 


(19) Water flows in laminar flow along a horizontal channel which has a width 
of 4 m. It is observed that the depth of water in the channel is 3 cm and that 
small gas bubbles, which form at the floor of the channel, take 220 seconds to 
reach the surface after detachment. It is also observed that the bubbles reach the 
water surface 7.32 m, measured in the direction of flow, from the point of 
formation. Determine the rate of flow of water along the channel. 


Answer: 0.004 m3s-! 
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Dimensional 
analysis 


Introduction 

Dimensional analysis is a particularly useful way of obtaining relationships 
between variables that predict the occurrence of natural phenomena useful for 
scaling up models. The techniques involved, also to be found under the titles of 
similitude, theory of dimensions, theory of similarity and theory of models, are 
based on the physical relationship between the variables being required to be 
dimensionally perfect. By organization of the variables, it may be possible to 
determine a fundamental relationship between them. 

Dimensional analysis is encountered in many other branches of engineering. 
In terms of fluid mechanics, problems involving the fundamental dimensions of 
mass, length and time can be applied to all forms of fluid resistance, flow 
through pipes, through weirs and orifices. The techniques require identification 
of all the requisite variables by which a phenomenon is affected; missing or 
wrong variables will lead to incorrect conclusions. The requirement for 
forming dimensionless groups is that they should be independent of one 
another within the set, but all possible dimensionless groups outside the set can 
be formed as products of powers of the groups within the complete set. 

There has long been a preoccupation with dimensional analysis employing 
Lord Rayleigh's method of indices of 1899 and the TI theorem of Buckingham 
of 1915 (although first stated in 1892 by A. Vaschy). There are limitations with 
these theories, however, largely associated with the need to identify the vari- 
ables responsible for predicting a particular phenomenon together with the 
interpretation of the dimensionless groups. A number of techniques have s 
quently been developed that better recognize the physical features of a 
particular situation. Nonetheless, dimensional analysis reasoning remains a 
useful tool in formulating models, particularly where the laws governing a 
particular phenomenon are not known. 





e- 
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4.1 Flow through an orifice 


Deduce by the method of dimensions, an expression in dimensionless terms for 
the rate of flow Q of a liquid of density p and viscosity \t in a pipe of diameter d 
with an orifice of diameter d,, and pressure drop Ap. 


Solution 

The Rayleigh method of dimensional analysis was proposed in 1899 by British 
physicist Lord Rayleigh (1842-1919) to determine the effect of temperature on 
a gas. The method, which has found wide application in engineering, involves 
forming an equation in which A is some function of independent variables, A4, 
Ap, ... etc, which are dimensionally consistent in the form 


=Kat ae 
A SKATSAS S) 
where k is a dimensionless constant. The dimensionless groups are obtained by 


evaluating the exponentsa, b. ... etc, and grouping those variables with the same 
power. For the orifice in a pipe, the relationship is therefore assumed 


Q - kp" did pay 
where the corresponding fundamental dimensions are 
DT = MET) LPL! (ML)! (MET 
Equating indices for M, L and T, respectively 
O=a+dt+e 
3=-a+b+c-3d-e 
-12-2a-e 


The governing variables are pipe diameter and viscosity. A rearrangement of 
the indices in terms of c and e gives 
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The power relationship is therefore 


e 
s 
Di og ud d 
Q « Apid'p | — Fa 
dAp?p? o 
The constant k therefore corresponds to a form of Reynolds number and ratio of 
pipe to orifice diameter. Thus 


Ap 
P 





Q= 


This result is expected as the rate of flow depends on the velocity through the 
orifice. Since the pressure drop is related to head by 
Ap =psH 


then the rate of flow can alternatively be expressed as 
Q - kd? gH 


Note that comparing this equation with that for flow through an orifice (see 
Problem 5.4, page 122), it can be deduced that 


Can 
4 


where Cy is the discharge coefficient and a, is the area of the orifice. Much 
experimental work has been done correlating Cy with Reynolds number for 
different ratios of pipe to orifice diameter. For Reynolds numbers above 
10,000, however, Cy approaches values between 0.6 and 0.65 irrespective of 
d/dy. 
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4.2 Flow over notches 
Show by the method of dimensions that the flow of liquid Q over a notch can be 
given by the dimensionless groups 


u o 


per? 


Q-kg^H 








Ps 
where H is the head, p is the density, t the viscosity and G the surface tension of 


the liquid, g is the gravitational acceleration and k is a constant. 


Solution 

For problems involving flow over notches it is known that the experimental 
coefficients are not constant for all heads but vary with fluid density, viscosity, 
and surface tension. Assuming notches to be geometrically similar — that is, 
they have the same characteristic dimensions where the angle 6 is the same for 
all notches tested — the power relation is assumed to be 


Q - ki" gh tpo* 
then 
Q « H" ghu*po* 
The fundamental dimensions are 
DT? «(y qr? y aec !r y anr?y armr?y 
Equating the indices for M, L and T, respectively 
O=ctdte 
3=a+b -c -3d 
-|12 -2b -c -2e 


The governing variables are viscosity and surface tension. Obtaining the 
indices in terms of c and e gives 
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Therefore 


c 


Lg f 
oneni] 
pg?H? | (P8H? 





The coefficient of discharge for notches therefore varies with fluid density, 
viscosity, head and surface tension. That is, the constant k is a function of 





u o 
pei? per? 


This is expected since the flow of fluids which have a free surface is influenced 
in some ways by surface tension and viscosity. This is true for open channel 
flow, flow over notches and weirs, and multiphase fluid flow. The precise 
nature of the influence is not known and generally is only appreciable for condi- 
tions of shallow flow or low head over obstructions, Surface tension effects are 
therefore ignored for most practical purposes and are otherwise incorporated 
into the coefficient of discharge. It should be noted that surface tension effects 
do not appear in problems concerning the single-phase flow of liquids through 
pipes. 
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4.3 Scale-up of centrifugal pumps 
Show by the method of dimensions that the scale-up of centrifugal pumps can be 
based on the dimensionless groups 


P Q gH pND 
a a o XC 
pN ND? N?. rn 











where P is the power required for pumping. N is the rotational speed and D the 
diameter of the impeller, H is the head developed, Q is the flow delivered, p is 
the density and \ the viscosity of the fluid, and g is the acceleration due to 
gravity. 


Solution 
The scale-up of centrifugal pumps is a common challenge in the process 
industries. Experiments may be carried out on a laboratory scale or pilot plant 
scale as a basis for designing a full-size plant. Alternatively, data may be avail- 
able from an existing full-size plant and may need conversion to another full- 
size plant of different capacity and operating conditions. The variables of 
importance can be combined to give dimensionless groups of power, capacity 
and head coefficient. This is based on the fact that the power required for 
pumping is a function of the head developed and yolumetric flowrate as well as 
pump size in terms of impeller diameter and its rotational speed, and the proper- 
ties of the fluid in terms of density and viscosity. For the purpose of the analysis 
H is combined with g where the power relationship is assumed to be 

P - kQ" (HY N*Dp*y/ 
then 

P <Q" (gH) N°D4 pw! 
In terms of the fundamental dimensions M, L and T 

MÈT? =(BT y qr Party MLO) (MTS 
Equating the indices for M, L and T, respectively 

l=e+f 
2=3a+2b+d-3e-f 
-3=-a-2b-c-f 


The governing variables are flowrate, gravitational acceleration, head (gH) and 


104 





DIMENSIONAL ANALYSIS 


viscosity. Therefore rearranging in terms of the indices a, b and f 
e=l-f 
€—3—-a-2b-f 
d=2-3a -2b +3e +f 
=2-3a -2b +3(1-f)+f 
=5 -3a -2b -2f 
Hence 
P. oc Q" (gHY N74 729-f p5-3a-2b-3f gl ff 
or rearranged into dimensionless groups 


P_ p| Q eH pND? 
pND ` |ND? N?D? y 








where the power coefficient is a function of the capacity coefficient Co, head 
coefficient Cj, , and a form of Reynolds number, respectively. 


a f 
P «pN3p*| | [-88 
ND? J Uv? p? 


where 


LE. 
pND 

















Co = el (see Problem 9.7, page 249) 
ND^ 
and 
gH 
Cu = w2p2 
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4.4 Frictional pressure drop for turbulent flow in pipes 
Show, by the Buckingham TI method of dimensional analysis, that the frictional 
pressure drop Ap, fora fluid of density p and viscosity y flowing with a velocity 
v through a circular pipe of inside diameter d, length L and surface roughness € 
can be given by the dimensionless groups 


Bop gh Bw 
pv? d’ d` pyd 
Solution 


The Buckingham TI method of dimensional analysis proposed by Edgar Buck- 
ingham (1867—1940) in 1915 states that if, in a dimensionally consistent equa- 
tion, there are n variables in which there are contained m fundamental 
dimensions, there will be n —m dimensionless groups. Buckingham referred to 
these groups as IT), Ij, ete, such that 


F(T, Ty, ... 11, —m) =O 





In this case, the functional dependency of frictional pressure drop for the fluid 
is 

App =f(L, € p, d, p, v) 
Thus, the functional relation is 

F(Apy, L, & u, d, p. v)=0 
Any dimensionless group from these variables is therefore of the form 

Tl- Af, IP, e^ p, d*, pf, v* 
for which the fundamental dimensions of ITare therefore 

Te (uL IT? Y qo oy ub! Y wf o? y or y 
For the seven variables with three fundamental dimensions (M, L and T) it is 
expected that there will be four dimensionless groups. The procedure is to 
select repeating variables by choosing variables equal to the number of funda- 
mental dimensions (three) to provide a succession of TI dimensionless groups 


which have dimensions which contain amongst them all the fundamental 
dimensions. It is usually suitable to select p, v and d. Thus 
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TI, -p?lyPlqclap, 
Th =p va? 
TI, 2p, 48g 
Ty 2 p^*y qc 
That is 
0 - (ML? y"! qur y qo! aur!r?) 
0 -(ML? y? Lr o)? qu) 
0 - (ML? y? ur ^ y? qu) 
0 - (ML? y * (or qo arr) 


Treating each separately using the Rayleigh method gives 


n- PUE 
p. 
n2 
d 
£ 
n, == 
vd 
u 
Nya 
= pyd 
The four dimensionless groups may be therefore written as 
App Ly [bie E 
p? d' d' pvd 
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4.5 Scale model for predicting pressure drop 

in a pipeline 

A half-scale model is used to simulate the flow of a liquid hydrocarbon in a 
pipeline. Determine the expected pressure drop along the full-scale pipe if the 
average velocity of the hydrocarbon in the full-scale pipe is expected to be 1.8 
ms~!, The model uses water where the pressure drop per metre length is noted 
10 be 4 kNm7. The respective densities of water and hydrocarbon are 1000 
kgm} and 800 kgnr^, and viscosities are 1x1073 Nsm-? and 9x10 Nsmr7. 


Solution 
By dimensional analysis it can be shown that the pressure coefficient is a func- 
tion of Reynolds number of the fluid in the pipe. That is 


Arp (n 
pv? pvd 


Both the model and full-scale pipe are dynamically similar since the resistance 


to flow is due to viscosity. The Reynolds number for both model and full-scale 
pipe must therefore be the same. That is 


H Jmodet \  Jfuit scate 


Rearranging in terms of velocity through the model is therefore 


5 ) 
M Jatt scate ( d fult scale 
d model 


model =F 
B L. 





Since the Reynolds number is the same in both model and full-scale pipe, then 
so too is the pressure coefficient. That is 
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2 2 
PV^ Jmodet VP Jfull scale 


Rearranging, the pressure drop per unit length in the full-scale pipe is therefore 


(pv gy cale 

Apg — AD model em 
(PY model 

800 x 18? 


1000 x 3.27 


= 4000 x 
=1012 Nm? 


The pressure drop per unit length in the full-scale pipe is expected to be 1.012 
kNm. 
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Further problems 
(1) State Buckingham's IItheorem. 


(2) Outline the procedure used in the Rayleigh method of dimensional 
analysis. 


(3) Show that by the method of dimensional analysis the flow over a rectan- 
gular weir can be given by the dimensionless groups 


2 -,; +, 


3 E 2 
Bg?H? — (pgiH? P8H 











where B is the breadth of the weir, H is the head, p is the density, t the viscosity 
and o the surface tension of the liquid, g is the gravitational acceleration and k 
is a constant. 


(4) Rework the problem of dimensionless groups for flow over a V-notch in 
which the angle 6 is not necessarily the same for all notches to yield 











2 E 9 8 


nz pg? H? pg? 


where H is the head, p is the density, u the viscosity and o the surface tension of 
the liquid, g is the gravitational acceleration and k is a constant. 


(5) Show, by the method of dimensions, that the viscous resistance of an oil 


bearing is dependent on the linear dimensions of the bearing d, the viscosity of 
the oil u, speed of rotation N and pressure on the bearing p, and is given by 


R fe 
pna? 43 


(6) Show by dimensional analysis that a form of frictional pressure drop for 
the flow of a fluid through a pipe is 


see 





pv H 


where p is the density, t. the viscosity and v. the velocity of the fluid, and d is 
the diameter of the pipe. 
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(7) Show by dimensional analysis that a form of frictional pressure drop for 
the flow of a fluid through geometrically similar pipes can be expressed by 
Aos (prd L 
p? u d 


where p is the density, 11 the viscosity and v the velocity of the fluid, and d is the 
diameter and L the length of pipe. 


(8) Show by the method of dynamical analysis that the volumetric flow of a 
fluid through a central circular orifice located in a pipe can be expressed as 


oua e 
p 


where p is the density of the fluid, A is the area of the orifice, k is a coefficient 
which depends on the pipe and orifice dimensions and the Reynolds number, 
and Ap is the pressure drop across the orifice. 


(9) Show that the rate of flow of a liquid of kinematic viscosity v over a 90° 
V-notch can be given by 


Q Hig? 
aad £ 
Hig v 


where H is the head and g the acceleration due to gravity. 


(10) For liquids flowing along pipelines above a critical velocity, show that 
pressure drop due to friction per unit length is given by 


2 
Any iva eu) 
p as 

where p is the density and v. the velocity of the fluid, d is the diameter of the 
pipe and Re is the Reynolds number. Hence show that the friction factor f in the 


formula AfLv? [^gd for frictional head loss is a function of the Reynolds 
number. 
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Introduction 

The measurement of flow of process fluids is an essential aspect of any process 
operation, not only for plant control but also for fiscal monitoring purp: A 
wide variety of flowmeters is available and it is important to select correctly the 
flowmeter for a particular application. This requires a knowledge and compre- 
hension of the nature of the fluid to be measured and an understanding of the 
operating principles of flowmeters. 

Before the.advent of digital control systems which collect and store flow 
information, the rate of flow of fluids was usually measured by instruments 
(flowmeters) using the principle of differential pressure. Today, there is a wide 
variety of flowmeters available. For single-phase, closed pipe flow, flowmeters 
are broadly classified into those which are intrusive and those which are 
non-intrusive to the flow of the fluid. Collectively, the classifications include 
differential pressure meters, positive displacement meters, mechanical, 
acoustic and electrically heated meters. 

Differential flowmeters indirectly measure velocity, and therefore flow, of a 
fluid by measuring a differential head. Consequently, they are also known as 
head or rate meters with the main group of meters being venturi, orifice, nozzle 
and Pitot tubes — although there are others. Such meters are based on the prin- 
ciple that when liquid flows through a restriction, its velocity increases due to 
continuity. The increase in kinetic energy evolves from the reduction in pres- 
sure through the restriction and it is this relationship which allows a measure- 
ment of pressure drop to be related to the velocity and therefore flowrate. The 
relationship is derived from the Bernoulli equation and is essentially an appli- 
cation of the first law of thermodynamics to flow proce: 

Venturi meters are simple fluid flow measuring devices which operate by 
restricting the flow of fluid, thereby increasing velocity and consequently 
reducing pressure at the point of restriction. By measuring the differential pres- 
sure drop at the point of restriction the flowrate can be readily determined. The 
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device, which has no moving parts, consists of a rapidly tapered section to a 
throat and gentle downstream expansion section. This design prevents the 
phenomenon of separation and thus a permanent energy loss. For a well- 
designed venturi, the discharge coefficient should lie between a value of 0.95 
and 0.98. If the discharge coefficient is not known or is unavailable for a 
particular venturi, a value of 0.97 may be reasonably assumed. Although the 
venturi meter can be installed in any orientation, care must be taken to use the 
appropriate equation for flow derived from the Bernoulli equation. 

The orifice meter is cheaper to manufacture than the venturi meter, but has 
the disadvantage of a higher permanent energy loss. The device consists of a 
plate or diaphragm ideally manufactured from a corrosion and erosion-resistant 
material, positioned between two flanges in a pipeline. The plate is often 
centrally drilled, although eccentric and segmental diaphragms are used. The 
meter operates by increasing the velocity of the fluid as it flows through the 
restriction and measuring the corresponding differential pressure across the 
device. There is, however, a region of high turbulence behind the orifice giving 
rise to a high permanent energy loss and consequently the discharge coefficient 
is considerably less than that obtained using the venturi meter. A value of 0.6 is 
frequently used for high flowrates although the coefficient varies with the size 
of orifice relative to the pipe diameter and rate of flow. 

Of the mechanical flowmeters which function on simple fluid flow princi- 
ples, the rotameter operates using a fixed differential head but a variable flow 
area. It consists of a vertical tapered tube through which a fluid flows upwards 
and where the elevation of the float contained within the tapered tube provides 
an indication of the rate of flow. 
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5.1 Pitot tube 


A Pitot tube is used to determine the velocity of air at a point in a process venti- 
lation duct. A manometer contains a fluid of SG 0.84 and indicates a differen- 
tial head reading of 30 mm. Determine the local velocity in the duct if the 
density of air is 1.2 kgm and may be assumed constant. 


Flow — 2 





Solution 

The Pitot tube, named after the eighteenth century French engineer Henri de 
Pitot (1695-1771) who invented it, is a device used to measure the local or 
point velocity of a fluid in a pipe or duct. It can also be used to determine 
flowrate by measurements of local velocities in the cross-section of the pipe or 
duct and is particularly useful therefore where the velocity profile is irregular. 
The device operates by measuring the difference between the impact and static 
pressures in the fluid and usually consists of two concentric tubes arranged in 
parallel; one with an opening in the direction of flow, the other perpendicular to 
the flow. Applying the Bernoulli equation, where there is virtually no head loss 
due to the proximity of the two points 


v? 
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The tube is also assumed to be small with respect to pipe size; otherwise, 
corrections for disturbances and flow area reduction are required. Since the 
fluid at point 2 is stationary, p> therefore corresponds to the impact pressure. 
Rearranging 


2(p2 — Pi) 


For the differential manometer, the density of the manometric fluid, p,, , is 
considerably greater than the density of air, p. The approximate differential 
pressure is thus 


P2 -P1 =Pm8H 
=840 x g x 003 


~i 





=247 Nm? 


The velocity of the air in the duct is therefore 


2x247 
12 


= 2029 ms“! 





The velocity is found to be 20.29 ms-!. Note that the density of the air can be 
readily determined applying the ideal gas law for a known barometric pressure 
and temperature. It may also be noted that the actual velocity measured by the 
Pitot is given by 


v =C\2¢H 


where C is a coefficient which is approximately unity for large pipes and 
smooth Pitots but is appreciably less for low Reynolds number flow. 
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5.2 Pitot traverse 

Determine the rate of flow and average velocity of air in a process vent pipe of 
50 cm diameter for which local readings of velocity from a Pitot tube are 
recorded below. 


Radius r (m) 0 005 010 O15 020 40225 025 
Velocity v (ms!) 19.0. 18.6 17.7 16.3 14.2 12.9 0 


a 


20v r (ms) 
s 


5 
0 
0 oos 010  0I5 020 025 
RADIUS r (m) 
Solution 


The distribution of velocity across the process vent pipe using a Pitot tube, 
known as a Pitot traverse, can be used to determine flowrate either numerically 
or by graphical integration. The elemental flowrate for the device at a radius r, 
recording the local velocity v , , is therefore 


dQ =v, 2nrdr 
The total flow may therefore be determined graphically from a plot of 2nv , r 


versus r. Thus 


Radius r (m) 0 0.05 0.10 (0.15 020 0225 025 
Velocity v (ms-!)19.0 18.6 17.7 16.3 142 12.9 0 
2nv,r(m’s!) 0 5.8 nI 15.4 17.8 18.2 0 
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From the plot, the area under the curve is found to be 2.74 m?s-!. The average 
velocity is found from the total flow across the flow area. That is 





The average velocity is found to be 13.94 ms-!. Note that the average velocity 
is not the sum of the velocities reported divided by the number of readings. 

This procedure is applicable only to symmetrical velocity distributions. For 
unsymmetrical flow or flow in non-circular ducts, the procedure can be modi- 
fied to evaluate the total flow over the flow section. For rectangular ducts this 
involves dividing the cross-section into regular sized squares and measuring 
the local velocities at these points. Alternatively, log-linear for circular cross- 
sections or log-Tchebychev positions for both circular and rectangular cross- 
sections can reduce the computations involved but require accurate positioning. 
Numerical velocity-area integration techniques are, however, now preferred to 
graphical techniques where the accuracy of the final result depends on the 
number of measurements made. Ideally, for circular cross-sections there should 
be about 36 points with 6 on each equi-spaced radii and not less than 12 with 3 
on each of 4 equi-spaced radii. 
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5.3 Horizontal venturi meter 

A horizontal venturi meter with a discharge coefficient of 0.96 is to be used to 
measure the flowrate of water up to 0.025 m?s-! in a pipeline of internal diam- 
eter 100 mm. The meter is connected to a differential manometer containing 
mercury of SG 13.6. If the maximum allowable difference in mercury levels is 
80 cm, determine the diameter of the throat and the shortest possible overall 
length of the meter. 


——— 


Flow —> 1 


n 





Solution 

This fluid flow measuring device, first named after the Italian physicist 
Giovanni Battista Venturi (1746-1822) by Clemens Hershels in 1886, consists 
of a tapered tube which constricts flow so that the differential pressure 
produced by the flowing fluid through the throat indirectly gives a measure of 
flowrate. The rate of flow can be determined by applying the Bernoulli equation 
at some point upstream of the venturi (point 1) and at the throat (point 2). For a 
horizontal venturi 


y v 
Bi gi abs V Vos 


pg 2g pg 2g 
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From continuity for an incompressible fluid 
avı =ax2 


Substituting for v and rearranging 


2 
merata A 
pg 2g |\ a2 


Rearranging, the velocity in the pipe v, is therefore 





The actual flow through the venturi incorporates a coefficient of discharge Cy 
to allow for frictional effects and is defined as the ratio of actual to theoretical 
flow. Therefore 


Q=Cyav, 





=Cya 
p(B 


where is the ratio of pipe diameter to throat diameter and a is the flow area of 
the pipe. For a pipe diameter of 0.1 m the pipe flow area is 


) 


nx 01? 
Qe EXUP. 
4 


=785 x 10° m? 
Rearranging, the flow equation in terms of B is 


1 
a 


B= 201 = P2) 5 


fé) 
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For the mercury-filled differential manometer the differential pressure is 
P1 ~P2 =(PHe —P)gH 
= (13,600 —1000) x g x 08 
— 98885 Nm? 


Bis therefore 





b= 2x 98,885 = 
1000 x| 0025 
0.96 x 785 x 10> 
= 2087 
E 
d; 


The diameter of the throat is therefore 


That is, the throat diameter is found to be 48 mm. The dimensions of the venturi 
are important to minimize permanent energy losses. The recommended dimen- 
sions of a venturi suggest an inlet entrance cone to have an angle of 15°-20° 
with an exit cone angle of 5°-7.5°. The throat length is 0.25 to 0.5 
pipe-diameters with tapping points located between 0.25 and 0.75 
pipe-diameters upstream. In this case, the shortest possible overall length of 
venturi is therefore an enirance cone of 7.1 cm length (20^), a throat of 2.5 cm 
(0.25 pipe-diameters) and an exit cone of 19.7 cm (7.5°) giving an overall length 
0f 29.3 cm. 
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5.4 Orifice and venturi meters in parallel 

An orifice plate meter and a venturi meter are connected in parallel in a hori- 
zontal pipe of inside diameter 50 mm. The orifice has a throat diameter of 25 
mm and discharge coefficient of 0.65 while the venturi has a throat diameter of 
38 mm and discharge coefficient of 0.95. Determine the proportion of flow 
through either meter. 


e o 





Venturi 





Solution 

It is assumed that the liquid flows along a frictionless horizontal pipe and that 
losses due to fittings can be neglected. For liquid to flow across both instru- 
ments, the pressure drop across them must be the same. Due to the different 
characteristics of these two meters, the flow through the orifice should be less 
than that through the venturi since the overall pressure drop for a given flow is 
higher through an orifice than a venturi. Thus for both meters, the flowrate is 


(see Problem 5.3, page 119) 
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and ratio of pipe to venturi throat area in terms of diameter is 


E 


Since the differential pressure across both meters is the same, then 


2 
os (a Y a) {e aY 
Cato } |do Cat) J Uv 


That is 


Gs] » sen] ^ 


which reduces to 
Q, -40, 


From continuity, the total flow Q is the sum of the flow through the meters. That 
is z 


Q =Q, +0, 
=Q, +40, 
-50, 


Thus 20% of the flow passes through the orifice meter while 80% of the flow 
passes through the venturi, Note that for the same throat size, the difference in 
volumetric flowrate is entirely due to the permanent frictional losses imposed 
by both meters. 
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5.5 Venturi meter calibration by tracer dilution 

A venturi meter with a throat diameter of 6 cm is used to measure the flowrate 
of water of density 1000 kgm-? along a horizontal pipeline with an inside diam- 
eter of 100 mm. The flowmeter is calibrated using a solution of salt water with a 
salt concentration of 20 gt"! and added continuously upstream of the meter at a 
rate of 3 litres per minute. Determine the discharge coefficient of the meter if a 
pressure drop across the throat of 3.7 kNm- is recorded and a diluted salt 
concentration of 0.126 gl-! is found by analysis downstream. 


Salt Sample point 


Flow Q, concentration C, 
y. Concentration C, 
Water 


flow —e —e Flow 0+0, 


9; 





@) 


Solution 

The calibration of flowmeters by the method of tracer dilution involves the 
addition of extraneous material whose presence can be quantitatively deter- 
mined by an analytical technique. The extraneous material may or may not be 
already present as an impurity. In this case, the flowrate of water through the 
meter is related to diluted salt concentration. From a material balance on the salt 


Q,C, 7 (Q; *Q,)C, 


whereQ; andQ,, are the rates of flow of water and salt solution, andC,, andC,, 
are the added upstream and diluted downstream salt concentrations. Rear- 
ranging, the flowrate of water is therefore 


Ca -Co 
Qi «o Es 


.. 0003 s 20 -0126 
60 0126 











=789 x10 ms! 


124 





FLOW MEASUREMENT BY DIFFERENTIAL HEAD 


The total flow Q through the meter during the calibration is therefore 
Q =Q; +Qa 


=789x10 +5x105 
27895 x 10 ? ms! 
The rate of flow through a horizontal venturi meter is given by 


2Ap 
p(B 


where D is the ratio of pipe to throat diameter 


Q=Cya 





) 


= 167 


Rearranging in terms of discharge coefficient 





_ 7895 x 10? 
zx0P . [ 2x3700 
4 1000 x (1.674 





= 0.962 


The coefficient of discharge is found to be 0.962. 
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5.6 Differential pressure across a vertical venturi meter 
A process liquid of density 850 kgm~ flows upward at a rate of 0.056 n?s-! 
through a vertical venturi meter which has an inlet diameter of 200 mm and 
throat diameter of 100 mm, with discharge coefficient of 0.98. Determine the 
difference in reading of two pressure gauges located at the respective tapping 
points a vertical distance of 30 cm apart. 








Solution 
Applying the Bernoulli equation between the upstream position (1) and the 
throat (2) 


2 2 
v v, 
AEA E A 
pg 2g ps 2g 
From continuity for an incompressible fluid 


aw; =a 
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Rearranging in terms of differential pressure 


7 


2 

v a 

Pi-P2 = pa} att aD 522 2. 
2¢|\a 


The velocity is related to flowrate by 
Q =Cyav 


and for the circular cross-section 


E 


therefore 
2 4 
p}( 2 d, 
-p»--b|c—l||z| -1i[t2&(G2 -z) 
Pı -P2 ds S 8(z2 —z| 
2 
4 
2 
| a Es -1 2x g x03 
2 MES 01 
ù 4 
= 23,594 Nro? 


The differential pressure is found to be 23.6 kNm™?. Note that a differential 
manometer would not record a difference in level of manometric fluid when 
there is no flow through the venturi whereas two independent pressure gauges 
would record a static pressure difference equal to 


Pi -pa2 7petzo -zi) 
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5.7 Flow measurement by orifice meter in a vertical pipe 
Oil of density 860 kgm™ flows up a vertical pipe section of diameter 225 mm. A 
manometer filled with fluid of density 1075 kgm is used to measure the pres- 
sure drop across an orifice plate with a throat diameter of 75 mm. Determine 
the flowrate of oil if the deflection of the manometer fluid is 0.5 m. Assume a 
discharge coefficient of 0.659 for the orifice. 


n 


— Orifice 





Flow 


Datum 
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Solution 
Applying the Bernoulli equation between points 1 and 2 
2 2 
v v. 
Le ee E 
pg 28 pg 2g 


From continuity for an incompressible fluid 
aw =a 


Therefore substituting for v and rearranging 


3 
pry 
Di ~P2 + Paley ue fae 


For the manometer 
Py - P2 7pg(za — 21) + (Pm - 8H 
The theoretical velocity through the pipe is therefore 





Note that this expression does not contain terms in z. The velocity and therefore 
flowrate is independent of the orientation of the pipe. The actual flowrate is 
then 


Q - C;ay, 


2 
= 0.659 x eee x0523 


— 0014 m?s7! 





The rate of flow is found to be 0.014 m? 
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5.8 Variable area flowmeter 

A rotameter used to measure the flow of water consists of a float with a mass of 
30 g set in a tapered glass tube 20 cm in length. The tube has an internal diam- 
eter of 22 mm at its base and 30 mm at its top end. Determine the rate of flow 
when the float is at mid-height in the tube. The rotameter has a coefficient of 
discharge of 0.6 and the density of the float is 5100 kgm. 





Flow 


Solution 

Often referred to as a rotameter, the variable area flowmeter is used to measure 
the rate of flow of a fluid by virtue of the elevation of a solid float within a 
vertical tapered tube. The tube is transparent and is usually made of glass while 
the float is made of metal, ceramic or plastic and is usually ‘bomb’ shaped, in 
that it has a cylindrical body with a cone-shaped bottom and short flat top piece. 
Some floats are grooved which encourages the float to spin thereby improving 
stability due to a gyroscopic effect. This device has effectively a fixed perma- 
nent pressure drop and a variable flow area. As the fluid flowrate is increased, 
the float moves up the tube until an equilibrium position of forces is reached. 
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From Archimedes' principle, this is a position where the upthrust is equal to the 
weight of fluid displaced. That is 


ApA; -pyVrg -pVyg 
where Ap is the cross-sectional area of float, py is the density of the float, p is the 
density of fluid and V, is the volume of float. Rearranging in terms of Ap then 


ap = ser 79) 
Ar 


where Ap is related to flowrate by the Bernoulli equation. At the elevation of the 
float 


2 2 
BoE ME ox I Me 


Pg 28 pg 2g 


where v, and v,, are the velocities of the fluid in the tube and annulus, respec- 
tively. Therefore 
Ap bo; 9D 
2 
Since the velocity of the fluid through the annulus is significantly greater than 
that in the tube, this approximates to 
2 


m urs 
ee 


Combining both equations for Ap, the velocity of the fluid through the annulus 
is therefore 


. BVgstpr - p) 
a pA; 


Introducing a discharge coefficient C; to allow for losses due to friction, the 
actual mass flowrate is therefore given by 


m - Capa,v, 


2Vrg(py —p) 
-Cqpa, E B 
nar 


2Vreppr —p) 

3 1 SP y 

s ja, eno p 
lr 
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If the density of the fluid is not influenced by temperature or composition, then 
the square root term remains constant. Further, if the discharge coefficient does 
not vary greatly, then there is an almost linear variation between mass flowrate 
and the area of the annulus. To relate the area of the annulus to the float position 
above the bottom of the tube, the float is assumed to sit perfectly in the bottom 
of the tube of internal diameter d; when there is no flow. From geometry, the 
internal diameter of the tube d for any float position above the bottom H, for an 
angle of taper 6, is therefore 


d=d, EL 


for which the area of the annulus, a, , around the float is 


aa Tw -dp 


x e 2 
=" (a, +2 tan®) —a 
( $ mS) J 


=2| a2 +4dyH wan ® + 4H? tan( 9 )| a2 
4 2 2 


This approximates to 
e 
a, 7 Td,H tan — 
2 
The mass flowrate is therefore 


m=Cynd,H un? 





=C4H tan 5 8V,gp(py —p)r 


That is, the mass flowrate is approximately linear with float position in the tube. 
In practice, a scale is marked on the tube and the rotameter is supplied with a 
calibration curve for a particular fluid, temperature and float. In this case, the 
float has a volume 
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Vass 
P 
. 003 
5100 
= 588 x 10° m? 


The half-angle of taper is related to top and bottom diameters separated by a 
height H Tube 


tan 8 = db cdi 
2 2H 
_ 003-022 
2x02 
= 002 


The mass flowrate through the rotameter for a float at the mid-point in the tube 
(10 cm) is therefore 


m =06 x 01x 002 x 8 x 588 x 10° * g x 1000 x (5100 — 1000) x 2 
= 0.0925 kgs! 


corresponding to a volumetric flow of 5.55 litres per minute. 
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5.9 Rotameter calibration by venturi meter 

A horizontal venturi meter with a throat diameter of 2.5 cm is used to calibrate 
a vertical rotameter in a pipeline carrying water. The pipeline has an inside 
diameter of 5 cm. The rotameter has a cone angle of 2° and uses a float with a 
volume of 100 cm? and density 8000 kgm~. Determine the flowrate of water 
and the discharge coefficient for the rotameter if a differential pressure reading 
of 1.2 kNm~ is recorded for the venturi when the float is at an elevation of 15 


cm. The discharge coefficient for the venturi is 0.97. 


Flow Q 





Solution 
The flow of water through the horizontal venturi is 


Q=Cya 





p(p 
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Therefore 


Q - 097 1X 005* x [.2x1200 
4 


1000x 2* -1) 
7762x107 ms! 


The flow through the rotameter is given by 
CH. 6 
Q- Moe 2 BV spp, -pm 


Rearranging, the discharge coefficient is therefore 
— —PÀ 
H8 
tan — 48V, = 
p "8p y —p)n 


762x107 


z 
915 xn? x y8 x 1x 10™ x g x 1000 x (8000 —1000) x x 


1000 
-070 


Ca = 





The coefficient for the rotameter is found to be 0.70. 
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Further problems 

(1). Water flows down a vertical tapering pipe 2 m long. The top of the pipe 
has a diameter of 10 cm and the diameter of the bottom of the pipe is 5 cm. 
Determine the difference of pressure between the top and the bottom ends of the 
pipe if the flow rate is 1 m5min-!. 


Answer: 1766 Nm? 


(2) The velocity of water in a pipe with a bore of 250 mm is measured with a 
Pitot tube. The difference in head at the centre of the pipe is found to be 10 cm 
of water. Determine the flowrate of water per minute if the average velocity of 
water is two-thirds the velocity at the centre. The coefficient of the Pitot tube 
may be taken as unity. 


Answer: 0.0458 m 





(3) A liquid of density 800 kgm-? flows through a horizontal pipe with an 
inside diameter of 150 mm under a pressure of 400 kNm-?. Assuming no 
losses, determine the flow when the pressure at a 75 mm diameter reduction is 
200 kNm2. 





Answer: 0.102 m? 


(4) Starting with the Bernoulli equation, show that the volumetric flowrate of 
an incompressible fluid through a horizontal venturi meter where the throat and 
upstream positions are connected by a U-tube containing a fluid of density p j, 

can be given by 





where Cy is the discharge coefficient, a is the pipe flow area, p is the fluid 
density, D is the ratio of the pipe diameter to throat diameter, g is the gravita- 
tional acceleration and H is the difference in levels of the fluid in the U-tube. 


(5) A venturi meter of inlet diameter 10 em and throat diameter 5 cm is 
installed in a vertical pipe through which water flows upwards. The meter is 
calibrated where it is noted that for a flow rate of 11.5 litres per second, the 
difference in readings of the pressure gauges connected to the inlet and throat is 
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20 kNm~. If the difference in the height of the two tapping points is 30 cm, 
determine the discharge coefficient for the venturi meter. 


Answer: 0.971 


(6) A venturi meter is installed in a horizontal pipe with an inside diameter of 
100 mm carrying an organic solvent of density 1200 kgm~. The only equip- 
ment available for measuring pressure differences is a manometer containing 
mercury of SG 13.6, with a maximum allowable difference of levels of 80 cm. 
Assuming a discharge coefficient is 0.97, determine the maximum flowrate of 
solvent which can be measured. 


Answer: 0.0251 m3s-! 


(7) The flow ina pipeline of diameter d is measured by an orifice meter with 
an orifice diameter of 0.5d. The discharge coefficient is 0.62. Calculate the 
throat diameter of a horizontal venturi meter which would give the same pres- 
sure difference at the same rate of flow in the pipe. The discharge coefficient for 
the venturi meter may be taken as 0.96. 


Answer: 0:405d 


(8) Discuss briefly the merits of both the orifice plate meter and the venturi 
meter as differential pressure meters. 


(9) Describe the types of orifice plate which are available and their 
application. 


(10) A horizontal venturi meter with a 5 cm diameter throat is used to measure 
the flow of slightly salt water in a pipe of inside diameter 10 cm. The meter is 
calibrated by adding 2 litres per minute of 0.5 molar sodium chloride solution 
upstream of the meter and analysing a sample of water downstream from the 
meter. Before the addition of the salt solution, a I litre sample required 10 cm? 
of 0.1 molar silver nitrate solution in a titration. After addition, a 1 litre sample 
required 24.4 cm? of 0.1 molar silver nitrate. Determine the discharge coeffi- 
cient of the meter for a differential pressure across the venturi of 17.4 kNm^? . 





Answer: 0.96 


(11) Identify suitable types of tracer substances with their appropriate detec- 
tion systems which can be used to calibrate flowmeters. 
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(12). Anorganic liquid of density 980 kgm flows with an average velocity of 
3 ms-! along a 75 mm diameter horizontal pipe. A restriction is placed in the 
pipe in which there is a 50 mm diameter opening for the water to flow. The pres- 
sure of the water in the 75 mm pipe is 150 kNm-?. Determine the pressure at the 
restriction. 


Answer: 144.2 kNm? 


(13) The difference in head registered between two limbs of a mercury gauge, 
with water above the mercury and connected to a horizontal venturi meter, is 20 
cm. If the venturi has a pipe and throat diameter of 15 cm and 7.5 cm, respec- 
tively, determine the discharge through the meter assuming a coefficient of 
discharge of 0.97. 


Answer: 0.031 m?s-! 


(14). The position of a float, of density 3000 kgm, in a rotameter for a partic- 
ular flow is 15 cm above the bottom of the tapered tube. If the float is replaced 
with an identical size of float but with a density of 2500 kgm™, then the posi- 
tion of the float in the tube is increased by 2.5 cm. Determine the density of the 
fluid. 


Answer: 1115 kgm^? 


(15) Highlight the advantages and disadvantages of using rotameters for the 
measurement of fluid flow. 


(16) A uniform pipeline of circular cross-section is to carry water at a rate of 
0.5 m3s-!. If an orifice plate is to be used to monitor the flowrate in the line, 
determine the differential pressure reading across the plate at the design 
flowrate. The pipe has an internal diameter of 45 cm and the orifice has a 
concentric hole of 30 cm. If the differential pressure is to be measured using a 
mercury-filled manometer, determine the difference in level between the two 
legs. The coefficient of discharge for the plate may be taken as 0.6. 








Answer: 55.8 kNm-?, 45.1 cm 


(17) An undergraduate laboratory experiment involves a venturi meter which 
is used to measure the rate of flow of water through a horizontal pipe with an 
internal diameter of 5 cm. A calibration is provided for the meter and is given 
by 


Q =0096 VH 
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where Q is the flowrate in litres per second and H is the manometric head in a 
water/mercury manometer measured in millimetres. Determine the throat 
diameter of the venturi if the coefficient of discharge may be assumed to be 
0.97. 


Answer: 16.9 mm 


(18) A Pitot tube is used to measure the velocity of water at a point in an open 
channel. If the Pitot is connected to an open air-filled inverted manometer, 
determine the velocity of the water for a difference in liquid levels in the 
manometer legs of 10 mm. 


Answer: 0.44 ms-! 


(19) A Pitot traverse is used to record the local velocities of air across the 
entrance of a circular duct of diameter 91 cm. The Pitot is connected to a water- 
filled differential manometer and measurements of the difference in levels 
between the legs are given below: 


Distance from pipe wall, cm 0 5 10 20 30 40 45.5 
Difference in levels, mm 0 7 10 12 14 15 15 


If the air has a static pressure of 101.3 kNm~ and temperature of 20°C, deter- 
mine the rate of flow and the average velocity. The mean molecular mass of air 
is 29 kgkmol-! and the Universal gas constant is 8.314 kJkmol-! K-!. 


Answer: 7.5 m3s~!, 11.5 ms“! 


(20) Outline the advantages and disadvantages of Pitot tubes for the measure- 
ment of fluid velocity and flow. 


(21) A venturi meter measures the flow of water in a 100 mm inside diameter 
horizontal pipe. The difference in head between the entrance and the throat of 
the meter is measured by a U-tube, containing mercury (SG 13.6) with the 
space above the mercury in each limb being filled with water. Determine the 
diameter of the throat of the meter such that the difference in the levels of 
mercury shall not exceed 300 mm when the quantity of water flowing in the 
pipe is 10 kgs~!. Assume the discharge coefficient is 0.97. 


Answer: 39 mm 


139 


























Tank drainage 
and variable 
head flow 


Introduction 

The rate at which fluids freely discharge from tanks and vessels through 
orifices and connecting pipes is dependent on the pressure or head within the 
tank and on frictional resistance. The shape, size and form of the orifice through 
which the fluid discharges, and the length and diameter of pipe and fluid prop- 
erties, also influence the rate of discharge due to the effects of frictional resis- 
tance. As the fluid passes through the orifice issuing as a free jet, it contracts in 
area, thus further reducing the rate of discharge. The contraction is caused by 
the liquid in the vessel in the vicinity of the orifice having a motion perpen- 
dicular to that of the jet and exerting a lateral force. The section of the jet at 
which the streamlines become parallel is known as the vena contracta, and 
there is no further contraction beyond this point. Owing to the reduction in both 
velocity and flow area of the jet, the actual rate of flow is much less than the 
theoretical prediction: the relation between them is known as the coefficient of 
discharge. 

The simplest method of determining the coefficient of discharge is to 
measure the quantity of fluid (liquid) discharged for a constant head and divide 
by the theoretical discharge. In terms of energy conversion for a freely draining 
tank, the potential energy is converted to kinetic energy or, in the head form, the 
static head is thus converted to velocity head. It is therefore deduced that the 
theoretical velocity of the jet is proportional to the square root of the head, 
described in what is called the Torricelli equation after the Italian mathemati- 
cian Evangelista Torricelli (1608-1647). The difference between the actual and 
theoretical velocities, known as the coefficient of velocity, is due to friction at 
the orifice and is small for sharp-edged orifices. The coefficient of velocity 
varies depending on the size and shape of the orifice as well as the applied head. 
Nevertheless, typical values are in the order of 0.97 and may be found experi- 
mentally for a vertically-mounted orifice in the side of a tank by measuring the 
horizontal and vertical co-ordinates of the issuing jet’s trajectory. 
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For tank drainage problems which do not involve simultaneous inflow to 
maintain a constant head above the orifice, there is a continuous loss of capacity 
and consequently of head. The rate of discharge is therefore not constant ‘but 
variable, with the head being dependent on the geometry of the tank. The time 
taken for tank drainage can be considered by applying an unsteady state mass 
balance over the tank and equating the rate of discharge with the loss of 
capacity from the tank. Analytical solutions can be deduced for most tank 
configurations in terms of tank geometry, for simultaneous inflow, flow into 
adjacent vessels through openings, submerged orifices (sluices), and for flow 
through connecting pipes with laminar flow. Problems in tank drainage which 
involve the flow through connecting pipes with turbulent flow are not, 
however, mathematically straightforward. 

In simple treatments of tank drainage problems, it is necessary to employ 
pertinent assumptions, These include a well-ventilated tank where the applied 
pressure (usually atmospheric) is the same at both the free surface of the liquid 
in the tank and at the jet, where there is no free vortex formation and where the 
discharge coefficient is constant. In practice, however, each of these has an 
influence and should be allowed for where appropriate by more detailed 
calculations. 
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6.1 Orifice flow under constant head 

Awater cooling tower receives 5000 m? of water per day. If a water distribution 
system is to be designed to cover the area of the tower using 1 cm diameter 
orifices in the base of the ducting in which the water will be at a depth of 20 cm, 


determine the number of orifices required. A discharge coefficient of 0.6 is 
assumed for the orifices. 


Solution 
To determine the rate at which the liquid discharges through the orifices, it is 
necessary to determine the velocity of the liquid in the free jets. Assuming that 


there is no pressure change across the free jet, the static head is directly 
converted to velocity head as 


2 





Rearranging, the theoretical velocity through the orifice, v, is therefore 
v 2 42eH 


and is known as the Torricelli equation. The actual velocity v ac , however, is 
less than this theoretical prediction due to permanent and irreversible energy 
losses and can be found by introducing the coefficient of velocity C, , which 
represents the ratio of the actual to theoretical jet velocity. That is 


y 
= Mace 
CG = 
T 


As liquid emerges from the orifice, the streamlines converge to form a vena 
contracta just beyond the orifice, at a distance of half the orifice diameter. The 
cross-sectional area of the vena contracta can be found by introducing a coeffi- 
cient of contraction C, , and is the ratio of actual to theoretical flow area. The 


coefficient of discharge Cy is therefore defined as the product of C, and C, 
where 


Qaa C, Cao 28H. 
7 Cua, 28H. 
The total flow through all the orifices is therefore 
Q - nCqa, 28H. 


where n is the number of orifices. Rearranging 
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E Q 
CgagV2gH 
"2000. 
3600 x 24 


Lo 360x224 
2 
06x r0 x fix g x02 


= 620 





‘A total of 620 orifices is required. Note that at the vena contracta the stream- 
lines are no longer converging as the flow passes through the orifices, but are 
parallel, The internal pressure in the jet will have decreased from its upstream 
pressure to atmospheric which is imposed on the free jet. At this point, the 
velocity across the jet is essentially uniform and equal to the total head. 


144 





TANK DRAINAGE AND VARIABLE HEAD FLOW 


6.2 Coefficient of velocity 

To determine the coefficient of velocity of a small circular sharp-edged orifice 
in the side of a vessel, the horizontal and vertical co-ordinates of the trajectory 
of the jet were measured for a head of 20 cm. The horizontal co-ordinate from 
the vena contracta was found to be 86 cm whilst the vertical co-ordinate was 96 
cm. Determine the coefficient of velocity. 





Solution 

The coefficient of velocity for sharp-edged orifices located in the side of a tank 
can be readily determined from the trajectory of the issuing free jet, in which 
the horizontal distance travelled by the jet measured from the vena contracta is 


X =V act 


through a vertical distance 
er 
zig 
ne 8 
Eliminating time f, the actual velocity of the jet is 


Vact = 
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The theoretical velocity is given by Torricelli’s equation as 
» - igi 
The coefficient of velocity is therefore 


V act 











24yH 
086 


“2x J096 x 02 


= 0981 


The coefficient of velocity for the orifice is found to be 0.981. Note that the 
relationship between velocity and head which was first developed by Torricelli 
in around 1645 is a simple treatment of the Bernoulli equation in which static 
head is converted to velocity head. Since the static head varies uniformly with 
depth in the vessel, it is possible to predict the trajectory of the free jet from 
sharp-edged orifices for other elevations on the side of the vessel. 

The velocity of the discharging jet sumed to be uniform at the vena 
contracta (see Problem 6.1, page 143). This is true for sharp-edged orifices but 
not for rounded orifices. Rounded orifices have the effect of slightly reducing 
the coefficient of discharge. Significantly, the coefficient of contraction is more 
marked with sharp-edged orifices. Precise measurements of the area of the vena 
contracta, however, are not always possible. Better estimates can usually be 
obtained from estimates of the coefficient of discharge and coefficient of 
velocity such that C, 2C, /C, . 
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6.3 Drainage from tank with uniform cross-section 

A cylindrical tank of diameter 1 m mounted on its axis contains a liquid which 
drains through a 2 cm diameter hole in the base of the tank. If the tank origi- 
nally contains 1000 litres, determine the time taken for total drainage. The 
discharge coefficient for the hole may be taken as 0.6. 








Solution 

Consider a tank of uniform cross-section A, , through an orifice of area ag, 
located in the base of the tank. An unsteady state mass balance on the tank for 
no inflow of liquid relates the rate of flow from the tank through the orifice to 
the change of capacity of the tank 


Cyd, 29H — -A S. 


dr 
Rearranging, the time taken for tank drainage from a head H to H5 can be 
found by integration 


i My, 
fa -— m [nan 
o0 — Cano 28 h, 


Completing the integration 


fac. Bly 
Caa \8 


-H 





uc 
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For total drainage, H =0. Therefore 


A E 
Cady V 8 


The initial liquid level above the orifice, H ;, is obtained from the volume of the 
tank V, , where d, is the diameter of the tank 














us s 
t 
_ 4x1 
axe 
=127m 
Therefore 
nx? 
P: 4 x xia 
ogx "X007 V 8 
4 
-2120s 


Total drainage is found to take 35 minutes and 20 seconds. 
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6.4 Tank drainage with hemispherical cross-section 
A hemispherical tank of 4 m diameter contains a liquid and is emptied through 
a hole of diameter 5 cm. If the discharge coefficient for the hole is 0.6, deter- 
mine the time required to drain the tank from an initial depth of 1.5 m. 











Solution 


Applying an unsteady mass state balance on the tank, the rate of discharge from 
the hole is equal to the change in capacity in the tank 


Caas JIH =-A H 


The depth of liquid in the tank does not change uniformly with cross-sectional 
area of the tank. The depth H is related to the surface area of the liquid in the 
tank where 


(see Problem 6.3, page 147) 


A =m? 


The distance x is related to depth H by Pythagoras 


R =x? +y? 
where 
R=y+H 
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Therefore 
X =R -q-Hy 
=2RH - H? 
The balance therefore becomes 
H, _H2 
-n (QRH -H Jan 


t 
je 2 ——— : 
0 Cade V28 A, H? 


ao i 3 
=— 7 — | QRH? -H?)AH 


C4aoN28 jy, 
Therefore on integration 


n 4 RH? -H3 E ZU 
t8 —E gn? -u2)-2ui - 
Lm DOT Ws 


To empty the tank, H3 =0. Therefore 


x (45,3 _2 J 
peat (Sey? =2 92 
mm Vs 
— 
og x X005. y Jag 


3 i 
4 2x15? -2 x15? 
3 5 
= 2286 s 


The time to drain the tank is found to be approximately 38 minutes. Note that if 
the vessel were initially full and is completely emptied, then 


H; =R 
H, =0 


to give 


5 
2o dMnR 
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6.5 Tank drainage with cylindrical cross-section 

Derive an expression for the time to discharge liquid from a depth H in a cylin- 
drical tank of radius R and length L positioned on its side through an orifice of 
area a, located on the underside of the tank. 





f row ou 


Solution 
The time taken for the tank to drain can be determined from 


Caa, gH =-A “ 
Ir 


The area of the liquid surface varies with depth where 
A -2xL 


(see Problem 6.3, page 147) 


The distance x is be related to liquid depth by Pythagoras 
R= ty? 

where 
R=y+H 

Therefore 


x =V2RH — H? 
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The balance equation is therefore 
2,dH 
=-2V2RH -H^L—— 
Ca, 428H E 
Rearranging, the total time for discharge is found by integrating 


t H, 1 
fur = 4 _ f on - uy an 
0 8 H, 





to give 


à 4 
AL (aR Hy)? -QR-H,)*) 


= 
3C yay 28 


Note that if the vessel was half full and was completely emptied through the 
orifice then 


H,=R 
H, =0 


then the time for drainage would be 


3 
4LR?* 


3 yay Pe 
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6.6 Drainage between two reservoirs 

A reservoir beneath a small forced convection water cooling tower is 4 m long 
and 2 m wide. Before recirculation, the water flows into a smaller adjacent 
reservoir of dimensions 2 m by 2 m by way of a submerged opening with a flow 
area of 100 cm? and discharge coefficient of 0.62. During normal operation the 
level in the larger reservoir is 30 cm above that in the smaller reservoir. Deter- 
mine the rate of flow between the reservoirs during normal operation and if the 
recirculation is halted, determine the time taken to reduce the difference in 
levels to 10 cm. 








Solution 
For steady operation with a continuous flow between reservoirs the rate of flow 
through the submerged opening is given by 


QO =Cyay2gH 
= 062 x 001 x 42g x 03 
- 0015 m^! 


Thus, the flow under steady conditions is found to be 0.015 m3s~!, When the 
recirculation is halted the drainage from the larger tank to the smaller tank 
proceeds and the difference in levels H changes by increments of dH while in 
the larger and smaller tanks by dh, and dhp, respectively 


H -dH - H —dh, -dh 
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The incremental loss in capacity of the larger tank is equal to the incremental 
gain in capacity in the smaller tank. That is 


Adh; = Azdhz 
The incremental difference in levels is therefore 
A 
dH -dh 1 —- 
^ 


The change in capacity of the larger tank is given by the unsteady state equation 
(see Problem 6.3, page 147) 


Cua gH =-A, A 


Substituting for dh, 


-Á 
Ca 28H gp (dH) 
A, \ de 
l+ 
A2 


Rearranging 
T "m H. zi 
far = AL f H aH 


0 coi yi y 
A2 
Integrating gives 
1 1 
2A? - H3) 
A 
Cyay2g| 1+ — 
ehe) 


l L 
2x4x2x(03? -012) 


Tmonne Bz fi x2) 
062x001x 42g x|l* —— 
ve ( ous] 


t= 


=45s 


The time taken is found to be 45 seconds. 
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6.7 Tank inflow with simultaneous outflow 

Derive an expression for the time taken for the liquid in a tank to reach a new. 
level if the liquid drains from an orifice while there is a constant flow of liquid 
into the tank. 








| Flow out 


Solution 
The unsteady state mass balance on the tank is 


dH 


Q =C4ao 28H +A E (see Problem 6.3, page 147) 
Ir 


whereQ is the rate of flow into the tank and dH /dr is the rate of change in liquid 
level. For convenience, k is used to group the constants 


k=Cyay 2g 


Therefore 
okla sai 
dt 


To determine the head of liquid in the tank at any given time, it is therefore 


necessary to integrate between the limits H ; and H5 
1 Hy 
for=aj a 


0 H, Q —kH? 
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To simplify this integration the substitution is used 
1 
u=Q -kH? 
Rearranging 
H= Q -u 
k 


and differentiating with respect to u 
_ -UQ -u)du 


dH z 


The integration becomes 


ju porta 


k? 
ert) 


Completing the integration gives 





sd 


i 
Q-kH; 
1 


Q -kH} 


"em 
p= Algot *KHi-HD 
n 


or in full 


-24 k vef -Cua, 2a Hs 


Po 
ET Ss £ dov 2 Cla JDg(H2 - H2) 
(Cqa, A28 Y Q -Caa, 428 Hy | ° um 
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6.8 Instantaneous tank discharge 

A concrete tank is 15 m long by 10 m wide and its sides are vertical. Water 
enters the tank at a rate of 200 litres per second and is discharged from a sluice, 
the centre line of which is 50 cm above the bottom of the tank. When the depth of 
water in the tank is 2 m, the instantaneous rate of discharge is observed to be 
400 litres per second. Determine the time for the level in the tank to fall 1 m. 


Solution 

For a tank of uniform cross-section which receives a steady flow of liquid at a 
rateQ and is allowed to drain freely through an opening or sluice, the time taken 
to alter the level can be shown (see Problem 6.7, page 155) to be 


ki 1 1 
—— Hj - Hi) 
p kH? 

Q-kHi 





The discharge coefficient and the area of the sluice are not provided. However, 
the rate of flow through the sluice Q, is related to head by 


L 
9, - ur 
or 
, 


k-Q,H ? 


When the level was 2.0 m the flow was 0.4 m?s-!, That is a depth of liquid 
above the centre line of the sluice, H , of 1.5 m. Thus 


m 
k-04x15? 
5 

-0326 mis! 


The final depth, H ? , is 0.5 m above the centre line of the sluice. The time taken 
for the level to fall is therefore 
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i 1 
*KHi-Hp 
i 1 1 
= 2X15%10 |O. jog |22032605" |; 0326x (05? 15%) 
0326 02 -0326x1 





-1536s 


The time is found to be approximately 25!; minutes. 
Note that although the area of the sluice and discharge coefficient are not 
provided, the constant k is used to substitute 


k=Cyay2g 


Assuming a reasonable discharge coefficient for the sluice of 0.62, the area of 
the sluice a can be found by rearranging 


k 
= 
Ca 2e 
_ 0326 

062 x 42g 
-0119 m? 


This would correspond to a sluice of square cross-section with sides measuring 
34 cm. Care should be taken, however, when considering discharges through 
large orifices where the head producing flow may be substantially less at the top. 
than at the bottom. 
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6.9 Instantaneous tank inflow with outflow 

A tank of uniform cross-section is provided with a circular orifice 50 mm in 
diameter in the bottom. Water flows into the tank at a uniform rate and is 
discharged through the orifice. It is noted that it takes 90 seconds for the head 
in the tank to rise from 60 cm to 70 cm and 120 seconds for it to rise from 120 
cm to 125 cm. Determine the rate of inflow and the cross-sectional area of the 
tank assuming a discharge coefficient of 0.62 for the orifice. 





TIME t 


Solution 
An unsteady state mass balance on the tank which receives a flow of water at a 
constant rate Q is 


Q -Cya, gH. A z (see Problem 6.7, page 155) 





Rearranging 
dH .Q Cua, 28H. 
d A A 
1 
=2 He 
A A 
where 
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k =Cyay\2g 


2 
= 062% OO x Pe 


s 
2539x107 mis! 


For an average head of 0.65 m, the change in level is 0.1 m over 110 seconds. 
That is 
01 dH 


90 dr 


1 
.Q 539x10? x065* 
(A A 
For an average head of 1.225 m, the change in level is 0.05 m taking 120 
seconds. That is 
005 dH 
120 dt 





i 
Q 5.39% 10 x 1225? 
A A 


Solving the simultaneous equations 
A =2.34 m 
Q = 694x107 mis! 


The area of the tank and flowrate are found to be 2.34 m? and 6.94 litres per 
second, respectively. 
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6.10 Tank discharge through a horizontal pipe with 
laminar flow 

A viscous Newtonian liquid of density 1100 kgm~ and viscosity 0.08 Nsm~ is 
fed to a process from a vessel of diameter 1.2 m through a 3 m length of hori- 
zontal pipe with an inside diameter of 25 mm attached near to the base of the 
vessel. If the initial level of liquid in the vessel is 1.5 m above the pipe and the 
flow through the pipe is laminar, determine the time to feed 0.75 m? of the liquid 
10 the process. 








Solution 
Neglecting entrance and exit losses, an unsteady mass state balance for the 
vessel is 
ieee 
dr 


If the head in the vessel is used to overcome the viscous resistance of the pipe, 
the average velocity of the flowing liquid is therefore given by 





m Apd? 
Sa (see Problem 3.7, page 73) 
_ pgHd? 
32L 
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Therefore 


apeHd" aH 
3L dr 


The total time taken for the level to fall from a head H ; to H is therefore found 
from 





where the final level in the vessel, H5, is found from the total volume 
discharged to the process 





av 
Ho =Hy-— 
] 
=15 ~4%075 
mx 12 
= 084 m 


The time taken is therefore 


mx 12? 
-32 x 0.08 x 3x 4 084 
be 1s 
x 0025" 5 1100. g x 0025? 
= 15218 


The time is found to be approximately 25! minutes. 
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Further problems 
(1) State the assumptions upon which simple tank drainage problems are 
based. 


(2) A cylindrical vessel mounted vertically on its axis has a cross-sectional 
area of 1.2 m? and contains water which is allowed to discharge freely through 
an orifice, with a cross-sectional area of 100 cm?, positioned in the side of the 
tank 4 m below the surface. If the tank is open to atmosphere, determine the time 
taken for the level to fall by 2 m. The discharge coefficient may be taken as 0.6. 


Answer: 53 s 


(3) A rectangular orifice in the side of a tank is 1 m broad and 50 cm deep. 
The level of water in the tank is 50 cm above the top edge of the orifice. Deter- 
mine the flow through the orifice if the coefficient of discharge is 0.62. Note 
that the velocity through the orifice may not be regarded as constant as the 
variation in head at different depths of the orifice will be considerable. 


Answer: 1.18 m?s7! 


(4) A rectangular orifice in the side of a large water tank has a breadth of 1 m. 
and depth of 0.5 m. The water on one side of the orifice is at a depth of 1 m 
above the top edge; the water level on the other side of the orifice is 25 cm. 
below the top edge. Determine the rate of discharge if the discharge coefficient 
may be taken as 0.62. 


Answer: 1.49 m3s-! 


(5) A spherical storage vessel 3 m in diameter contains a process liquid at 
half capacity. Determine the time taken to drain the vessel to a depth of 1 m 
through an orifice with a diameter of 25 mm at the bottom of the vessel. The 
discharge coefficient of the orifice is 0.6. 


Answer: 2341 s 
(6) Derive an expression for the time to drain an open hemispherical vessel 


of radius R through an orifice of area a, if the vessel has an inflow of liquid of 
constant rate Q. 


(7) Show that the time taken to drain a spherical vessel of radius R 


completely from full through an orifice of area a located at the bottom of the 
vessel is 
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5 
147R? 


dc UEM 
15Cqa 2g. 


where Cy is the coefficient of discharge and g the acceleration due to gravity. 


(8) A cylindrical vessel of diameter 3 m vertically mounted on its axis 
contains a liquid at an initial depth of 4 m and is allowed to discharge to a 
similar nearby vessel of diameter 4 m by way of a pipe with an inside diameter 
of 100 mm. Determine the final depth in both vessels and the time to reach this 
condition. Ignore any losses due to friction. 


Answer: 1.44 m, 502 s 


(9) Arectangular tank 6 m long by 2 m wide is divided in two parts by a parti- 
tion so the area of one part is twice the area of the other. The partition contains a 
sluice of area 100 cm2. If the level of liquid in the smaller division is 2 m above 
that of the larger, determine the time to reduce the difference in level to 50 cm. 
Assume a discharge coefficient of 0.6. 


Answer: 142 s 


(10) Two identical open tanks of cross-sectional area 4 m? are connected by a 
straight length of pipe of length 10 m and internal diameter 5 cm. A viscous oil 
with a viscosity of 0.1 Nsm-? and density 900 kgm^ is initially at a depth of 2m 
in one of the tanks while the other tank is empty. If a valve in the connecting 
pipe is fully opened, determine the time for the difference in level to fall to 5 mm. 
Neglect pipe entrance and exit losses and assume the connecting pipe is initially 
full of oil. 





Answer: 5100 s 


(11) In an experiment to determine the coefficients of contraction, velocity 
and discharge for a circular orifice of 9 mm diameter, water was discharged 
through the orifice mounted vertically in the side of a tank. A constant head of 
water of 1.2 m was maintained above the centre line of the orifice. The issuing 
jet was found to strike a target plate a horizontal distance of 850 mm from the 
vena contracta a vertical distance of 155 mm below the centre line of the 
orifice. Determine the values of the coefficients for the orifice if the measured 
discharge from the orifice was 91 litres in 470 seconds. 


Answer: 0.636, 0.985, 0.627 
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(12) Show that the time taken for liquid of density p and viscosity uto drain 
with laminar flow from tank of diameter d, and depth H through à vertical pipe 
of inside diameter d and length L, attached to the underside of the tank is given 





(13) Show, where the area of the tank is not appreciably greater than the area 
of the orifice, that the theoretical velocity of a jet of liquid v, fora constant head 
H, flowing from an orifice in the side of a tank, can be given by 





where a is the area of the orifice and A is the cross-sectional area of the tank. 


(14) The bottom of a process vessel has a conical section and contains a liquid 
which is required to be drained through an opening with a diameter of 25 mm. If 
the liquid in the section is initially at a depth of 1.5 m above the opening corre- 
sponding to a diameter of the section of 1.73 m, determine the time to drain the 
section completely. State any assumptions used. 


Answer: 883 s 


(15) A process vessel with a uniform area of 1.2 m? receives a liquid at a 
steady rate of 0.04 m?s-! and is simultaneously discharged through an opening 
at the bottom of the vessel. When the depth of liquid was 0.5 m the instanta- 
neous rate of discharge was noted as being 0.03 m?s-!, Determine the time 
taken for the level to rise by 20 m. 


Answer: 70 s 


(16) Water flows at a steady rate of 360 mth! into a vertical-sided tank of 
area 10 m2. The water discharges continuously from a sluice of area 0.0564 m?. 
If the initial level in the tank is 2.5 m above the sluice, determine the final depth 
after 5 minutes assuming a discharge coefficient of 0.6 for the sluice. 


Answer: 0.6m 
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Open channels, 
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Introduction 

Unlike the single-phase flow of fluids through pipes, the flow of liquids 
through open channels or flumes is characterized by a free surface normally at, 
or near, atmospheric pressure. Rivers and artificial canals are examples of open 
channels. Open channel flow also applies, however, to the flow of liquids 
through pipes which are not run full, as in the case of sewerage pipes. Weirs are 











vertical obstructions which lie across open channels, whilst notches are open- 
ings, normally rectangular or triangular, cut in such weirs. In prz the liquid 





in the channel builds up behind the weir until flow over the weir occurs with an 
equilibrium head of liquid being a measure of flowrate. 

There are numerous applications of open channel flow, notches and weirs. 
Although primarily encountered in the process effluent industry, there are many 
other applications in the chemical, power and oil industries. The perforated 
trays of distillation columns, for example, involve the open flow of liquid 
through which vapour is bubbled up from the beneath, with the depth of 
liquid maintained by a weir over which the liquid discharges into the 
downcomer. 





Many sophisticated mathematical methods and procedures have been devel- 
oped and applied to various types of open channel flow. Traditional methods 
are still widely used, however, based on the pioneering work of the French 
engineer Antoine Chézy in 1775 and the experimental work of the Irish engi- 
neer Robert Manning in 1891. In its simplest form, the channel is a 
have a gentle slope in which the rate of flow under the influence of gravity is 
balanced by frictional losses and such that the streamlines in the liquid run 
parallel to the floor (or bed) of the channel. Under conditions of uniform flow, 
there is therefore no change in velocity in the flowing liquid (no acceleration) 
and the depth is constant. Unlike single-phase pipe flow, the pressure condi- 
tions within the liquid are determined by hydrostatic principles and constant 
atmospheric conditions at the free surface. 





sumed to 
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The geometry of channels varies from wide and shallow to narrow and deep 
with side walls that may be vertical, inclined at an angle or curved. As open 
channels have a free surface, open channel flow therefore differs from pipe 
flow in that the cross-section of flow is not constrained. For increased rate of 
flow there is an increased depth. This affects not only the flow area but also 
wetted surface (both bed and wall) of the channel. An important parameter used 
to characterize the geometry of channels is the mean hydraulic depth defined as 
the ratio of flow area to wetted perimeter. 

For problems involving open channels, determining the rate of flow for a 
given depth of liquid is straightforward. Arriving at the depth for a given 
flowrate, however, usually requires trial and error approaches since the depth 
influences both the flow area and the wetted surface. The common method is to 
guess a value for the depth and calculate the corresponding flowrate, repeating 
the procedure a sufficient number of times to arrive at the answer. 

It should be noted that many of the equations associated with open channel 
flow are empirically based. Many of the available data have been obtained 
using water as the liquid medium. Empirical relationships used to predict the 
flow for other liquids should therefore be applied with caution. 
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7.1 Chézy formula for open channel flow 
Derive the Chézy formula for the uniform flow of a liquid along an open 
channel inclined with a slope i. 





F- mgsinü 


Solution 
Consider a liquid of density p flowing in an open channel of uniform cross- 
sectional area A, inclined at a small angle 0. Under steady uniform flow, the 
gravitational force perpendicular to the base of the channel is balanced by the 
viscous resistance forces acting in the opposite direction. The mass of the liquid 
in the channel can be expressed in terms of liquid density and volume 

F = mg sin 6 

= pALg sin 0 

The resistance forces which retard flow are obtained by considering the mean 
wall shear stress t, of the liquid in contact with the wetted channel wall. 

F =t„PL 
For equilibrium 

PALg sin 9 - t, PL 
Ifthe mean wall shear stress can be related to the kinetic energy per unit volume 


7: 
2 


p? 
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where f is a form of friction coefficient. Noting also that for small angles, sin 8 
is approximately equal to tan @ or the slope of the channel i, and defining the 
ratio of flow area A to wetted perimeter P as the mean hydraulic depth 


m=— 
P 


the average velocity of the liquid flowing in the channel can therefore be 
expressed as 


f 


That is, for a given friction factor and mean hydraulic depth the velocity is 
proportional to the square root of the slope i. Defining also 


c - [2 
f 
The flow through the channel may therefore be given by 


Q - CAN mi 


This is known as the Chézy formula where C is the Chézy coefficient. Experi- 
mental values and several correlations are available, the most common of 
which is the Manning formula given by 
1 
gm 

n 
where n is a dimensional roughness factor, the magnitude of which depends on 
the type of surface. Examples of typical roughness factors are given on page 
292, Note that as the roughness of the channel increases, the value of n also 
increases, reducing the value of the Chézy coefficient and therefore flow. 
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7.2 Flow in a rectangular open channel 


Water for cooling is delivered to a cooling tower along a rectangular concrete 
channel 1 m wide. Determine the rate of delivery for uniform flow at a depth of 
30 cm. The slope of the channel is 1:1000 and the roughness factor for the 
concrete is 0.014 m"! 3s, 


SS 


03m 
e im 


Solution — 
The delivery of water through the open channel is given by the Chézy formula 


Q =CAVmi (see Problem 7.1, page 169) 


where iis the slope, A is the flow area and C is the Chézy coefficient given by 
the Manning formula (see Problem 7.1, page 169). The discharge is thus 


1 
6 
Q =" Ami 
n 
aA ui 
n 
For the dimensions of the channel, the flow area is 
A=1x03 
-03m 
and wetted perimeter is 


P=2x03+1 
-16m 
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The mean hydraulic depth is therefore 





203 
16 
701875 m 
Thus, the rate of delivery is 


03 1 
= —— x 0.1875? x 40001 
g 0014 


20222 ms! 


The delivery of water through the channel is found to be 0.222 m?s-!. 
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7.3 Depth of flow in a rectangular channel 

A concrete-lined rectangular channel 12 m wide has a slope of I in 10,000. 
Determine the depth of water flowing in it if the volumetric flowrate is 60 m3s~!, 
Obtain the Chézy coefficient from the Manning formula taking the roughness 
factor as 0.015 m7", 


H 
SS 
I 12m 

Solution 
The flow through the channel is given by the Chézy formula 

Q=CAVmi (see Problem 7.1, page 169) 
where i is the slope and C is the Chézy coefficient given by the Manning 
formula 

1 
cam 


where n is the roughness factor. The flow area of the channel in terms of depth is 





A=12H 
and the wetted perimeter 
| P=12+2H 
The mean hydraulic depth (MHD) is therefore 
A 
m= 
P 
__12H 
1242H 
_ 6H 
6+H 
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An analytical solution for depth H is not possible. It is therefore necessary to 
use a graphical or trial and error approach illustrated below. 


Head Area Perimeter MHD Flowrate 
H.m A, m? P, m m,m Q, ms! 
40 48 20 24 57.36 
4.1 49.2 20. 244 59.38 
42 50.4 204 247 61.39 
4.13 49.56 20.26 245 59.98 


The depth is found to be approximately 4.13 m. 


62 


[2] 


59 


58 


FLOWRATE Q, m` s" 


57 


40 4| 4143 42 
HEAD H, m 
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7.4 Economical depth of flow in rectangular channels 

A rectangular channel is to be designed for conveying 300 m? of water per 
minute. Determine the minimum cross-sectional area of the channel if the slope 
is 1 in 1600 and it can be assumed that 


v 270 mi 


where v is the velocity of water in the channel, m is the mean hydraulic depth 
and i is the inclination of the channel. 


Y 
H 
S 
B 
pur] 
Solution 


For a constant flow through the channel the maximum flow occurs when the 
wetted perimeter is a minimum. The wetted perimeter is 


P=B+2H 
and flow area 
A=HB 


In terms of depth, the wetted perimeter is therefore 


paAson 
H 


Differentiating the wetted perimeter with respect to depth and equating to zero 





=0 
Therefore 
A-2H? 
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Thus 

B=2H 
That is, for the minimum cross-sectional area, the breadth of the channel is 
twice the channel depth. The mean hydraulic depth is therefore 


m 


M 
vi» 


v 
m 
n 


j 


> 
z 


w/z 


From the equation provided for velocity, the flow through the channel is 
Q -70A dmi 


In terms of depth H, the flow can therefore be expressed as 


uon? fi 
pe 2 


Rearranging in terms of channel depth 


3 
H =| 22 
140 \i 
300 P 
=j 160 x 42 x 1600 
140 
-16m 


The flow area is therefore 
A=2H? 
=2x 167 
=512 m? 


The minimum flow area is found to be 5.12 m?. 
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OPE 


7.5 Circular channel flow 


Determine the depth of flow of water in a sewer pipe of diameter 0.9144 m and 
inclination 1:200 when the discharge is 940 m3h~!. The Chézy coefficient may 
be taken as 100 m!?s~!, 


Vus 


Let r be the radius, 0 the half angle subtended at the centre by the water level 
and d the depth. From geometry 


coi ge eH 

R 
from which the angle 0 may be obtained in radians. The area of wetted 
cross-section is given by 


2 
A = POR ~ R° sin 80s 8 


- e(o E s) 
2 


and wetted perimeter is 


P=2R0 





for which the mean hydraulic depth is 
A 


m--— 
P 
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The flow through the channel is given by the Chézy formula 
Q - CAY mi (see Problem 7.1, page 169) 


where C is the Chézy coefficient and i is the channel slope. A graphical or trial 
and error approach is used for guessed values of depth: 


H.m cos 0 @ rad Am? — P.m m,m Q. n?h-! 
0.1 0.781 — 0.674 0.039 0.616 0.063 248.7 
0.15 0.671 — 0.834 0.070 0.763 04092 543.2 
0.2 0.562 0973 0.106 0.890 0.119 932.2 
0.25 0.453 1101 0146 1.006 0.145 14129 


0.201 0.560 0.975 0.107 0.891 0.120 940.0 


The depth in the channel is found to be 0.201 m. 


1500 
Ta 1000 
940 
S 
5 
z 
z 500 
0.201 
o 
o 01 02 03 
DEPTH H, m 
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7.6 Maximum flow in circular channels 
Determine the theoretical depth of | liquid in a closed circular channel of radius 
R for maximum velocity. 


Solution 

The velocity of liquid in a circular channel depends on the depth. As the 
velocity is proportional to the mean hydraulic depth, its maximum value may be 
obtained by 


pee a HP 9 
dO ' de 


where from geometry, the wetted area of flow is 


A=R20-R2 D 


(see Problem 7.5, page 177) 
i a | 





2 
and the wetted perimeter is 
P=2R0 


Differentiating A/P and equating to zero 
2R? (1 — cos 20) ea» - E 


Therefore 
20 = tan20 
The solution to which is 
20 =2575° 
The depth for maximum velocity is therefore 


2575° 
H =R - Rcos 2575 





-R-*062) 
=162R 


The depth for maximum velocity is 0.81 times the channel diameter. 
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7.7 Weirs and rectangular notches 

The weir above a downcomer in the rectifying section of a distillation column 
has a breadth of 1.2 m over which liquid flows with a head of 18 mm. Determine 
the rate of flow if the coefficient of discharge for the weir is 0.56. 








Flow 
— 





Solution 

To determine the flowrate of liquid over weirs (and notches) consider the flow 
of liquid along a horizontal plane at depth z below the free surface. It is assumed 
that the upstream velocity, sometimes known as the velocity of approach, is 
negligible in comparison to the flow over the weir. At a point just beyond the 
weir, the velocity of the liquid will have increased to some value, v, in which the 
static head is converted to velocity head 


vi 

ge 

2g 

Therefore 
v z42gz 


Thus, the velocity is proportional to the square root of the head. To determine 
the rate of flow, consider an elemental depth dz of the flow over the weir 
between depths z and z *dz. Since the elemental flow area is Bdz, the theoretical 
elemental flowrate dQ is therefore 


dQ = 2gz Bdz 
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The total flowrate of liquidQ is found by integrating over the depth of the weir 
or rectangular notch H. That is 


H 
[ao - 2s [zia 
0 0 


to give 
z 
9-25 PEL 


Tn practice, there are permanent energy losses due to edge effects giving rise to 
a flowrate which is lower than the theoretical. Introducing a discharge coeffi- 
cient C, , the actual flowrate is 


2 3 
Q - 5 Cub e H? 
In this case 
3 
Q - 5x 056x 12x J2 x 0018? 


7479 x 10? ms! 


This corresponds to a rate of flow of liquid over the weir of 17.25 m?h-', 

Note that the discharge coefficient takes into account the contraction of the 
overflowing jet as well as the effects of viscosity and surface tension (see 
Problem 4.2, page 102). Little is known of the separate influence of these two 
fluid parameters except that they become appreciable when the head on the 
weir and size of the weir decrease (see Problem 7.9, page 184), 
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7.8 Depth of a rectangular weir 

A hydroelectric power station is located on a dammed reservoir and is largely 
used to generate electricity at times of peak demand. In addition to the water 
which passes through the turbines, a flow of water continuously bypasses the 
turbines and is discharged from the reservoir over a rectangular weir into the 
river below. When the turbines are operating at peak demand, the bypass flow 
is ata minimum of 34 m3h~! and the fall from the reservoir to the river is 5.5m. 
When the power station operates at low demand, the flow of water bypassing 
the turbines is at a maximum of 217 m?hr^! and the fall is 8.5 m. Determine the 
height of the weir crest above the surface of the river. 





River bed 


Solution 
The flow over the weir of height H , into the river below for both cases is 


à 
2 


Qi - cub e s -Ho) (see Problem 7.7, page 180) 


and 
2 3 
Q) = Fr (Hj -HgY 


where B is the breadth of the weir. Assuming a constant discharge coefficient 
C, then Q, is related to Q5 by 


3 
9, (m Ho 
Q (Ho -Ho 
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Rearranging in terms of Hy 


2 
Q Hg 
H3 a -Hi 
Ho= 2 


2 
2; _, 
Q> 
2 
s5x(27) -85 
- 34 
2 
ar) 
34 
=427 m 


The weir is found to be a height of 4.27 m above the river bed. 
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7.9 Instantaneous flow through a rectangular weir 

Water used as a cooling medium in a process plant is stored in a large open 
tank which has a rectangular overflow weir. When the tank is overfilled the 
water cascades over the weir. Determine the time taken to lower the level of 
water to a head of 1 cm above the weir crest if it is noted that it takes 600 
seconds for the level to fall from an initial head of 8 cm to 6 cm over the weir 
crest when the inflow to the tank has ceased. 


Solution 
The elemental rate of flow through the weir is 


à 
do = E By2gH 2dr (see Problem 7.7. page 180) 


and the discharge from the tank is 
dQ =-AdH 


Therefore 
z 
ŽcaBJagH ?dr - -AdH 


The discharge over the weir is therefore given in general terms by 


kH? =—dH 


dr 
For an average head of 0.07 m then 


3 
kxoo7? =| 002 
600 


Thus 


k 218x107 m-I2s-1 
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The total time r to drain over the weir from an initial head H , to a final head Hy 
is found by integrating with respect to H 





2 3 xi 
=——— x (001 ? -008 7) 
18x107 


=7182s 
The time is found to be approximately 2 hours, Note that it would take an infi- 
nite length of time to empty the tank, H =0. This equation therefore does not 


hold for very small values of H, because a layer of water adheres to the weir 
crest due to surface tension. 


TIME: o6 
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7.10 Flow through a triangular notch 

A weak solution of caustic soda is fed to an absorber to remove sulphur dioxide 
from a process gas. The distributor to the absorber consists of channels with 
90° V-notches through which the caustic soda solution discharges. Determine 
the rate of flow through each notch if the head above the root (bottom) is 3 em. 
A discharge coefficient of 0.62 is assumed for the notches. 





Solution 

Triangular notches are generally used to measure small rates of flow in rivers 
and channels and can be used to channel flow for distribution purposes as illus- 
trated in the question. For the general case for flow over a triangular notch, let 
the notch semi-angle be 0 for which the elemental area lies at some depth 
between z and z-dz below the free surface. The theoretical elemental rate of 
flow, dQ, through the elemental area, is therefore given as 


dQ = /2g22(H — z)tan 0dz 


Introducing a discharge coefficient C ;, the total actual flowrate of liquid over 
the notch is therefore obtained by integrating from z=0 at the free surface to 
<=H at the bottom (or root) of the notch. That is 


Q H i 3 
fao =C42 42g tan 0 | (Hz? -2° )az 
0 0 
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Integrating with respect to z gives 


3 $ 
Q - C4 Dg 2 tan dins -k 


H 


0 
8 i 
= i 42 tan 9H? 
If the semi-angle is 45° (90° notch), the equation for liquid flowrate reduces to 
8 i 
QoS Cd \2g H? 
8 5 
- is (052 2e x 003? 
2228x107 ns! 


This corresponds to a rate of flow through each notch of 0.822 m3h-!, 

Note that, as with weirs and rectangular notches (see Problem 7.7, page 180) 
and discharging orifices (see Problem 6.2, page 145), the discharge coefficient 
C, allows for energy losses and the contraction of the stream cross-section. The 
actual discharge through the triangular or V-notch is therefore found my multi- 
plying the theoretical flow (or discharge) by Ca: 

In this analysis the velocity of the liquid approaching the notch is assumed to 
be considerably less than the rate of flow through the notch. The kinetic energy 
of the approaching liquid is therefore neglected. That is, the velocity through 
the horizontal elemental area across the notch is assumed to be dependent only 
on the depth of flow below the free surface. If, however, the velocity of 
approach can not be considered to be negligible, as may be the case for a narrow 
channel with a single notch, the total head producing flow will be increased by 
the kinetic energy of the approaching liquid such that the total flow. through the 
notch is found from 


2 3 2 e 
id Ei 
2g 


-3c, gno Hall | - 
15 2g 


where v; is the average velocity of the approaching liquid in the channel. 
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7.11 Tank drainage through a V-notch 
A sharp-edged V-notch in the side of a rectangular tank measuring 4 m long by 
2 m broad gives a calibration: 


s 
Q=15H? 


where Q is measured in cubic metres per second and H is measured in metres. 
Determine the time to reduce the head in the tank from 15 cm to 5 cm if the 
water discharges freely through the notch and there is no inflow to the tank. 


Solution 
Consider the level of water to be some depth H above the bottom of the notch. 
A small quantity of flow dQ would reduce the depth by a level dH in time dr. 
Then 
5 

dQ -15H'?dr 
with a change in capacity of the tank of 

dQ - -AdH 


The time to lower the level in the tank from an initial depth H ; to a final depth 
Hh is therefore found from 


0 H, 
Integrating 
ZAgu-R 8 
tss 2H à 
ais p 
2x4x2 -i -} 
=——— x (005 ? -015 ? 
3x15 f ; 
2257s 


The time is found to be approximately 4!; minutes. 
Note that the constant 1.5 is equal to 


152 Ca ie (see Problem 7.10, page 186) 


corresponding to a discharge coefficient C, of 0.635. 
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7.12 Flow through a trapezoidal notch 
Deduce an expression for the discharge through a trapezoidal notch which has 
abase Banda head H, and the sides of which make an angle 010 the vertical. 





Solution 
The flow can be considered as the flow through a triangular notch (see Problem 


7.10, page #86) and a rectangular weir (see Problem 7.7. 
i -7, page 180 
respective discharges are Pag ) where the 


8 $ 
Qnorch = gg C4 V28 tan OH? 
2 3 
Qweir = 3C Bog H? 
Combining the discharges 
8 5 3 
Srapnach = Ca I in 8H. Lc, p agir! 
8 a 
== Cy 2g H [ H tang + 3B 
15 4 


This equation may also be obtained from first principles for an elemental flow 
through the trapezoidal section and integrating over the total depth of flow H. 
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Further problems 
(1) Determine the theoretical depth of liquid in a closed circular channel for 
maximum discharge. 


Answer: 0.95 d 


(2) Aconcrete-lined rectangular channel is to be built to transport water from 
a nearby river to a power station as a supply of make-up water for the cooling 
towers. Determine the minimum depth of the channel if the rate of flow is not 
expected to exceed 650 m?min-! in a channel 5 m wide with a fall of 1 in 5000. 
Assume the Manning formula can be used to determine the Chézy coefficient, 
where the roughness factor for the concrete may be taken as 0.014 m-!3s, 


Answer: 2m 


(3) Explain why there should be a depth, less than full capacity, at which a 
pipe will carry water at a maximum rate of flow and suggest a way in which this 
can be determined analytically. 


(4) Determine the angle of the walls and the relationship between breadth 
and depth of a channel of trapezoidal cross-section necessary to deliver 
maximum flow. 


Answer: 30°, B =2/V3H 


(5) Water from a reservoir discharges over a rectangular weir of breadth 10 
m into a sluice below. When the reservoir is at maximum capacity, the rate of 
discharge over the weir is 5000 m?h-7! and the distance from the surface of the 
reservoir down to the sluice is 2 m. Determine the height of the weir crest and 
the discharge coefficient if the drop to the sluice is 1.9 m for a rate of discharge 
over the weir of 1000 m3h-!, 


Answer: 1.85 m, 0.81 
(6) Experimental data are recorded for the discharge through a 90° triangular 


notch for different heads. Plot these data on log-log paper and determine the 
coefficient of discharge for the notch. 


Flow,m?h-! 0.30 2.99 6.95 970 1302 16.94 
Head, m 0.00 0.05 0.07 0.008 0.009 0.10 


Answer: 0.63 
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(7) Derive an expression for the discharge through a weir where the area of 
the channel through which the water approaches the weir is such that the 
velocity of approach is not insignificant. 


(8) Derive an expression for the rate of flow of liquid over a rectangular weir 
of breadth B in terms of head H. 


(9) Determine the discharge coefficient for a rectangular weir placed in a 
channel 1 m wide of inclination 1:2500 with a Chézy coefficient of 66 m!/2s7! 
ifthe depth in the channel is 50 cm and the height above the weir crest is 32 cm. 


Answer: 0.62 


(10) Determine the rate of flow through a 60° V-notch if the head of liquid 
above the root is 5 cm. A discharge coefficient of 0.6 may be assumed. 


Answer: 1.65 m3h~! 


(11) Deduce an expression for the flow of liquid through a triangular notch of 
angle 0. 


(12) Show that the head of liquid over a weir on a sieve tray in a distillation 
column is given by 
n zons(Q 
B 
where Q is the rate of flow (m3s~!) and B is the breadth of the weir (m). 


(13) Determine the coefficient of discharge for the weir in Further Problem (12). 
Answer: 0.56 


(14) An aqueous solution of sodium hydroxide is fed continuously to a packed 
column to remove hydrogen chloride from a process gas. The liquid distributor 
above the packing consists of channels and the aqueous solution overflows 
through 60 identical 60° V-notches. Determine the head above the root of the 
notches in the distributor if the liquid flow to the column is 58 m3h-!. A 
discharge coefficient of 0.62 for the notches may be assumed. 





Answer: 4 cm 
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Pipe friction and 
turbulent flow 


Introduction 


Most applications concerned with fluids in pipelines involve turbulent flow. 
Great efforts have been made by scientists and engineers over the years to 
develop empirical relationships that predict the nature of turbulent fluid flow in 
both rough and smooth-walled pipes. But in spite of this, no exact solutions are 
available which will predict precisely the nature of the turbulent flow. 

French engineer Henry Philibert Gaspard Darcy in 1845 and German engi- 
neer and scientist Julius Weisbach in 1854, after much experimental work, first 
proposed.a friction factor equation that expressed the pressure loss in a piping 
system in terms of velocity head. In 1911, Paul Richard Heinich Blasius, a 
student of Ludwig Prandtl (1875-1953), demonstrated that for smooth-walled 
pipes the friction factor was dependent only on the Reynolds number. Blasius 
produced the first plot of its kind with friction factor versus Reynolds number 
for the empirical relationship 


uf 
f =0079Re * 





valid for Reynolds numbers between 4x10? and 1x105. Three years later, 
British engineers Sir Thomas Ernest Stanton (a former assistant of Osborne 
Reynolds) and J.R. Pannell established the friction factor/Reynolds number 
relationship to be independent of the fluid. Over the following years several 
other relationships were also established of the form 





f =a +bRe" 


a b c Validity 
Lees(1924) — LSxI03 — 0.152 — -035 4x 108 <Re<dx105 
Hermann (1930) 1.35x10 0.099 — -03 4x103<Re<2x106 
Nikuradse (1932) 8x10* 0.055 — 0.237. 4x10?<Re<3.2x106 
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To predict the friction of fluids in smooth pipes with turbulent flow, Prandtl 
developed a number of empirical models based on boundary layer theory. 
mixing length and wall effects. These models led to Prandtl’s universal resis- 
tance equation for turbulent flow in smooth pipes given by 


1 


=4log jo(Revf) -04 





and valid for Reynolds numbers between 5x10? and 3.4x10°. 

In 1858, Darcy conducted numerous detailed experiments involving pipes 
made from various materials. In 1932, German engineer Johannes Nikuradse 
analysed Darcy's data and noted that for turbulent flow in rough-walled pipes. 
the friction factor varied only slightly with Reynolds number, showing a 
decrease in friction factor with increasing Reynolds number. Above a certain 
limit, the friction factor was found to be independent of Reynolds number. 
Nikuradse also conducted his own work using artificially roughened pipes. 
Gluing carefully graded sand on to the inside of smooth pipes, he established a 
scale of relative roughness in which the diameter of sand particles in the pipe, €, 
is taken as the ratio of the uncoated pipe diameter, d. Using the data of Niku- 
radse in 1921, the Hungarian-born aerodynamicist Theodore von Kármán 
(1881-1963), a co-worker of Prandtl, developed an empirical relationship for 
turbulent flow through rough-walled pipes in the modified form 


=2ogu( £ )+ 228 





In 1939 C.F. Colebrook, working in collaboration with C.M. White, devel- 
oped a mathematical function which agreed closely with values obtained for 
naturally rough commercial pipes. This involved using both Prandtl's smooth 
pipe law of friction and combining it with von Kármán's fully rough pipe law of 
friction into a single empirical expression in the form 


-L 228 -4log TE A515 ) 


F d Redf 


It was Louis F. Moody who. in 1944, first presented a composite plot of fric- 
tion factor with Reynolds number which included the straight line relationship 
for laminar friction factor, smooth pipe turbulent friction factor and rough pipe 
turbulent friction factor, together with the concept of relative roughness. 
Known as the Moody plot (see page 293), it continues to be a valuable tool for 
evaluating friction factors for pipe flow in which it is evident that in the 
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turbulent region, the slope of the curve relating friction factor to Reynolds 
number decreases with increasing Reynolds number, becoming independent of 
Reynolds number at high values. 

Nomenclature and scientific reasoning are both major difficulties encoun- 
tered in the field of pipe friction. There are several widely-accepted approaches 
used by industry, each having a validity which — irrespective of procedure — 
provides identical final answers. To the unwary, the approaches have the 
propensity to confuse. As a comparison, the friction factor developed in 1893 
by the American engineer John Thomas Fanning (1837-1911), and known as 
the Fanning friction factor f, is related to other commonly-encountered forms 
and symbols for the friction factor where 

S =C =fr skato aR 
4 4 p? p? 


A common problem in calculating pressure loss through pipes in which the 
flow of fluid is known to be turbulent but of unknown relative roughness and/or 
flowrate, and thus unknown Reynolds number, is that the Fanning friction 
factor cannot be readily obtained. In such cases it is usual to assume either a 
value of 0.005 to simplify the problem, or to use this as a basis to begin an itera- 
tive procedure converging on a solution. 
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8.1 Economic pipe diameter 


Explain what is meant by economic pipe diameter. 


Minimum 
cost Capital 
charges 
5 
e 
o 
Operating. 
Costs. 
EPD 
DIAMETER 
Solution 


The economic pipe diameter is the diameter of pipe which gives the minimum 
overall cost for any specific flowrate. In general, however, pipes are determined 
by selecting a reasonable fluid velocity which provides a reasonable pressure 
drop and is virtually independent of diameter. It is possible, however, to deter- 
mine the economic pipe diameter based on economic principles by which 
capital charges for pipes and valves are balanced against pump investment and 
operating costs. An economic optimum can be readily determined where the 
cost of pipe material appropriate to handle a particular fluid is correlated with 
pipe diameter. Once a calculation has been made, it is necessary to determine 
the Reynolds number to ensure that flow is turbulent and that other essential 
requirements are met, such as net positive suction head in the suction lines to 
centrifugal pumps. 

The minimum investment is calculated for expensive or exotic pipe mate- 
rials such as alloys, pipelines larger than 300 mm in diameter and carbon steel 
lines with a large number of valves and fittings. The pipe scheduling is selected 
by determining either the inner or outer diameter and the pipe wall thickness. 
The minimum wall thickness is a function of allowable stress of the pipe mate- 
rial, diameter, design pressure, and corrosion and erosion rates. 

In the case of highly viscous liquids, pipelines are rarely sized on economic 
considerations. 


196 





PIPE FRICTION AND TURBULENT FLOW 


8.2 Head loss due to friction 


Derive an equation for the pressure and head loss due to friction for a fluid 
flowing through a pipe of | length L and inside diameter d. 








Solution 
Let the pressure difference or drop be solely due to friction manifest as a wall 


shear stress, T, . For steady state conditions, a force balance on the fluid in the 
cross-section of the pipe is 


2 2 
20 ond nd? 
Py Pe 47 Ap; a 
= ty nld 
The pressure drop due to friction is therefore 
L 
Apy -4t, d 


If the wall shear stress is related to the kinetic energy per volume, then 


f 
2 





t, 





p? 


where f is the Fanning friction factor. The pressure drop due to friction may 
therefore be expressed as 
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or in head form 


AfL v? 
H; == 
f d 2g 


This is known as the Fanning or Darcy equation. It was much earlier, however, 
that Darcy in 1845 and Weisbach in 1854 first proposed the friction factor equa- 
tion after much experimental work, giving rise to the head loss due to friction in 
the form 


v? 


Hy 
J 2g 


RE 


This is known as the Darcy-Weisbach equation where the friction factor is 
related to the Fanning friction factor by 


A-Af 
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8.3 General frictional pressure drop equation 

applied to laminar flow 

Deduce a relationship between the friction factor used in the Darcy equation 
and Reynolds number for a fluid with fully developed laminar flow. 


Solution 
The average velocity of a fluid with laminar flow through the pipe is given by 
ys Lg (see Problem 3.7, page 73) 
8u L 


Rearranging, the pressure drop along the pipe is therefore 








The Darcy equation for frictional pressure drop along a pipe is 


fov? L 
AEn 


Combining the two equations in terms of pressure drop 


2pvL 3ouLv 





d d? 
which, in terms of the Fanning friction factor, reduces to 
_ 164 
sie prd 
af 
Re 


That is, for laminar flow the friction factor is inversely proportional to 
Reynolds number and when the log of Reynolds number is plotted against the 
log of the friction factor (Moody plot) gives a straight line with a gradient of 
—1. Note that, using the Darcy-Weisbach equation, this is equivalent to 
jat 
Re 
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8.4 Blasius' equation for smooth-walled pipes 
Ethylbenzene, with a density of 867 kgm and viscosity of 7.5x10 Nsnr?, is 
10 be transferred at a rate of 12 mh! through a smooth-walled pipeline 200 m 
long under the action of gravity. Determine the internal diameter of the pipeline 
if there is a fall in elevation along the pipeline of 10 m. Assume the Blasius 
equation can be applied to determine the friction factor. 


Solution 
The head loss due to friction is equal to the static head. That is 


H=H; 


= Af y? 
d 2g 


where the friction factor is given by the Blasius equation 


ty 
f =0079Re 4 








and where the Reynolds number expressed in terms of volumetric flowrate is 
Re = 02 
ndu 
Therefore 
e 
4 
4x 0.079 ipo L 42 
H= Tdp nd? 
2dg 
zi 
=i 
4x 0079 x| #2 2 a 
1 nyu x 
2g 
k 
4x 867 x i j“ 4x D. Š 
4x 0079 x| 3600 | x 200x| —3600 | x q7: 
xx75x107 E 
107 - 
2g 
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This reduces to 


19 
144x10° =d * 


log, 144 x 10° =P bee, d 


log, d - -2985 


Solving, the diameter of the pipeline is found to be 0.0505 m (50.5 mm). A 
check for Reynolds number gives 


Re = PQ 


zdu 


4x867x 2 
3600 


/ RX 00505 x 75 x 107 
-97153 


This value lies between 4x10? and 1x105, confirming the validity of using the 
Blasius equation for smooth-walled pipes. 
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8.5 Prandtl’s universal resistance equation for 
smooth-walled pipes 

A benzene mixture, with density of 873 kgm and viscosity of 8.8«10~4 Nsm~, 
is fed from an open overhead tank to a process through a smooth pipe of 25 mm 
bore and length 18 m. If the process operates at a pressure of 55 kNm- above 
atmospheric, determine the minimum allowable head of the liquid surface in 
the tank above the process feed point if the flow is not to fall below 40 kemin-. 
Allow for entrance and exit losses and apply Prandtl's universal resistance 
equation for turbulent flow in smooth pipes. 





Open tank 


Solution 

Assuming that the velocity in the overhead tank is small in comparison to the 
velocity in the pipe (v;=0), then applying the Bernoulli equation between 
points | and 2 





© 873 x nx 0025? 


=155 ms! 
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For Reynolds number 
"E 
u 
_ 873 x 155 x 0.025 
(0 88x10 
= 38441 


which satisfies the criteria for Prandtl’s equation (5x10°<Re<3.6x10°), The 
Prandtl equation is 
1 





p 748e Vf) -04 


= 4log 938441 Vf ) -04 


By trial and error 


f LHS RHS 
0.004 15.81 13.14 
0.005 14.14 13.33 
0.006 * 4291 13.49 
0.0055 13.48 13.42 


That is, a friction factor of approximately 0.0055. The total head loss is due to 
pipe friction, entrance and exit losses 


2 2 2 

v v v 
Hg aa 952 41922. 
d 2g 2g og 





Rearranging the Bernoulli equation, the minimum allowable head is therefore 
2 


2 
= vy Y 
H=P2 P1 4 "2 4 "2 (Mh is 
pg 2g 2g\d 








3 2 2 
258x105 , 155? | 155 | (4x 00055 x18 | i. 
873xg 28 — 2g 0025 


-866m 


The minimum allowable height of liquid in the tank above the process feed 
point is 8.66 m. 
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8.6 Pressure drop through a rough-walled 

horizontal pipe 

Two ethanol storage tanks are connected by 100 m of straight pipe. Both tanks. 
are open to atmosphere and the connecting pipe has an inside diameter of 50 
mm and relative roughness of 0.002. Determine the pressure drop over the 
length of pipe if the flowrate is 15 m?h-! and estimate the difference in level of 
ethanol between the two tanks. There are no extra fittings in the pipe but the 
entrance to the pipe and the exit from the pipe should be taken into account. The 
density of ethanol is 780 kgm and its viscosity is 1.7*10-3 Nsm~. 





Solution 
The average velocity of ethanol through the pipe is 





which corresponds to a Reynolds number of 
rud 
u 
_ 780 x 2.12 x 005 
17x10? 


= 48,635 
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The relative roughness for the pipe, €/d, is 0.002. From the Moody plot (see 
page 293) the friction factor is 0.0065. The head loss due to friction along the 
pipe is therefore found from 


2fpv?L 
MU 
_ 2x 0.0065 x 780 x 2.12? x 100 


005 
-91146 Nm? 


The difference in levels between the two tanks can be determined by applying 
the Bernoulli equation at the free surface of both tanks from which the head loss 
is therefore equal to the static head. That is 


H-H, 


= Hp +H exi +H entrance 





212? _ (4x 00065 x 100 
=x Ons +15 


=1225m 


The pressure drop due to friction is 91.1 kNm-? and the difference in levels is 
1225 m. 
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8.7 Discharge through a siphon 

A liquid of density 1150 kgm is siphoned from an open vat using a tube of 
internal diameter 25 mm and length 40 m. The tube rises vertically from its 
upper end a distance 10 m to the highest point. The discharge end is 2 m below 
the upper end and is open to atmosphere. Determine both the minimum height 
of liquid allowable in the vat above the open end, and the rate of flow. The 
minimum allowable pressure of the liquid is 18 kNm-?. Assume a Fanning fric- 
tion factor of 0.005. 





Solution 

Siphoning is a useful technique for decanting a liquid from a vessel where there 
may be a layer of sediment which must not be disturbed. It involves the transfer 
of liquid to another vessel at a lower elevation by means of a pipe or tube whose 
highest point is above the surface of the liquid in the upper vessel. The device is 
started by filling the siphon tube with liquid by applying pressure on the upper 
surface or suction at the outlet. Liquid then continues to rise to the top point and 
discharge continuously under the influence of gravity to the lower vessel. A 
limitation, however, is that the pressure at the highest point must not fall below 
the vapour pressure of the liquid. Should this occur, vapour will be released 
from solution (boiling) and cause a break in the liquid stream. The siphoning 
action will also cease if the pipe entrance is no longer below the upper liquid 
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surface. Applying the Bernoulli equation between the free surface (1) and the 
top point in the tube (2), then 





PL.asPR..1.H, 


pe pg 2g 


where the head loss due to friction and entrance loss is 





vi 2 
Hy chr. p05 
d 2g 2g 
2 
24x00 x10, v? o. v? 
0025 2g 2g 
=g 
2g 


The Bernoulli equation is therefore 


3 E 
JUNSPUIU" vm 18x10 "2 siosgs t 
1150x g 1150x g 2g 2g 


which reduces to 
2 
vo 
H =262+95— 
2g 
Applying the Bernoulli equation between the liquid surface and discharge point 
which are both at atmospheric pressure 


"5 
H+2=2+H, 
28 


where the losses are due to tube friction, entrance and exit losses 


2 
LALY Qs. ag 
d 2g 2g 2g 


4 2 
-.4x0005x40. v 15x“ 
0025 2g 2g 


52 
23351— 
28 
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Therefore 


2 2 jd 
262495!— 42-21. 4335! — 
2g 2g 2g 
Solving, the velocity through the tube is found to be 1.9 ms~!. The corre- 
sponding flowrate is therefore 


2 
o-n 


mx 0.0257 
-———x 
4 


19 


7933x107! ms! 


and the minimum depth in the vat is 
2 


V 
H 22624952 
2g 


2 
=262 +95x p 
2g 

=4.37 m 


The minimum allowable height of liquid in the vat is 4.37 m. The flowrate in the 
siphon tube is found to be 9.33x10-4 m?s-!. Note that this is a problem 
involving variable head flow. The flowrate is therefore not constant in the tube 
and is at a maximum when the level in the vat is at its highest position above the 
open end. 
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8.8 Flow through parallel pipes 

A liquid flows through a short pipe which branches into two parallel pipes A 
and B each with a length of 50 m and with inside diameters of 25 mm and 50. 
mm, respectively. The ends of the pipes are connected. together by another short 
pipe. Determine the flow through each pipe if they have a drop in elevation of 
3 m. Assume a constant Fanning friction factor in both pipes of 0.005. 


L=50m i 
Pipe A —7Q,  d,-s0mm 


Lu ere ae 


Flow —> 


Pipe B —90, — d,-25mm 


Solution - 
The head loss due to friction in pipe A is 


EA 
fd 3g 


24X0005x 50. VÀ 
0025 2g 
2 
= 4024 
2g 


and for pipe B is 


2 
xdi ts. 
d 2g 


LAx0005x 50. v? 


005 2g 
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As the pressure is the same at the end of each pipe then 


2 2 
v v 
40-4 =20 8 
2g 2g 
=H 

=3m 


from which the velocities v 4 and v are found to be 1.21 ms^! and 1.71 ms! 
respectively. The rate of flow through the pipes is therefore 


2594 x 107* ms! 


and 
nd? 
B 
ia 
Qn 3-778 
_ nx 005? 


x171 


2336x107 ms! 


The rates of flow through parallel pipes A and B are therefore 2.14 m?h-! and 
12.08 m?h-!, respectively, corresponding to a total flow of 14.22 m?h-'. 


210 





PIPE FRICTION AND TURBULENT FLOW 


8.9 Pipes in series: flow by velocity head method 

Two water reservoirs are connected by a straight pipe 1 km long. For the first. 
half of its length the pipe is 12 cm in diameter after which it is suddenly reduced 
to 6 cm. Determine the flow through the pipe if the surface of the water in the 
upper reservoir is 30 m above that in the lower. Assume a friction factor of 
0.005 for both pipes. 





Solution ` 
From continuity, the velocity in both sections of pipe of circular cross-section is 
related to flow for the incompressible fluid as 





Basing the calculations on the smaller diameter pipe, the head loss at the 
entrance to the larger pipe is 
2 2 


v 
05—- 2003125 m (see Problem 2.6, page 47) 
g 


eL. 
28 
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The head loss in the 12 cm inside diameter pipe due to friction is given by 


2 
Mb 
d 2g 


or in terms of the smaller 6 cm inside diameter pipe 


2 
LL Afta v5 


16d 2g 


_4x0005x500 73 
16x012 2g 
2 
v 
252-23 
28 
The entrance head loss at the sudden contraction 7. between pipes is found 
using the chart provided in Problem 2.6 on page 47, and uses 
2 
V 
H, =k d 
2g 


in which k is found by interpolation to be 0.44, noting that the ratio of pipe areas 
ag /ay is equal to 0.25. That is 
" 
H, =044— 
28 


The head loss due to friction in the 6 cm inside diameter pipe is 


2 
= {ln "2 


per 
_ 4% 0005 500 "3. 
006 2g 

2 


y. 
=1667-> 
2g 
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The head loss at the exit is 
v? 
dis D (see Problem 2.5, page 44) 


The total head loss H between the two reservoirs is therefore 
x 
H = F (003125 + 52 + 044 +1667 +10) 
& 


2 

v 
=17342 
2g 


This is equal to the static head between the reservoirs. That is 


2 

v 
30217342. 
2g 


Rearranging, the velocity in the smaller pipe is therefore 


3 





n= 


1734 


=184 ms! 


The rate of flow between reservoirs is thus 
2 
2 
m—y. 
[^ 4 
xx 006? 
4 


X184 


52x10? m^! 


M 


The rate of flow is found to be 5.2x10-3 m3s-! or 1.872x10-? m?h-!. 
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8.10 Pipes in series: pressure drop by equivalent 

length method 

Kerosene, with density of 815 kgm™ and viscosity of 7x10-3 Nsm~, flows 
through a pipe system at a rate of 12x10? kgh"!. The pipe system consists of 
50 m of 50 mm bore straight pipe with two 90° elbows, followed by a reducer to 
a 38 mm bore section 50 m in length with two further 90? elbows. The pipe has 
an absolute wall roughness of 0.02 mm. Determine the pressure loss due to fric- 
tion in the system if the equivalent length of each elbow is 40 pipe diameters and 
the reducer is equal to 0.2 velocity heads. Assume that the von Kármán equa- 
tion for rough-walled pipes applies. 


Reducer 
38 mm. 





Flow —> 


Solution 
The velocity in the 50 mm bore section is 


3 
4x 2x10 
" 3600. 
815x xx 005? 


-208 ms“! 


214 





PIPE FRICTION AND TURBULENT FLOW 


Applying the von Kármán equation 


1 d; 
-= =2log (a +228 
vi 7 
50 
-2x]logjg| — |4 228 
ei aap ) 
gives a friction factor of 0.0121. The frictional head loss through the pipe and 


two 90° elbows expressed in equivalent length is given by 


Silea v? 


50 mm d 2g 


5 4x 00121 x (50 +2 x 50x 005) . 208? 
005 2g 


=1164 m 
The velocity in the 38 mm bore section is 


. 3 
4x 12x10" 
js 3600 
815 x mx 0.038" 


= 361 ms"! 


The friction factor is obtained from the von Kármán equation 


5 


-L = 210g 19( 22 )+228 
dn e 
38 
-2x lo —— {+228 
(as 


to give a value of 0.0128. The head loss due to friction in terms of equivalent 
length is therefore 


m _ 4X 00128 x (50 + 2 x 50 x 0038) 3.617 
ee Mies 
0038 2g 


=4815m 


215 















































FLUID MECHANICS 


The head loss at the reducer, H , , based on the smaller bore pipe is 


=013 m 
The total head loss is therefore 
H, =1164 +4815 + 013 
=5992 m 
The pressure loss due to friction is therefore 
Aps =P8H, 
=815 x g x 5992 
= 479.07 x 10° Nm? 


The total pressure drop due to friction through the pipe system is 479 kNm. 
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8.11 Relationship between equivalent length and 
velocity head methods 


Determine the equivalent friction factor relating the head loss for flow of a fluid 
around a 90° elbow which is expressed as both 1.2 velocity heads and 60 pipe 
diameters. 


Solution 
The head loss for the elbow expressed in terms of velocity head is 


y? 
H, =12— 
2g 


and expressed in terms of equivalent length, the equation for frictional head loss 
for the elbow is 


Heeg v 
d 2g 





H, 


where the equivalent length of pipe is 
Leq 7604 


Note that the equivalent length of fitting is that length of pipe which gives the 
same pressure drop as the fitting. Since each size of pipe (or fitting) would 
require a different equivalent length for a particular fitting, it is usual to express 
equivalent length as so many pipe diameters which is therefore independent of 
pipe. Therefore 


2 
H, =240f — 
2g 

Combining the two equations for head loss 


i2 =240f 
2g 2g 


Therefore 
12 
f^ 
= 0.005 
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That is, the equivalent Fanning friction factor is found to be 0.005. This rela- 
tionship between velocity head and equivalent length of pipe holds well for 
high Reynolds numbers. At low Reynolds numbers in the turbulent region the 
two methods deviate due to the influence of pipe surface roughness. In 1943, 
American hydraulics engineer Hunter Rouse gave a limiting equation that 
distinguished between the transition region and the fully rough regime where 


100 


D) 


Beyond the Reynolds number predicted by this equation, the friction factor 
becomes essentially independent of the Reynolds number. The Rouse limit line 
can be included on the Moody plot (see page 293). 





Re = 
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8.12 Flow and pressure drop around a ring main 

Water is supplied to a small laboratory four-side ring main ABCD at A and B, 
with flowrates 3x10 m3s-! and 0.5x10-3 m?s-!, respectively, and is drawn at 
rate C and D with the flowrate at D being 0.9x10~3 m3s~!, Determine the rate 
of flow from B to C and the lowest pressure drop in the main. The connecting 
pipes each have lengths of 10 m and internal diameters of 25 mm. The Fanning 
friction factor is assumed to be 0.005. 


H 





Solution 


The pressure drop around the main is 


Apap * Apgc 7 Apap * Appc 


where 
2 
Apap = t 
— 2fL[ 43x10? -Q) 
d ad? 


_ 2x 0.005 x 1000 x 10x 16x 107 x (3 -Q* 
7x 0025? 


-16600(3 -Q» 
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Likewise 
Apap = 16,6000" 
Ap pc = 16,6000 +05) 
Appc = 16,600(21 -0) 
where Q is expressed in litres per second, Therefore 
Q? «(Q «05? -G-0) «Q1-0 
which reduces to 


112Q -1316 20 


Solving, the flowrateQ is 1.175 litres per second. The rate of flow from B to C is 
therefore 1.675 litres per second. The corresponding pressure drops through the 
pipes are therefore 


Apap 7555 kNm ? 
Apag 7228 kNm ? 
Appc 7464 kNm? 
Appc 7143 kNm? 


with the lowest pressure drop in the pipe connecting C and D. Note that where 
the friction factors, lengths and diameters are not equal, a quadratic equation in 
terms of flowrate arises. 
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8.13 Tank drainage through a pipe with turbulent flow 
Water flows from an open cylindrical tank of cross-sectional area 4 m? through 
a 20 m length of horizontal pipe of 25 mm inside diameter. Determine the time 
10 lower the water level in the tank from 2 m to 0.5 m above the open end of the 
pipe. Assume a friction factor of 0.005 and allow for entry and exit losses. 





Solution 

Atmospheric pressure is exerted both at the free surface in the tank and at the jet 
issuing from the open pipe. Applying the Bernoulli equation where it is reason- 
able to assume that the velocity of the water in the tank is far less than the 
velocity in the pipe, the static head is therefore 


2 
H= +H, 
2g 


where the head loss for pipe friction, exit and entrance loss is given by 
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Rearranging 
28H 


4f. 55 
d 


v= 


Applying an unsteady state mass balance over the tank, the change in capacity 
of the tank is equal to the rate of flow through the pipe. That is 


ee 
dt 
Substituting for pipe velocity and rearranging, the total time to lower the level 


from H to H3 is 


5 a Waasu T 
dis — 4 f Han 
p i 


Integrating 


afs | 
mu (HÊ -H 


; 
ie i Hh 
afi 


2x4x PX OO ED 425 - 
5 E x0 -08) 
xil SJE 


=7912s 


The time taken is found to be 2 hours and 12 minutes. 
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8.14 Turbulent flow in non-circular ducts 

Hot water, at a temperature of 80°C with a corresponding density of 972 kgm? 
and viscosity 3.5x10~ Nsm~, flows at a rate of 50 m3h! through an item of 
process plant which consists of a horizontal annulus consisting of two concen- 
tric tubes. The outer tube has an inner diameter of 150 mm and the inner tube 
has an outer diameter of 100 mm. Determine the pressure loss due to friction 
per unit length if the surface roughness of the tubing is 0.04 mm. Use the Moody 
plot to obtain the friction factor. 


i 
NSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSASSSSSNN 
Flow D, — = Py 


SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSASSSSAS 


SSSSSSSSSSSSSSSSSSSSSSSSSSAASASSSSNASNSS 
m ^ B aes n 
— +e w 4 
SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 


Solution 
For the flow through a non-circular pipe or duct, consider a horizontal duct of 
any shape with uniform cross-section. A force balance on the fluid in the 
cross-section gives 
(Pp, —p2)a = Appa 
= TPL 
where P is the wetted perimeter of the wall, Lis the length of element andais the 


uniform cross-sectional area. Relating the wall shear stress to friction factor in 
the form 


w 


m 
Ty cav 





_ fov? LP 
2a 


py 


Comparing this to the equation for pressure drop due to friction in a circular 
pipe 





2fp 
d 





App = 
it is noted that 4a/P has replaced d. Therefore defining the term deq . known as 
the equivalent hydraulic diameter, as 

eq P 


the frictional pressure drop equation for non-circular pipes can be expressed in 
the form 





d, 


eq 


The flow area for the tube annulus is 


a- ida -d?) 
- T5015 -01)) 
4 


=982 x107 m? 
and the wetted perimeter is 
P — nd, dj) 
= nx (015 01) 
=0785 m 
The equivalent hydraulic diameter is therefore 


_4x982x 10 
Sp OTS 


=005m 
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The relative roughness based on the equivalent hydraulic diameter is therefore 


The average velocity through the annulus is 


Q 


yest 
a 


50 
= 3600 
982 x 10° 


=141 ms! 


The Reynolds number is 
pe = tes 
p 
_ 972x141x005 
O 35x10% 
= 195,788 


From the Moody plot (page 293), the Fanning friction factor is 0.0052. The 
pressure loss per unit length due to friction is therefore 


App 2v 
L deg 
_ 2x 00052 x 972x141? 
005 
-402 Nm? m^! 


That is, the pressure drop due to friction is found to be 400 Nm? per metre 
length of tube. 
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8.15 Head loss through a tapered section 


Water flows through a conical section of pipe which narrows from an internal 
diameter of 30 cm to an internal diameter of 10 cm. The total length of the 
section is 3 m. Determine the frictional head loss for a flow of 0.07 m? if the 
Fanning friction factor is 0.005. Ignore entrance and exit losses and effects due 
to inertia, 


L (3m) 


e—a 


Not to scale 


Flow Ine ee 
3 D cerae 


Solution 
The head loss due to friction is given by 


avt 


H 
Jd 28 


or in terms of flowrate 


4Q 
AfL | nd? 


d 2g 


.3NL9* 
T. gd: » 
From geometry, the diameter of the section is related to length by 


d —0.3 - 00666 L 








where diameter d and length L are measured in metres. The total frictional loss 
is therefore obtained by integration over the length of the section. 
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That is 
H, 23 
fa, = 22° je 
S mg (03 00666) 


The integration is simplified using the substitution 
u —0.3 — 00666L 
for which the differential is 
du = -0.0666dL 
The integration required is therefore 
T "EE ur n 
00666 n° g 


Thus, integrating with respect to u gives 


c» ut 
Hy = 
no g 4 
That is 
s 
ses ~32f0? | (0.3 00666) ^ 
1 90666772 =4 y 


00666 x T? x g 4 
-03m 


.32x0005x007 . [s -00666 x 3)* -037* | 


The frictional head loss through the section is found to be 0.3 m. 
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8.16 Acceleration of a liquid in a pipe 

Two large storage tanks containing a process liquid of density 1100 kgm™ are 
connected by a 100 mm inside diameter pipe 50 m in length. The tanks operate 
with a constant difference in level of 4 m and flow is controlled by a valve 
located near the receiving tank. When fully open, the valve has a head loss of 13 
velocity heads. If the valve is suddenly opened, determine the time taken for the 
process liquid to reach 99% of the final steady value. Assume a friction factor 
of 0.006 and neglect entrance and exit losses. 








Solution 

The process liquid is assumed to be inelastic so that there is no pressure wave. 
The average steady velocity of the liquid under full flow may be determined 
from the case where, for the freely discharging liquid, the static head is 
balanced by the frictional resistance to flow due to the pipe wall and the valve as 





2 2 
v v 
Heigl 
d 2g 2g 
Rearranging 
y= 2gH 
1 
afk 413 
d 
m 2xgx4 
- 4x 0006 x 50 
— +n 
01 
2314 m^? 
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Therefore 


v, 2177 ms! 
For the accelerating liquid, the inertia head is also included where it is assumed 


that the liquid is inelastic. Thus 
2 3 
paypi Lue 
dag 2g gd 








Rearranging 
28H — 2, 2L ode 
3L as ME igi 
d 
sieh 2x50 dv 
i= ee 
450006 x 50 5 d 
01 
=v? +4 
dr 
Rearranging 
Adv. 
dt = ——; 
v? -y? 








oot) Ag 1 e 
vilvi tv vi 7v 


Integrating between the limits of no velocity (stationary liquid) and 99% of the 
final steady value for velocity, the time taken is 


+y 0 
E v 
OPES 

vy vpn» lo 
22 uis, [1775099 x 177 
Tum" m -ossxim 


-598s 





The time taken is found to be approximately 6 seconds for the liquid to have 
reached 99% of its final velocity. 
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Further problems 

(1) Determine the pressure drop due to friction along a | km pipeline of 
inside diameter 254 mm used for drawing sea water at a rate of 500 m?h-! into a 
chemical plant for use as a process cooling medium. The friction factor is 0.005 
and density of sea water 1021 kgm^?. 


Answer: 302 kNm? 


(2) Determine the pressure drop due to friction along a pipe of inside diam- 
eter 100 mm, 500 m in length, carrying a process liquid of density 900 kgm-^ at 
a rate of 46.8 m?h7!. The Fanning friction factor is assumed to be 0.005. 


Answer: 123.3 kNm? 


(3) Using the Darcy equation, determine the pressure drop along a 10 m 
length of pipe of inside diameter 25.4 mm carrying a process liquid of density 
1010 kgm^? and viscosity 0.008 Nsm-? at a rate of 0.015 m?min-', 


Answer: 1.91 kNm^? 


(4) Determine the pressure drop due to friction along a | km length of pipe of 
inside diameter 200 mm carrying a process liquid at a rate of 125 m?h-!. The 
process liquid has a density of 950 kgm™ and viscosity of 9x10~ Nsm? and 


the pipe wall has an absolute roughness of 0.04 mm. 


Answer: 48.8 kNm-? 








(5) Determine the head loss due to friction in terms of the number of velocity 
heads for a flowing fluid in a pipeline of length 2.5 km and inside diameter 100 
mm with a Fanning friction factor of 0.006. 


Answer: 600 v? /2g 


(6) An experimental test rig was used to determine the parameters a and b 
relating the friction factor to Reynolds number in the form 


y -aRe" 


The rig consists of 10 m of smooth-walled pipe with an inside diameter of 16 
mm and uses water as the liquid medium. Results from two trials gave frictional 
pressure drops of 13.9 kNm- and 20.0 kNm-? for flows of 0.94 m?h7! and 1.16 
mh"! , respectively. Determine the values fora and b. 


Answer: a = 0.079, b = -0.25 
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(7) A distillation column tower operating at a gauge pressure of 700 kNm-? 
receives a hydrocarbon feed of density 780 kgm™ and viscosity 0.6x10-4 
Nsm- at a rate of 486 m?h-! from a neighbouring column operating at a gauge 
pressure of 1.2 MNm~. The connecting pipework consists of 40 m of 25 cm 
inside diameter pipe with an absolute roughness of 0.046 mm and has a total of 
seven elbows and a single control valve. If the level of hydrocarbon mixture in 
the base of the feed column is steady at 4 m above the ground and the feed point 
of the receiving column is at an elevation of 15 m above the ground, determine 
the pressure drop across the control valve, Allow for entrance and exit loss 
and a loss for each bend equivalent to 0.7 velocity heads. 





Answer: 387 kNm? 


(8) Show that the head loss due to friction, Hy, for the flow of fluid in a 
circular pipe is given by 


where f iş the friction factor, L is the length of pipe, v is the average velocity, d. 
is the diameter of the pipe and g the gravitational acceleration. 





(9) Asmooth pipe of internal diameter 150 mm is used for transporting oil of 
density 854 kgm-? and viscosity 3.28x10-3 Nsm~2. The pipe has a length of 1.2 
km and an elevation of 15.4 m. If the absolute pressures at the lower and upper 
ends are 850 kNm-? and 335 kNm^2, respectively, determine the rate of flow 
assuming the Blasius equation applied for smooth pipes. 


Answer: 0.0445 m3s~! 


(10) Explain what is meant by the term equivalent head in the evaluation of 
energy losses across pipe fittings. 


(11) Water at a rate of 1 m3s~! and with viscosity 1x10-3 Nsm~? is to be trans- 
ferred from an open reservoir to another one 60 m below it through a smooth 
pipe 2400 m long. Determine the internal diameter of the pipe using the Blasius 
equation for smooth pipes. 


Answer: 0.48 m 


(12) A pipeline, 500 m long and with an internal diameter of 150 mm, 
connects two large storage tanks. If the difference of liquid level in the two 
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tanks is 30 m, determine the free flow through the pipe. Ignore entrance and exit 
losses but assume a friction factor of 0.01. 


Answer: 0.021 m3s-! 


(13) Draw a diagram, with a brief description, to illustrate how the friction 
factor varies with Reynolds number. 


(14) Describe the principle of the siphon and include a sketch. 
(15) Explain the limiting factors of a siphon used to transfer a liquid. 
(16) A four-side ring main ABDC is supplied with water at A for distribution in 


a process plant at points B, C and D at rates of 0.25 m?s-!, 0.1 m3s~! and 0.05 
ms-!, Determine the flowrate in pipe AB if the pipe details are: 


Pipe Length, m Diameter, m 
AB 1000 0.5 
BC 1500 0.3 
CD 500 0.3 
DA 1000 0.4 


The Fanning friction factor is assumed to be 0.005 in each of the pipes. 


Answer: 0.274 m*s-! 


(17) Determine the pressure drop along a 40 m length of pipe in terms of 
velocity head and equivalent length approaches for the flow of water at a rate of 
20 mh, if the pipe contains a gate valve (open) and three 90° elbows, The 
internal diameter of the pipe is 100 mm, the velocity head values for the valve 
and elbows are 0.15 and 0.7 velocity heads and the equivalent length of pipe are 
7 and 35 pipe diameters, respectively. The Fanning friction factor is 0.006. 


Answer: 2.96 kNm^? and 3.07 kNm? 
(18) Explain why calculations of the pressure drop due to friction for flow of a 


liquid through a pipe using the velocity head and equivalent length of pipe 
approaches do not necessarily provide the same result. 


(19) Water is supplied to a three-side ring main ABC with a flow of 0.1 m?s-! 


at A and is drawn from B and C at rates of 0.06 m3s~! and 0.04 m?s^, respec- 
tively, The internal diameters of the pipes are: AB 50 mm, AC 38 mm and BC 38 
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mm, and each has a length of 50 m. Determine the rate of flow in pipe BC if the 
friction factor is the same and constant for each pipe. 


Answer: 0.007 m3s-! 


(20) If, in the previous problem, the friction factor is given by the Blasius 
formula for smooth-walled pipes (f. - 0079Re ^), comment on the effect this 
would have and determine the flow in pipe BC. 


Answer: 0.0061 m3s~! 





(21) A fermented broth of density 990 kgm-? is decanted from a fermentation 
vessel to a receiving vessel by siphon. The siphon tube is 10 m long with a bore 
of 25.4 mm. Determine the rate of transfer if the difference in head between the 
two vessels is 2 m. Determine the pressure at the highest point of the siphon 
tube which is 1.5 m above the broth level and the tube length to that point is 4 m. 
Assume a friction factor of 0.004 and standard atmospheric pressure. 





Answer: 1.136x107? m*s-!, 74.3 kNm^? 


(22) A 0.155 m internal diameter pipeline is fed with water from a constant 
head tank system. The water is discharged from the pipeline 12.5 m below the 
level of the water in the head tank, The pipeline contains four 90° elbows, a 
fully open globe valve, a half open gate valve and a three quarters open gate 
valve. The respective equivalent length for elbow, open, half open and three 
quarters open valves are 40, 300, 200 and 40 pipe diameters, respectively. 
Determine the rate of flow if the Fanning friction factor is 0.025. 





Answer: 3.018x107? m: 


(23) Acylindrical tank with a diameter of 1.5 m mounted on its axis contains a 
liquid at a depth of 2 m. The liquid is discharged from the bottom of the tank 
through a 10 m length of pipe of internal diameter 50 mm at an elevation 1 m 
below the bottom of the tank. Flow is regulated by a valve. Determine the time 
to discharge half the liquid if the valve is open such that the energy loss through 
the valve is equivalent to 40 pipe diameters. The Fanning friction factor is 
assumed to be 0.008. 


Answer: 412 s 
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Pumps 


Introduction 

Pumps are machines for transporting fluids from one place to another, usually 
along pipelines. They are classified as centrifugal, axial, reciprocating and 
rotary and may be further grouped as either dynamic or positive displacement 
pumps. Dynamic pumps, which include centrifugal and axial pumps, operate 
by developing a high liquid velocity (kinetic energy) and converting it into 
pressure. To produce high rates of discharge. dynamic pumps operate at high 
speeds, although their optimal efficiency tends to be limited to a narrow range 
of flows. Positive displacement pumps operate by drawing liquid into a 
chamber or cylinder by the action of a piston; the liquid is then discharged in the 
required direction by the use of check valves. This results in a pulsed flow. 
Positive displacement pumps, however, are capable of delivering significantly 
higher heads than dynamic pumps. Rotary pumps are another form of positive 
displacement pump capable of delivering high heads; in this case fluids are 
transported between the teeth of rotating and closely meshing gears or rotors 
and the pump casing or stator. Unlike reciprocating pumps the flow is contin- 
uous, although rotary pumps tend to operate at lower speeds than dynamic 
pumps and their physical size is likely to be larger. 

In the decision-making process for the specification and selection of a pump 
to provide safe, reliable and efficient operation for a particular application, the 
general procedure follows many well-defined steps. Flow regulation and head 
requirements are the two most obvious fundamental factors that must be 
considered. Equally, however, the physical characteristics of the fluids — 
viscosity and lubrication properties, solids content and abrasiveness, as well as 
Corrosion and erosion properties — significantly influence the choice of pump 
type. A complete analysis of the system to discover the requirement for the 
number of pumps and their individual flow is also necessary. This analysis must 
include basic calculations for differential static head to provide an estimate of 
the necessary pump head. 
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Once an appropriate type of pump has been identified for a particular fluid, it 
is then necessary to establish the relationship between head and flowrate that 
can be delivered by the pump to match the system requirements in terms of 
static head and frictional losses. There are many additional calculations and 
checks that must be made before the analysis is complete. It is essential, for 
example, to ensure that pumps are correctly located to ensure avoidance of 
undesirable effects such as cavitation in centrifugal pumps and separation in 
reciprocating pumps. Clearly, not all pumping installations demand the same 
level of attention. Pumps required to deliver clean, cold water with modest 
delivery heads, for example, are less complex than those whose duties involve 
high temperature and pressure, high viscosity and abrasive liquids. 

In practice, there are numerous additional and important considerations that 
must be taken into account to ensure the safe, reliable and efficient installation 
and operation of a pump. It is important to consider maintenance requirements, 
implications and likelihood of leaks, availability of spares, economic factors 
such as capital and operating costs, safety for personnel and environmental 
considerations including noise. For the engineer new to pump specification and 
selection, the decision-making process may appear complicated, confusing and 
even, at times, conflicting. Experience and familiarity with pump selection 
eventually results in confidence, thereby reducing the effort required in the 
successful procurement, installation and operation of a pump. 
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9.1 Centrifugal pumps 
Describe the essential features and merits of centrifugal pumps and possible 
reasons for vibration during operation. 


Discharge 


j 







Volute 


Casing 


Impeller 
vane 


Impeller 


Solution 

Centrifugal pumps are by far the most commonly used type of pump in the 
process industries, because they are versatile and relatively inexpensive. They 
are very suitable for handling suspended solids, able to continue operating 
when the delivery line is blocked, have low maintenance costs and are easily 
fabricated in a wide range of corrosion-resistant materials. However, they are 
unable to develop high heads unless multiple stages are used and they also 
require priming by ancillary equipment. They offer reasonable efficiency over 
only a limited range of conditions and are not particularly suitable for very 
viscous fluids. 

The design of a centrifugal pump consists of a series of blades attached at the 
centre of a shaft known as an impeller and rotated at high speed inside a casing. 
Fluid is fed in axially at the centre (eye) of the impeller and is thrown out in a 
roughly radial direction by the centrifugal action. The large increase in kinetic 
energy which results is converted into pressure energy at the pump outlet either 
by using a volute chamber or a diffuser; the latter is more efficient but more 
expensive. There are considerable variations in impeller design, but almost all 
have blades which are curved, usually backward to the direction of rotation. 
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This arrangement gives the most stable flow characteristic. Centrifugal pumps 
do not operate by positive displacement and it is important to note that the head 
they develop depends not only on the size and rotational speed of the impeller 
but also on the volumetric flowrate. 

There are, in fact, many different types of centrifugal pump. Each is 
intended to perform a specific duty, such as handling various types of liquids 
including slurries, operating at high temperature or delivering high heads or 
flow rates. A variation of the single impeller pump shown is the two-stage 
centrifugal pump. This pump has two impellers mounted on the shaft such that 
the output from the first impeller is fed into the second impeller. By operating 
impellers in series, higher heads can be developed. Multi-stage centrifugal 
pumps have three or more impellers mounted in series on the same shaft and 
these pumps are therefore capable of producing considerably higher heads. 

In general, centrifugal pumps can be classified according to their construc- 
tion and layout in terms of impeller suction, type of volute, nozzle location, 
shaft orientation, bearing support and coupling to the driver. Single suction 
pumps, for example, have a suction cavity on one side of the impeller only, 
whereas double suction pumps have suction cavities at either side of the 
impeller. Single suction pumps, however, are subject to higher axial thrust 
imbalance due to flow coming in on one side of the impeller only. Pumps with a 
horizontal shaft are popular due to ease of servicing and maintenance, although 
pumps with a shaft in the vertical plane tend to be used where space is perhaps 
limited. Shafts that are unsupported by a bearing are known as overhung and 
can be fed exactly at the eye of the impeller. Shafts with bearing support on both 
ends such that the impeller is located in between the bearings, however, provide 
less shaft deflection although the shaft actually blocks the impeller eye. 

It is not uncommon for centrifugal pumps to be troubled by vibration during 
operation. This may be due to the destructive phenomenon of cavitation in 
which vapour is released from solution leading to a loss of discharge and deliv- 
ered head. The source of the vibration may, however, be due to operating the 
pump dry, against a blocked or closed delivery line. It may also be due to the 
shaft being out of alignment, an unbalanced impeller, bearing or seal wear, or 
the pump being incorrectly wired. Vibration may also be due to poor pipe 
support, the effects of solids, including foreign bodies, in the fluid, and air locks. 
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9.2 Centrifugal pump matching 

A centrifugal pump is used to transfer a liquid between two open storage tanks. 
A recycle loop mixes the contents of the feed tank and a restriction orifice in the 
recycle line is used to limit the liquid velocity to 2 ms~!. Determine the number 
of velocity heads across the restriction orifice when a valve in the transfer line 
to the receiving tank is closed. If the valve is fully open, show that the flowrate is 
approximately 1.11072 m?s7! when the level of the liquid in the receiving tank 
is 12 m above the centre line of the pump and the level in the feed tank is 5 m 
above the centre line. The pump characteristic is 


H =12-700 - 43000? 


All pipelines have an inside diameter of 100 mm and the suction, recycle and 
transfer lines have equivalent lengths of 2 m, 10 m and 20 m, respectively. A 
Fanning friction factor of 0.005 may be assumed. 


Not to scale 





‘Transfer 






Recycle - 
yi Restriction 


orifice 


Datum 





Suction 


Solution 
To determine the head loss across the restriction orifice consider the case in 
which the valve is closed. The rate of flow through the recycle line is therefore 





700157 n?s^! 
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where v, is the average velocity in the recycle line limited to 2 ms™! by the 
restriction orifice. 

From the pump characteristic, the delivered head by the centrifugal pump is 
thus 


H 212-709 43009? 
= 12 -70 x 0.0157 —4300 x 00157? 
=984 m 


At the duty point this is equal to the system characteristic. Applying the 
Bernoulli equation to the suction and recycle lines, the head loss is due to fric- 
tion, loss through the restriction orifice and exit loss from the pipe. That is 


SU. + Lr) +L VE yA 4 ONE 


H- 
d 2g 2g 2g 


where L, and L, are the suction and recycle equivalent lengths of pipes, 
respectively. Rearranging 


2gH  Af(L, * Ly) 


ks- rL] 
v2 d 
_ 2x gx 9844x0005 x(2 +10) _ | 
2 01 
2449 


That is, the head loss across the restriction orifice is equal to 44.9 velocity 
heads. 

For the case of the valve, being fully open allowing transfer of liquid to the 
receiving tank, the velocity in the transfer line, v;, is 


40, 
ad? 





v= 


_ 4x 1x10? 
nx 01° 


=14 ms! 


where Q, is the rate of transfer. Applying the Bernoulli equation to both 
suction and recycle lines, the system head is equal to the frictional head loss in 
the recycle line, suction line, across the restriction orifice and pipe exit. That is 
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H - Ae vr AfL, (v vp), v2 v? 
+ ee tol 
d 2g d 2g 2j 


Likewise, applying the Bernoulli equation to both the suction and transfer line, 
the system head is equal to the frictional head loss in the suction line, transfer 
line and pipe exit loss as well as the difference in static head between the liquid 
levels in the two tanks. That is 

M ons +i (v? tv) Hayr ve 


z +10 
` Og d 2g 2g 


Equating and rearranging, the velocity in the recycle line is therefore 


i 
a i 
2g(z, —Z5)+ v +1)? 
vp = 
JL, 
—ceb Koh 
3 1 
4x0. : 
2x gx (12-5) + ** mel x14? 
3 4 x 0005 
XIO donat 


01 


= 175 ms“! 
This therefore corresponds to a head of 


Ha Mle? AfL, (2 v2) 


HE aat 

d 2g d — 2g 2g 2g 
afari (Arth aye 
d 2g d 2g 


= 4x 0005 x2, 14? [eos 


175? 
01 2g 01 2, 


eun 


7749 m 
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The total flowrate delivered by the pump (through the transfer and recycle 
lines) is therefore 


Q =Q, +Q, 


2 
=o +v) 


_ 7x01? 





x (175 +14) 


= 00247 m’s™! 


which, from the pump characteristic equation, corresponds to a delivered pump 
head of 


H -12 -709 -43009? 
= 12 — 70 x 0.0247 — 4300 x 0.0247? 
=765 m 


The delivered head therefore closely matches the system head at the flowrate of 
1.1x1072 m?s-!, corresponding to the duty point. 

It should be noted that while the duty point corresponds to the maximum 
flow of liquid delivered by the pump for the system illustrated, it may not corre- 
spond to the most efficient operation of the pump (best efficiency point, or 
‘bep’). Also, to ensure that cavitation is avoided, details of the pump’s required 
net positive suction head (NPSHR) are also required (see Problem 9.5, page 
245). 
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9.3 Centrifugal pumps in series and parallel 
Highlight the benefits of using centrifugal pumps in series and parallel. 


Solution 
There are occasions when it is desirable to use centrifugal pumps to deliver a 
higher head than is normally possible for a fixed flowrate. This can be achieved 
by using two pumps in series. The total pump characteristic is the sum of the 
individual characteristics. Pumps in series are often used for high pressure 
applications such as boiler feed pumps. Special designs are often used — for 
example, running the impellers on a common shaft with the casing laid out so 
that the volute of one pump leads into the eye of the next. 

It is possible to deliver high flowrates for a given head by using similar 
pumps in a parallel arrangement where the pump characteristics are added 
together. 


Head 





Parallel ^ Q, 


Pump 1 


Head 


Pump 1+ Pump 2 











Flow 





Pump 1+ Pump 2 





Flow 
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9.4 Cavitation in centrifugal pumps 


Explain what is meant by cavitation in centrifugal pumps. 


Solution 

Vapour forms in any liquid when the pressure in the liquid is iess than the 
vapour pressure at the liquid temperature. The possibility of this happening is 
much greater when the liquid is in motion, particularly on the suction side of 
centrifugal pumps where velocities may be high and the pressure correspond- 
ingly reduced, Having been formed, the vapour bubbles travel with the liquid 
and eventually collapse with explosive force giving pressure waves of high 
intensity. This collapse, or cavitation, occurs on the impeller blades causing 
noise as well as vibration and erosion of the blades, which may eventually 
result in a typically pitted appearance similar to that of corrosion. Apart from 
the audible and destructive effect of cavitation, another sign is a rapid decrease 
in delivered head and pump efficiency. As a remedy, a throttle valve should be 
placed in the delivery line and when any symptoms of cavitation are observed 
the valve can be partially closed, thereby restricting the throughput and 
reducing the kinetic energy. 

Cavitation is more likely to occur with high-speed pumps, hot liquids and 
liquids with a high volatility (high vapour pressure). Problems can also be 
encountered with slurries and liquids with dissolved gas, and in particular in gas 
scrubbers which usually operate with liquids saturated with gas. Note that the 
related phenomenon of separation can occur in reciprocating pumps due to the 
acceleration of the liquid. This is more likely when the delivery lines are long 
and the corresponding kinetic energy is large. Apart from providing sufficient 
net positive suction head (see Problem 9.5, page 245), separation can be 
avoided by attaching air vessels to both the delivery and suction lines. 
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9.5 Net positive suction head: definition 


Explain what is meant by available and required net positive suction head 
(NPSH). 





Discharge 






Suction 


Solution 

To avoid problems of cavitation in centrifugal pumps it is necessary for the 
lowest pressure of the liquid in the pump, which is usually at the eye of the 
impeller, to exceed the vapour pressure of the liquid being delivered. Should 
the pressure fall below the vapour pressure of the liquid, local vaporization 
(boiling) is likely to occur although the precise mechanism of cavitation 
inception is not fully understood. Since the vapour pressure of liquids is a 
function of temperature, it may be the case that cavitation is more likely to 
occur during the hot months of summer than in cold winters. Atmospheric 
conditions may also affect the likelihood of cavitation occurring at low 
barometric pressures. The required net positive suction head (NPSH) or 
NPSHR is therefore the positive head necessary to overcome the pressure drop 
in the pump and maintain the liquid above its vapour pressure. It is a function of 
pump design, impeller size and speed, and the rate of delivery. To prevent cavi- 
tation from occurring this must be exceeded by the available NPSH or NPSHA 
and is a function of the system in which the pump operates. This includes the 
minimum working gauge pressure of the vapour on the surface of the liquids p, 
the minimum atmospheric pressure taken as 0.94 of standard atmospheric pres- 
sure p,,, the vapour pressure of the liquid at the maximum operating tempera- 
ture p, , the static head above the pump centre H, and the head loss due to 
friction and fittings H ; . The available NPSH or NPSHA at the pump suction is 
then 
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NPSHA -P*Pa Pv 4g Lg 
Ps 
and is effectively a measure of the total head available at the pump suction 
above the vapour pressure. To ensure that vapour bubble formation at the eye of 
the impeller — and thus problems of cavitation — are avoided, it is essential 
that the NPSHA is greater than the NPSHR. That is 


NPSHA » NPSHR 


In the case of reciprocating type positive displacement pumps, the same 
approach to calculating the available NPSH is used. This allows for the static 
liquid head at the pump inlet, minimum barometric pressure, minimum 
working pressure and vapour pressure of the liquid at the maximum operating 
temperature. It should be noted, however, that the maximum frictional head in 
the pump suction line occurs at the point of peak instantaneous flow and not the 
average flow. This is dependent on the configuration of the pump and has 
greatest effect for single cylinder (or simplex) pumps (see Problem 9.15, page 
265). Additionally, a head correction is required to allow for the acceleration of 
the fluid in the suction line during each pulsation cycle. 

As with reciprocating pumps, the available NPSH for rotary pumps which 
operate with pulsed flow requires correction for both frictional head loss in the 
suction line at peak instantaneous flow and head required to accelerate the fluid 
in the suction line during each pulsation cycle. 

Note that while the available NPSH is a function of the system in which the 
pump operates, the vapour pressure of the liquid being pumped is not always 
specifically included in the NSPHA calculation. In this case the vapour pres- 
sure head is included such that 


NPSHA > NPSHR + 2" 
pg 


Details of the required NPSH are usually supplied by the pump manufacturer. 
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9.6 Net positive suction head: calculation 1 

A centrifugal pump is used to deliver a liquid of density 970 kgm- from an 
open storage vessel at a rate of 5 mh. The storage vessel has a diameter of 3 
m and is initially at a depth of 2.5 m. The pump is located at an elevation of 3 m 
above the bottom of the vessel and the frictional head loss in the suction pipe is 
0.5 m. The vapour pressure of the liquid at the temperature of operation is 18. 
kNm^? and the NPSH is 5 m. Determine the quantity of liquid delivered and the 
time taken before cavitation occurs. Allow for worst case meteorological 
conditions. 





Solution 
Applying the Bernoulli equation between the free surface of the liquid (1) and 
the suction point (subscript s) in the pump 


2 


1 
SS Pr ep ey He aE: 
sg 





The minimum suction head before cavitation occurs is 


2 
Ps A NpsH- YE 4 Pe 
pe 2g pg 
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Thus, for cavitation not to occur 


y2 


2 
NPSH 284 Pe -P1 42,2, 8 a, 
28 ps ps 2g 
Rearranging, and noting that the pressure on the free surface is taken as 0.94 of 
the standard atmospheric pressure, the depth that can be pumped before cavita- 
tion occurs is 


Py Pi 


zı = NPSH + PI tz, +H, 
Pg 
3 
= 5 4 (8 -094 x 101.3) x 10 43405 
970x g 
=0.375 m 


corresponding to a volume of 


v ix 3? x (25 -0.375) 


= 1502 m? 


The time taken to reach this level is therefore 


The quantity of liquid delivered is therefore 15 m3, taking 3 hours before cavi- 
tation occurs. 
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9.7 Specific speed 


Explain what is meant by specific speed of a centrifugal pump. 


Solution 

The specific speed is useful for the selection and scale-up of centrifugal pumps 
and is a measure of pump performance in terms of pump discharge Q, head H 
and impeller speed N. It is always evaluated at the point of maximum effi- 
ciency, otherwise known as the best efficiency point. The general procedure is 
that once the head and the flow are established for a particular duty, the pump's 
specific speed can be determined to ensure the selection of the appropriate 
pump in terms of impeller type and characteristics with optimal operating effi- 
ciency. This is determined for a head coefficient C;; and capacity coefficient 
Co. which, for two geometrically similar pumps, are 


28H, 
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NID; 


(see Problem 4.3, page 104) 





From the head coefficients, the ratio of pump impellers is 


; 
Di Naz BUY 
Dı N\\H2 


and from the capacity coefficients, the ratio of discharges is 


9...) [Di 
Qj N2(D; 
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The ratio of discharges is therefore 
1 
Dien Ne (ee US 
Q2 N2|Ni\ Ao 


3 
-(ne Tp Y 
N, J (H2 
or rearranging 


L i 
NiO; _ N20; 


3 3 
HE. (MS 





and is called the specific speed, N,. Although the pump specific speed is 
usually described as being dimensionless, it has the dimensions of 247-32, 
The numerical value of the specific speed corresponds to the type of pump for 
required conditions of rotational speed, flowrate and head. As a guide 
N,«0.36 Multi-stage centrifugal/positive displacement 

0.36«N , «1.10 Single-stage centrifugal (radial) 

1.10€N , «3.60 Mixed flow pumps 

3.60<N , «5.50 Axial flow pumps 


The maximum efficiency for each pump type spans a narrow range of specific 
speeds with little overlap between pump types. The specific speed is therefore 
valuable to identify the pump type for a particular application. 








* 
- 
E 
5 
i 
a 
Maximum efficiency |. Maximum efficiency Maximum efficiency 
range of centrifugal | range of mixed-flow range of axial-flow 
pumps pumps pumps 
- 
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9.8 Net positive suction head: calculation 2 


A low level alarm is to be installed in a vessel containing warm water at a 
gauge pressure of 50 kNm'? to prevent cavitation at the design temperature of 
40*C. Determine the minimum level of the alarm above the centre line of the 
pump if separate tests show that the minimum NPSH for the pump is given by 


4 
NPSH-28N?H 


where N , is the specific speed of the pump at maximum efficiency and has a 
value of 0.14 (m? s3? ), and H is the differential head and is 30 m. The fric- 
tional losses in the suction line amount to 0.2 m. 


Solution 
The minimum net positive suction head is 


4 
NPSH = 28N3H 
4 
=28x 0145 x 30 
=611m 
The NPSH is given by 
NPSH =?1—P» 42) H, 
Ps 
where p,, is the vapour pressure of the water at 40°C (7380 Nm-?) and density, 
p. is 992 kgm^3. Rearranging 


zı = NPSH + PL ug 
pg 


5 p 
26114 038 50)x 10: +02 
992x g 


=193 m 


The position of the low level alarm is 1.93 m above the centre line of the pump. 
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9.9 Effect of reduced speed on pump characteristic 
The characteristic of a centrifugal pump is given by a manufacturer in dimen- 
sional form as 


H =40 -140Q — 42000? 


where Q is the flowrate (m?s-! ) and H is the head (m). Obtain an expression for 
the pump characteristic for the same pump operating at 75% of its test speed. 


Solution 
A centrifugal pump can operate over a range of flowrates corresponding to a 
range of delivered head, known as the pump characteristic. The shape of the 
curve depends on the layout of the impeller, blades and volute. In this case, the 
specific speed of the pump at full operating, N, and reduced speed, Nj, is 
N; =075N 
The head coefficient for the pump in both cases is therefore 
gH 
N?p? 
aH, 
252 
NID 


Cu = 





(see Problem 4.3, page 104) 





Therefore 
HN;D? 
pec 
Nip? 





_ HOISNY 
(0 Np? 
-05625H 
The capacity coefficient is 
edP 
ND? 
Qi 
ND? 


Co 
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Therefore 
ND? 
o, - 2D 
ND 
_ Q(075NyD* 
ND? 
-075Q 


The pump characteristic for the test conditions is therefore 


A, Qı Qi 
— 540 -140x 2L -4200x| 2L 
05625 075 "logs 


That is 
H, =225 -1050, - 42000? 


The pump characteristic for both full operating and test conditions is therefore 
DischargeQ. ms! — 0 0.001 0.02 0.03 0.04 0.05 0.06 0.08 
Diff. head-H, m 40.00 38.18 35.52 32.02 27.68 22.50 17.88 1.92 
Diff. head H,, m 22.50 21.03 18.72 15.57 11.58 6.75 2.13 - 


Both delivered head and discharge are therefore reduced at the lower 
impeller speed. Centrifugal pumps usually operate with a fixed speed. 
However, mechanically worn pumps may operate at reduced speeds and may 
therefore be unable to perform their required duty. 


HEAD (m) 





0 0.02 0.04 0.06 0.08 
FLOW (ms!) 
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9.10 Duty point and reduced speed of a 


centrifugal pump 
A centrifugal pump with an impeller speed of 1200 rpm has the following 
characteristic: 


Discharge Q, m?^s-! 0.000 0.002 0.004 0.006 0.008 0.010 
Head H, m 400 395 3&0 350 300 200 


The pump is used to deliver a process liquid along a pipe of length 100 m and 
inside diameter 50 mm and discharges from the end of the pipe at atmospheric 
pressure at an elevation of 10 m above the open feed tank. Determine the duty 
point and the speed of the pump that would result in a reduction in flow of 25%. 
Assume a friction factor of 0.005 for the pipe and neglect entrance and exit 
losses. 


Solution 

The duty point is the maximum flowrate that the pump can achieve for a 
particular power demand and forms the intersection between the pump and 
system characteristic curve (a plot of discharge against head). The system char- 
acteristic is usually parabolic, since the frictional head loss is proportional to 
the square of the fluid flowrate. Applying the Bernoulli equation, the head to be 
developed by the pump is 


2 
Ho +z -z +H, 


2g 

2 2 
ots oe 

2g d 2g 


or in terms of flowrate 


160? (agn 
H=22 -z, +— |= 
: ees i 


d 
iig. S5 4x0005x100 |, 
g x i x 005* 005 
-10-- 542 x 109 9? 
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The pump and system characteristic are therefore 

DischargeQ, m?s-! — 0.000 0.002 0.004 0.006 0.008 0.010 
Head (pump),m 400 395 380 350 300 200 
Head (system),m 10.0 122 187 295 447 642 


40 System 


34 


30 
Pump 
2) 


HEAD (m) 


10 


0 0.002 0.004 0.006 0.008 0,010 
FLOW (ms!) 


From the graph, the duty point corresponds to a flowrate of 0.0066 m?s7! and a 
head of 34 m. A 2556 reduction in flow is therefore 0.00495 m?s7! with a head 
of 23 m. Since the head coefficient for the pump at the two speeds is 


(see Problem 4.3, page 104) 





then the reduced speed is 


1 
HY 
Np =N,| 22 


i 
= 1200x{ 23 
34 


= 987 rpm 
The speed of the pump to reduce the flow by 25% is 987 rpm. 
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9.11 Power, impeller diameter, speed and 
delivered head 


A centrifugal pump is used to transfer water from a vessel at an absolute pres- 
sure of 50 kNm-? to another vessel at an absolute pressure of 150 kNm-2. The 
elevation of the water in each vessel above the centre line of the pump is 5 m 
and 15 m, respectively. The connecting pipework on the suction and delivery 
sides of the pump has equivalent lengths of 5 m and 20 m, respectively. In addi- 
tion, there is a control valve in the line which has a loss of 25% of the pump 
differential head; the flow controller maintains the flow at a rate of 0.05 m?s-!. 
Data from a similar geometry pump with an impeller diameter of 0.15 m and a 
rotational speed of 1500 rpm are 


Discharge Q, Is! 0 5 10 15 20 25 
Differential head Hm || 9.25. — &81 7.85 648 481 296 
Power input P, kW - 0.996 103 119 1.26 145 


For a liquid velocity of 2 ms~!, determine the differential head, the impeller 
diameter and rotational speed required for operation at maximum efficiency, 
and the power input to the pump. Assume a friction factor of 0.005. 


























Solution 
Applying the Bernoulli equation at point 1 and point 2 


PL vi +H =P2 
Ps ps 








= PUMPS 


where H; is the loss across the valve and frictional loss in the pipeline for 
which the diameter is 


e) 


 (4x005 Y 
nx2 


70178 m 





The suction line losses are 


AfL, v? 
d 2g 





Hy(s) = 


4x 0005 x5 _ 22 
SES a 


0178 2g 
=0114m 


and the delivery line losses are 


LAX0005x 20 . 
0178 2g 


= 0458 m 





Rearranging the Bernoulli equation 


A 
- v3 
H-P P2 uus + Hye) + Hyay +025 


ps 2g 
The differential head is therefore 


2 
| [pni "2 
= ———— veo gy dd yia T 
1-025| ps 2g 5779 rosy FAG) 





S 2 
E g|üs-0speloS 2c (5-6 d 0TTASE 0458 


a 1-025 1000 x g 2g 





-2796 m 
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The actual power input for the model pump is 


p -PSOH 
n 


Rearranging, the efficiency is therefore 


_ PsQH 
Venu 


and calculated from the pump data as 


Flowrate Q, ls-! 0 5 10 15 20 25 
Differential head H,m 9.25 8.81 785 648 4.81 2.96 
Power input P, kW - 0.96 1.03 1.19 1.26 145 
Efficiency n, 96 0 45 75 80 75 50 


The maximum efficiency is therefore 80% for which the speed is determined 
from the pump head and capacity coefficients (see Problem 9.7, page 249) 


E 


3 L 
= 1500 „ (2796 * (0015 Y 
60 648 005 

=41 rps 


The diameter of the pump impeller is therefore 


i 
s 
D=D, QN, 
QN 
i 
005 x 1300 * 
z2015x| ———À— 
0015 x 41 
-019m 
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The actual power is therefore 
p PSQH 
n 
 1000x g x 005x 2796 
080 
217143 x 10 W 


The power input is found to be 17.1 kW for a pump with an impeller of 19 cm 
and rotational speed of 2460 rpm. In practice, the NPSH calculation should also 
be completed to ensure cavitation is avoided. 

Note that at maximum efficiency, the specific speed is 


1 
N,Q? 
n, TEL 


By 


1500 


i 
x0015° 


i 


2 
6484 


" 


075 (mS 45-32) 


corresponding to a radial single-stage centrifugal pump (see Problem 9.7, page 
249). 
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9.12 Suction specific speed 

A centrifugal pump is to deliver 50 m?hr-! ofa process liquid. If the total head at 
the best efficiency point for the pump is 5 m, determine the suction specific 
speed as a dimensionless number for an impeller speed of 2000 rpm. 


Solution 

It has been the practice to define the specific speed for a pump in terms of its 
cavitation characteristics so that an impeller diameter D, rotating with a speed 
N, and delivery volumetric flowrate Q, at optimal efficiency would produce 
the same pump inlet conditions for a geometrically similar pump under similar 
operating conditions. Thus, the head coefficient Cj; and capacity coefficient 
Cp for the two geometrically similar pumps are 


gH 
Nip? 





(see Problem 4.3, page 104) 


ND? 


Eliminating D and rearranging in terms of the impeller speed 
L 3 
Ch (8H )* 
N=2 31 
402 
cho 


Defining the suction specific speed S,, as the dimensionless group 





where N is the impeller speed (rps),Q is the flow through the pump and H is the 
net positive suction head, then 
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The suction specific speed provides an indication of impeller specification as 
well as the likelihood of cavitation where as a guide 


$,«0.12 Pumps are free from cavitation. This is a useful limit 
for pumps used to deliver very corrosive gas-free liquids 
where materials of construction which are resistant to the 
destructive effects of cavitation cannot be used 
(see Problem 9.4, page 244). 


0.12<S,, <0.4 The pump inlet bore may be the same as the bore of the 
pump inlet connection. An inlet reducer is therefore not 
required, It is also recommended that bends connected to 
the inlet have a radius in excess of five times the inlet bore. 


0.4<S,, <0.7 Materials of construction should be resistant to corrosion 
as pumps operate in this range with incipient cavitation, 
This range is typical for commercial pumps handling clean 
liquids. 


$,20.7 Applies only to pumps which are used to handle clean non- 
corrosive liquids. As cavitation is likely, it is essential that 
pump suction calculations and pump layout are examined 
closely. Inlet pipes should be straight and free from 
obstructions, and inlet isolation valves (such as ball 
valves) should not restrict flow. 
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9.13 Reciprocating pumps 
Describe, with the aid of a diagram, the operation of reciprocating pumps, 
their applications and limitations. 


| Discharge pipe. 


Pressure. 
valve 


YZ Piston te 


Suction 
valve 






fg 























Suction 


| pipe 





Solution 

In general, reciprocating pumps transport fluid by the action of a piston on the 
fluid contained within a cylinder. The piston is provided with a reciprocating 
motion by means of a connecting rod and crank where the discharge is 
dependent on the swept volume and the stroke frequency; the swept volume is 
the amount of fluid displaced on each stroke depending on the cross-sectional 
arca of the cylinder and the length of stroke. Check valves are needed to ensure 
that flow is in the correct direction. 

Although they account for a relatively small proportion of all the pumps 
used on process plants, reciprocating pumps are applied universally across a 
wide range of industries. They are particularly useful for delivering fluids into 
high pressure process streams. Indeed, reciprocating pumps are capable of 
delivering heads up to several thousand bar and are limited only by the 
mechanical strength of the pump and slight leakage. They are also self-priming. 
The flow is not constant but pulsed, although the use of gas (usually nitrogen) 
vessels and other pulsation dampers smoothes the flow. Reciprocating pumps 
provide unsatisfactory performance with very viscous liquids due to slow 
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check valve action, and with suspended solids due to attrition between the 
piston and cylinder. 

The volumetric efficiencies of reciprocating pumps are generally high and 
dependent on the stroke frequency and pressure — limited largely by check 
valve leakage caused by wear, corrosion or mis-seating due to trapped 
particulates. Another limitation is the phenomenon of separation in which the 
pressure at the face of the piston falls below the vapour pressure of the liquid. 
The liquid therefore becomes separated from the piston by vapour. Further 
problems can arise involving vapour or gas which enters the cylinder space 
from around a leaking piston seal. This results in a fall in the rate of discharge. 
audible noise and damage to the pump. 

In addition to air vessels on both the delivery and suction lines, the flow 
characteristic can be made smoother by the use of multi-cylinder pumps. The 
total flow is thus the sum of the instantaneous flows from the cylinders where 
odd numbers of cylinders (three or triplex, or five or quintuplex) produce the 
smallest pulsations (see Problem 9.15, page 265). 

Abrasive, toxic and corrosive liquids, including slurries and radioactive 
fluids, can be handled by a special type of reciprocating pump, the diaphragm 
pump, which consists of a flexible diaphragm fixed at the edges and flexed to 
and fro, The diaphragm pump, however, requires priming prior to operation. 

Another variation is the cylinder or ram pump in which a plunger replaces 
the piston and has a clearance between the cylinder wall. It is more expensive 
than the piston pump although leakage is easier to detect, the packing is easier 
to replace and it is more suitable for delivering high heads. 
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9.14 Single-acting reciprocating piston 

A single-acting reciprocating piston has a piston area of 100 cn? and a stroke 
of 30 cm and is used to lift water a total height of 15 m. If the speed of the piston 
is 60 rpm and the actual quantity of water lifted is 10.6 mh”, determine the 
coefficient of discharge and the theoretical power required to drive the pump. 


Solution 
The expected volume displaced in the piston cylinder is 


V = ALN 
-001x03x © 
60 


23x10? ms! 


The actual volume discharged is 


. 106 


act ^ 3600 
2294x107 ms! 


The coefficient of discharge is therefore 


Cy = Vact. 


V 


| 294x107 
3x10 
- 098 


and the power required is 
P - pgQH 


=1000x g x 18 515 
3600 


=433 W 


The coefficient of discharge and theoretical power required are 0.98 and 433 
W, respectively. 
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9.15 Discharge from reciprocating pumps 
Determine the ratio of peak to average discharge from single, two, three, four 
and five-cylinder reciprocating pumps. 





Solution 
Itis assumed that the piston in each of these pumps operates with perfect simple 
harmonic motion and that the volumetric efficiency and acceleration effects of 
the fluid are ignored. Thus, for a single-cylinder pump, the displacement of 
liquid is 
x-r-rcos 

= r(l -cos 8) 

= ríl -cos wr) 
The velocity of the piston is 

_ de 


v=o 


dt 
= or sin OF 


= wrsin ð 


The volume displaced on the pressure stroke (between 0 and m) is therefore 
x 
V - Aor [sin 640 
0 


= -Awr cos(n — 0) 
=2A0r 


where A is the area of the piston. The average flow over the cycle (between 0 
and 27) is therefore 
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av 
Q= 2n 11050, 
Triplex A 
_ 2Aar 9510, 
EET 
_ Aor 
x 
Since the peak flow occurs at 7/2, then 
x 1110, 
Qpeak = Aor sin z Quadruplex p 
= Awr 9.790, 
The ratio of peak to average flow is therefore 
Q peak . Aor 
Q Aor 
n 
024 
eu Quintuplex . o s 
0.960, 
= 314 
Likewise, for duplex (two-cylinder), triplex (three-cylinder), quadruplex 
(four-cylinder) and quintruplex (five-cylinder) pumps, the ratios are found to 
be 1.57, 1.05, 1.11 and 1.02, respectively. 
360 
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9.16 Rotary pumps 


Describe briefly, with a labelled diagram, the layout and operation of rotary pumps. 





Rotary pump (gear type) 


Solution 

In rotary pumps, portions of liquid are carried from the entry port to the delivery 
port in compartments within the pump formed between the rotor and the stator. 
They are manufactured to a very close tolerance to reduce back-leakage. They 
are therefore expensive to fabricate. The flowrate depends on the size and 
number of compartments as well as the frequency of rotation. They have the 
capability of delivering against high heads, are self-priming and provide 
smooth and continuous flow. They offer higher speeds than reciprocating 
pumps and are suitable for viscous liquids since no check valves are required. 
But they require careful maintenance and are unsuitable for liquids without 
lubricating properties such as water, and for slurries due to small clearances and 
possibility of friction. 

There are various forms of rotary pump. The Mono pump, for example, has a 
helical worm rotating inside a similarly shaped flexible stator. This type of 
pump can handle corrosive and gritty liquids, and is commonly used to feed 
slurries to filter presses, but it must never be allowed to run dry. The Roots 
blower is used for transporting large volumes of gas at low gauge pressure, 
using rotors in the form of cycloidal lobes. 

Note that it is essential to provide all positive displacement pumps (recipro- 
cating and rotary types) with a safety relief valve to avoid excessive pressure 
build-up in the event of the delivery line being inadvertently closed or blocked. 
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Further problems 
(1) Describe briefly, with the use of clearly labelled sketches, the essential 
features, uses and limitations of rotary pumps. 


(2) Explain the phenomenon of cavitation applied to centrifugal pumps and 
describe practical measures which can lead to its reduction. 


(3) Construct a flowchart for the decision-making process for the selection of 
centrifugal pumps. 


(4) Determine the additional available NPSH of a system for a pump used to 
transfer water that allows for worst case meteorological conditions. 


Answer: 0.62 m 
(5) Describe the essential features of the centrifugal pump; include a sketch. 


(6) A centrifugal pump is used to deliver a dilute solution of sodium 
hydroxide with a density of 1050 kgm-? from a large open buffer tank to a 
header tank 10 m above it. The velocity of the solution in the suction pipe to the 
pump is 1.5 ms~!. The suction pipe has an internal diameter of 3 cm and the 
delivery pipe has an internal diameter of 4 cm. If the frictional head loss is 1.5 
m, determine the theoretical power requirement for the pump. 


Answer: 126 W 
(7) Explain the significance of net positive suction head (NPSH). 


(8) Sulphuric acid with a specific gravity of 1.83 and a viscosity 0.04 Nsm-? 
is pumped into a large vessel through a 25 cm inside diameter pipe 30 m long 
inclined at an angle of 30° to the horizontal. Determine the power required to 
pump 0.064 m?s-! into the vessel. A friction factor of 0.008 may be assumed. 


Answer: 17.6 kW 


(9) Sulphuric acid with a specific gravity of 1.84 is to be pumped from an 
open tank to a process column at a rate of 25 kgs-!. The column operates at a 
pressure of 125 kNm-? and the acid is sprayed from a nozzle located 20 m 
above the acid level in the tank at a linear velocity of 3 ms~!. If the total addi- 
tional losses are equivalent to 2.75 m of water head, determine the power 
required to pump the acid if the pump has an efficiency of 65%. 


Answer: 9.01 KW 
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(10) 75% sulphuric acid, of density 1650 kgm-? and viscosity 8.6x 103 
Nsm^?, is pumped at a rate of 10000 kgh-! through an 800 m section of smooth 
pipe with an internal diameter of 50 mm with a rise of 12 m. If the pump effi- 
ciency is 50%, determine the power requirement. Assume the Blasius equation 
applies. 


Answer: 1.74 kW. 


(11) Liquid, of density 850 kgm^2, is pumped at a rate of 10 kgs-! from an 
open storage tank at atmospheric pressure to a distillation column which oper- 
ates at a pressure of 500 kNm-? above atmospheric pressure. The inside diam- 
eter of the pipe supplying the column is 10 cm, and the surface of the liquid in 
the tank is 5 m above the ground level and can be considered constant. The 
column feed point is 24 m above the ground. Determine the power required for 
the pump if the overall efficiency of pump and motor is 70%. Losses due to fric- 
tion amount to 3 m. 





Answer: 11.5 kW 


(12) The overhead condensate from a distillation column of density 600 
kgm and viscosity 0.6x107* Nsm-? is fed at a rate of 5 kgs~! from the conden- 
sate drum held at a gauge pressure of 12 bar to a pressurized storage vessel at a 
gauge pressure of 10.2 bar prior to further processing. The pipeline has a length 
of 2.5 km and an inside diameter of 7.5 cm with a wall roughness of 0.046 mm. 
If the difference in level of the condensate in the condensate drum and the 
discharge of the pipe into the storage vessel is 1 m, determine the power 
requirement of the pump and motor assembly if their combined efficiency is 
65%. Allow for entrance and exit losses and where the pipe contains fifty-six 
90° elbows of 0.7 velocity heads each and four open gate valves of 0.15 velocity 
heads each. 





Answer: 7.3 kW 


(13) Ethane is removed from an LPG mixture in a distillation column in which 
the ethane-rich overhead vapour is condensed, and some is returned to the top 
of the column as reflux. The top of the column operates at a gauge pressure of 
12 bar. The gauge pressure in the reflux drum is 11.5 bar. The difference in 
elevation between the top tray of the column and the level of liquid in the reflux 
drum is 10 m. The connecting pipework consists of 50 m of 0.254 m inside 
diameter pipe with 27 bends, two gate valves and one control valve. If the 
flowrate of refluxed liquid, of density 600 kgm-2, is 299 m?h-', determine the 
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power requirement of the reflux pump and motor assembly assuming an effi- 
ciency of 75%. Allow for entrance and exit losses where each bend and valve 
has a loss equivalent to 0.7 and 0.15 velocity heads, respectively. The constant 
friction factor may be taken as 0.006. 


Answer: 14.4 kW 





(14) An LPG mixture, of density 600 kgm-? and viscosity 6x107* Nsm-2, is 
pumped at a rate of 10 kgs“! from a storage vessel held at a gauge pressure of 
500 kNm~ to a distillation column for separation. The pipeline is 50 m in 
length and has a bore of 10 cm and surface roughness of 0.046 mm. If the gauge 
pressure at the feed tray in the column is 1.1 MNm™ and is located 20 m above 
the centre line of the pump while the surface of the LPG mixture in the storage 
vessel is 3 m above the pump, determine the power requirement of the pump 
and motor assembly assuming an efficiency of 60%. Allow for entrance and 
exit losses; the pipe has eight 90° elbows and three open gate valves of 0.7 and 
0.15 velocity heads each, respectively. 





Answer: 20.2 kW 


(15). Anethanol distillation column operating at a gauge pressure of 50 kNm-2 
receives a liquid feed of density 950 kgm-? and viscosity 9x10-4 Nsm-? at a 
rate of 39 m?h-! from an open storage vessel. The connecting pipework 
consists of 50 m of 10 cm inside diameter pipe with an absolute roughness of 
0.04 mm. The pipe has a total of 12 elbows and three gate valves. The level of 
liquid feed in the storage tank is steady at 4 m above the centre line of the pump 
and the feed point on the distillation column is at an elevation of 5 m above the 
pump. If the overall efficiency of the pump and motor is 70%, determine the 
power requirement of the pump. Allow for entrance and exit losses and where 
the loss for each bend and valve is equivalent to 0.7 and 0.12 velocity heads, 
respectively. 





Answer: 1.2 kW 


(16) A methanol-water liquid mixture is to be transported from a storage tank 
to a distillation column for separation. The mixture, of density 895 kgm", is 
stored within an enclosed tank maintained at an absolute pressure of 1 atmo- 
sphere and is to be pumped to the column at a rate of 11 kgs-! through a pipe of 
inside diameter 12 cm. The column, which operates at a gauge pressure of 50 
kNm-?, has a feed point 15 m above the level of the liquid mixture in the storage 
tank. If the overall efficiency of the pump and motor is 65% and the losses due 
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to friction amount to 4 m, determine the power requirement of the pump. 


Answer: 4.1 kW 


(17) Show that the work exerted on a fluid by a centrifugal pump is given by 


NE 
53.3 


where W is the work (Watts), N is the speed of the impeller (revolutions per 
minute) and F is the force (Newtons) on a torque arm measured from the centre 
line of the shaft which has a radius of 179 mm. 


(18) Experimental data were obtained from laboratory tests of a centrifugal 
pump operating at 2500 rpm using water as the test fluid. F is the force exerted 
on a torque arm measured from the centre line of the shaft which has a radius of 
179 mm. 


Flowrate,Q (x10-2 m?s7!) — 1.47. 1.31. L18 1.00 0.9] 0.7] 0.48 
Diff. pressure, Ap (kNm?) 18.8 36.0 44.8 61.5 63.6 71.5 797 
Force on arm, F (N) 46 41 40 37 33 30 24 


Determine the flowrate delivered at maximum pump efficiency. Comment on 
the general performance of the pump. 


(19) Water is to be transported from one tank to another using a 10 m pipeline 
of inside diameter 25 mm. A centrifugal pump is to be used to transport the 
water and is positioned below the first tank whose free surface is 2 m above the 
freely discharging outlet. Determine the duty point for the pump characteristic 
given below. The friction factor may be assumed to be 0.005. Allow for 
entrance and exit losses. 


Discharge, m?min-! 0.000 0.008 0.012 0.018 0.024 
Head, m 35 32 30 24 14 


Answer: 0.02 m?min-! 
(20) Acentrifugal pump is used to transfer a liquid of density 790 kgm? from 
an open storage vessel. The centre of the pump is located at an elevation of 2 m 


above the liquid in the vessel. Determine the available NPSH if the vapour pres- 
sure of the liquid at the temperature of operation is 25 kNm-? and the frictional 
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head in the suction pipe is 0.3 m. Allow for worst case meteorological 
conditions. 


Answer: 6.77 m 


(21) A liquid hydrocarbon liquid mixture, of density 538 kgm-? and viscosity 
3.04x10-* Nsm-2, is pumped from the bottom of a distillation column oper- 
ating at 120 KNm^? through a reboiler at a rate of 14x10? kgh-!. The centre line 
of the pump is a vertical distance H below the column take-off point and the 
suction line has a total length of H +3 metres and an inside diameter of 65 mm. 
The reboiler is at an elevation of H +0.5 metres above the pump. If the pressure 
drop through the reboiler is 14 kNm-, determine the differential head across 
the pump. Neglect kinetic energy head loss and losses due to fittings, and 
entrance and exit losses. Assume a constant Fanning friction factor of 0.004, 
The NPSH for the pump is 3 m. 





Answer: 30.11 m 


(22) A centrifugal pump delivers water at a rate of 0.02 m3s~! against a head 
of 10 m. If the impeller speed is 1500 rpm, determine the specific speed Ny. 


Answer: 0.63 m34s-32 


(23) Determine the suction specific speed S, for a centrifugal pump operating 
at an impeller speed of 40 rps delivering a flow of 110 m?h7! against a head of. 
14 m. 


Answer: 0.174 


(24) Determine the numerical value of k which relates the specific speed of a 
centrifugal pump, N,, to suction specific speed, S,,, as 


N 
at 


for a pump with a discharge measured in litres per minute. 


Answer: 176 
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Glossary of terms 


Absolute pressure 


Pressure measured above a vacuum. 


Absolute roughness 
Average height of undulations and imperfections on the inner surface of a pipe 
wall. 


Angular velocity 


Rotational.speed expressed in radians per second. 


Archimedes' principle 
States that the upthrust on a partially or totally immersed body is equal to the 
weight of liquid displaced. 


Average velocity 
Defined as the total flowrate per unit flow area. 


Boundary layer 
Region between a wall and a point in the flowing fluid where the velocity is at a 
maximum. 


Cavitation 
Destructive collapse of vapour in liquid in regions of high pressure. 


Centrifugal pump 
Mechanical device used to transport fluids by way of an enclosed rotating 
impeller imparting velocity to the fluid. 


Coefficient of contraction 
Ratio of area of vena contracta to orifice area. 
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Coefficient of discharge 
Ratio of actual flowrate to theoretical flowrate through an opening or 
restriction. 


Coefficient of velocity 
Ratio of actual to theoretical velocity for a discharging jet of fluid. 


Density 
Measure of the quantity of a substance per unit volume. 


Dynamic viscosity 
Phenomenon in which liquids shear between adjacent layers. 


Dynamical similarity 
Similarity between systems when applied forces have the same ratio. 


Elbow 
A curved 90° section of pipework. 


Entrance and exit losses 
Irreversible energy loss when fluids enter or leave pipes. 


Entrance transition length 
Distance between pipe entrance and point of fully developed laminar flow, 


Equivalent hydraulic diameter 

Used in place of diameter d in turbulent flow equations when the pipe or duct 
cross-section is not circular. Equal to four times the flow area divided by the 
wetted perimeter. Does not apply to laminar flow. 


Equivalent length 
The length of straight pipe that would give the same pressure drop as the fitting 
of the same nominal diameter expressible in pipe diameters of the straight pipe. 


Friction factor 
Empirical factor and is a function of Reynolds number and relative pipe 
roughness. 


Fluid 
A substance which offers no resistance to change of shape by an applied force. 
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Forced vortex 


Circular stream of liquid the whirl of which is caused by an external source such 
as a stirrer or impeller. 


Free jet 
Discharging liquid from a pipe or orifice which is no longer influenced by the 
applied head. 


Free surface 
Interface of a liquid with a vapour or gas above. 


Free vortex 
Streamlines that move freely in horizontal concentric circles with no variation 
of the total energy across them. 


Frictional head 
Head required by a system to overcome resistance to flow in pipes and associ- 
ated fittings. 


Gate valve 
Device used to regulate flow in a pipe, consisting of a vertical moving section 
across the flow area. 


Gauge pressure 
Pressure measured above atmospheric. 


Geometrical similarity 
Systems which have similar physical dimensions such that their ratio is a 
constant. 


Head 
Pressure of a liquid expressed as the equivalent height a column of the liquid 
would exert. 


Laminar flow 


Flow that occurs when adjacent layers of fluid move relative to one another in 
smooth streamlines. Flow occurs in pipes for Reynolds numbers below 2000. 


Mass flowrate 
Mass flow of a fluid per unit time. 
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Mass loading 
Rate of mass flow per unit area. 


Mean hydraulic depth 
Ratio of flow area to wetted perimeter in an open channel. 


Net positive suction head (NPSH) 
Total suction head determined at the suction nozzle of a centrifugal pump less 
the vapour pressure. 


NPSHA 
NPSH available is the absolute pressure at the inlet of a centrifugal pump and is 
a function of elevation, temperature and pressure. 


NPSHR 
The NPSH required by a centrifugal pump is a function of the pump character- 
istics and must be exceeded by NPSHA. 


Newtonian fluids 
Class of fluids in which viscosity is independent of shear stress and time. 


Non-Newtonian liquids 
Liquids in which the viscosity depends on shear stress and/or time. 


No-slip condition 
Assumption that the fluid layer in contact with a surface is stationary. 


Notch 
Rectangular or triangular cut in a weir. 


Open channel 
Conduit carrying liquid with a free surface. 


Orifice 
Small hole or opening. 


Orifice meter 
Device used to measure flowrate by the pressure drop through a small hole. 
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Piezometer 
Device used to measure pressure, consisting of an open vertical (glass) tube 
attached to the side of a horizontal pipe. 


Pipeline 
Long section of large bore pipe. 


Pitot tube 
Device used to measure velocity head of a flowing fluid. 


Power 
Work done per unit time. 


Pressure 
Force of a fluid applied over a given area. 


Pressure differential 


Difference in pressure between two points. 


Pressure head 


Pressure-volume energy of a fluid expressed in head form. 


Relative roughness 


Ratio of absolute pipe wall roughness to pipe inside diameter in consistent units. 


Reynolds number 
Dimensionless number expressing the ratio of inertial to viscous forces in a 
flowing fluid. 


Rotameter 
Instrument used to measure flowrate by the vertical elevation of a float in a 
tapered tube. 


Separation 

Phenomenon of fluid streamlines changing direction due to changes in 
boundary shape as a result of fluid inertia or velocity distribution near the 
boundary surface. 


Separator 
Vessel used to separate two liquids, usually as layers. 
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Shear stress 
Force of a liquid applied over a surface. 


Shear rate 
Change in velocity of a fluid perpendicular to flow. 


Sluice 
Submerged opening. 


Specific gravity 


Ratio of mass of a liquid to mass of water occupying the same volume. 


Specific speed 
Dimensionless group used to classify centrifugal pump impellers at optimal 
efficiency with respect to geometric similarity. 


Static head 
Potential energy of a liquid expressed in head form. 


Steady flow 
No change of fluid velocity with respect to time. 


Streamline 
An imaginery line which lies in the direction of flow. 


Siphon 
Transfer of liquid from one vessel to another at a lower elevation by means of a 
pipe whose highest point is above the surface of the liquid in the upper vessel. 


Transition flow 
Flow regime between laminar and turbulent flow (2000<Re<4000). Velocity 
fluctuations may be present and impossible to predict. 


Turbulent flow 

Flow regime characterized by fluctuating motion and erratic paths above 
Reynolds numbers of 4000. Flow occurs when inertial forces predominate 
resulting in macroscopic mixing of the fluid. 


Uniform flow 
No change in fluid velocity at a given time with respect to distance. 
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Velocity gradient 
See ‘Shear rate’. 


Velocity head 


Kinetic energy of a fluid expressed in head form. 


Vena contracta 
Minimum flow area formed downstream by a fluid flowing through a sharp- 
edged opening. 


Venturi meter 
Device used to measure flowrate by the pressure drop through a tapered section 
of tube or pipe. 


Viscosity 


See ‘Dynamic viscos 





Viscous flow 
See ‘Laminar flow’. 


Volute chamber 
Spiral casing of centrifugal pump surrounding impeller. 


Weir 
Vertical obstruction across a channel. 


Weir crest 
Top of a weir over which a liquid flows. 
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Nomenclature and preferred units 





Quantity Name 
Area 

Degree plane angle 
Density 

Force Newton 
Length metre 
Length kilometre 
Ma kilogram 
Mass flow 

Mass loading 

Pressure bar 
Pressure Pascal 
Pressure 

Radian plane angle 


Revolution angular displacement 
Shear stress 

Surface tension 

Time hour 
Time minute 
Time second 
Torque 

Velocity 

Viscosity 

Volume 

Volume litre 
Volume flow 

Work joule 
Watt power 
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Symbol 


m^ 


kgm? 
N 

m 

km 

kg 
kgs! 
kgms 
bar 
Pa 
Nm? 
rad 
rps 
Nm? 
Nm! 
h 

min 


Nm 





Type of 
SI unit 
Base 
Allowable 
Base 
Derived 
Base 
Allowable 
Base 

Base 

Base 
Allowable 
Derived 
Derived 


In terms of 
base SI 


T/180 rad 
kgms~? 


105 m 


10° kgm- 
kgm"! 
kgm 









Supplementary 


Allowable 
Derived 
Derived 
Allowable 
Allowable 
Base 
Derived 
Base 
Derived 
Base 
Allowable 
Base 
Derived 
Derived 


2n rad 














Useful conversion factors 


Quantity 


Acceleration 
cms? 
fts 





Area 
mm? 
cm? 
in? 
ie 
yd? 


Density 





Ibm.UKgal-! 
Ibm.USgal-! 


Energy 
kJ 

Nm 
Ibf-ft 


Force 
kN 
dyne 
lbf 
UKtonf 
UStonf 


Multiplier 


0.01 
0.3048 


0.000001 
0.0001 
0.0006451 
0.0929 
0.836 


1000 
1000 
16.018 
99.77 
119.83 


1000 
1 
1.356 


1000 
0.00001 
4.448 
9964 
8896 


SI or preferred unit 





N 
N 
N 
N 
N 
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Length Time 

mm 0.001 m 1 min 60 s 

cm 0.01 m I h 3600 s 
km 1000 m d 86.400 s 

in oe m Velocity (linear) 

ft i mn cms"! 0.01 ms! 
yd : m fis! 0.3048 ms! 
Mass aus " Velocity (angular) 

g in ha rps 0.10472 rads-! 
id iui i rpm 6.2832 ads 
oz(avdp 2: degs-! 0.017453 ads“! 
Ibm 0.4536 kg - : T 


Plane angle Viscosity x 











: 2 
deg) 0.01745 rad E d 
Power Ibf.s.i 6894.7 Nsm? 
1000 w cP 0.001 Nsm? 
1 w P 0.1 Nsm? 
d M Viscosity (kinematic) 
f fs! 0.092903 
Pressure cSt 0.000001 
Pa 1 Nm? St 0.0001 
kgui-s t l Nav? Volume 
mmH, 133.3 
ar 100,000 Nm m erre uf 
Std atm 101,300 Nm? Ld n T 
dynecm? 01 Nm ume a 
Ibfin- (psi) ^ 6895 Nm? m Q odii ma 
Pressure drop yd? 0.7645 m? 
Pam! 1 Nm?m! bbl(42USgal) 0.1589 m3 
psi.ft-! 0.022621 Nm?! UKgal 0.004546 n? 
Surface tension Usgal QS nm 
kgs! 1 Nm! Volume flow 
Ibf.f-! 14.59 Nm! mh"! 0.0002777 màs 
dyne.cm-! 0.001 Nm! fes 0.028317 màs! 
dyne.m-! 0.00001 Nn! bbl.h-! 0.0004416 mis 
UKgal.h-! 0.0045461 màs! 
USgal.h-! 0.003754 m3s-! 
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Physical properties of water 
(at atmospheric pressure) 


Temperature, Density, 


"E 

0 
5 
10 
15 
20 
25 
30 
35 
40 
45 
50 
55 
60 
65 
70 
75 
80 
85 
90 
95 
100 


kgm 
1000 
1000 
1000 
999 
998 
997 
996 
994 
992 
990 
988 
985 
983 
980 
978 
975 
972 
969 
965 
962 
958 


Viscosity, 
Nsm? 
0.00790 
0.00152 
0.00131 
0.00114 
0.00100 
0.00089 
0.00080 
0.00072 
0.00065 
0.00060 
0.00055 
0.00051 
0.00047 
0.00044 
0.00040 
0.00038 
0.00035 
0.00033 
0.00031 
0.00030 
0.00028 


Vapour pressure, Surface tension, 


Nm? Nm! 
611 0.076 
872 0.075 

1230 0.075 
1700 0.074 
2340 0.074 
3170 0.073 
4240 0.072 
8600 0.071 
7380 0.070 
9550 0.069 
12,300 0.068 
15,740 0.067 
20,000 0.067 
25,010 0.066 
31,200 0.065 
38,550 0.064 
47,400 0.063 
57.800 0.062 
70.100 0.061 
84,530 0.060 
101,300 0.059 


Losses for turbulent flow 
through fittings and valves 


Fitting or valve 

Elbow 45° standard 

Elbow 45° long radius 

Elbow 90° standard 

Elbow 90° long radius 

Elbow 90° square 

Bend 180° close return 

Tee, standard, along run 

Tee, standard, through branch 


Pipe entry sharp 
well rounded 
Borda 


Pipe exit sharp 


Gate valve fully open 
+, open 
!5 open 
!4 open 


Diaphragm valve — fully open 
?, open 
15 open 
!4 open 


Globe valve fully open 
15 open 

Angle valve fully open 

Butterfly valve fully open 


Loss coefficient, k 
0.35 

0,2 

0.75 

0.45 

1.3 

1.8 

0.4 

L5 


0.5 
0.05 
1.0 


0.5 


0.15 

0.9 

4.5 
20.0 


6.0 
9.0 
3,5 
0.5 
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Check valve, swing 
Pump foot valve 
Spherical plug valve fully open 


Flow meters orifice B=0.2 
B-0.5 
B=0.8 
nozzle B=0.2 
B=0.5 
B=0.8 
venturi 
turbine 
Cyclone 
290 


2-3.5 
1.5 
0.1 
2.6 
1.8 
0.6 
1.0 
0.55 
0.15 
0.15 
6.0 


10-20 


Equivalent sand roughness of pipes 


Type 
Glass 

Drawn tubing 

Steel, wrought iron 
Asphalted cast iron 
Galvanized iron, new 
Galvanized iron, old 
Cast iron 

Wood stave 
Concrete 

Riveted steel 


Equivalent sand roughness, £ (mm) 
0.0003 
0.0015 
0.0457 
0.1219 
0.1524 
0.2743 
0.2591 
0.1829 
0.3048-3.048 
0.9144-9.144 
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Manning coefficient for various 
open channel surfaces 


Type 
Welded steel 
Planed wood 
Unplaned wood 
Concrete culvert 
Unfinished concrete 
Cast iron 
Straight sewer 
Riveted steel 
Galvanized iron 
Earth channels clean 
stony 
unmaintained 
Natural streams clean 
stony 
deep pools 
Rivers no boulders 
irregular 
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Manning coefficient, n (m-/3s) 
0.012 
0.012 
0.013 
0.013 
0.014 
0.015 
0.015 
0.016 
0.016 
0.018 
0.035 
0.080 
0.030 
0.045 
0.100 
0.025 
0.035 


Moody plot 

















0025 
Lamar 
0.020 
0.05 
0.04 
0015 0.03 
0.02 
< 00s 3 
Zooo oo 3 
Š 0.009 008 5 
5 0.008 0.000 G 
£ 0.007 40004 £ 
z S 
Z 0.006 ow 2 
5 oo € 
E 00008 co 
u 0.0006 2 
£ oo 0.0004 E 
0.0002 $ 
0003 0000 2 
0.00005 
0.002 0.00001 
10 10* n 10° 107 108 


REYNOLDS NUMBER, Re 


Plot of friction factor versus Reynolds number developed by Louis F. Moody. 


Laminar flow (Re<2000) 
Transition region (2000<Re<3000) 
Turbulent region (Re>3000) 

For rough-walled pipes 
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(see Problem 8.11, page 217) 





Index 


A 
absolute pressure 1,2,5,6 
definition 275 
absolute roughness, definition 275 
acceleration 244, 246 
accumulators 4l 
angular velocity, definition 275 
ARCHIMEDES, 1,2 
Archimedes’ principle 25, 28, 29 
definition 275 
atmospheric pressure 1,2 
average velocity, definition 275 
axial pumps 235, 250 
B 
barometers 2 
bearings 93 
BERNOULLI, D. 35 
Bernoulli equation 42. 43, 45, 50, 


115, 119, 126, 129, 131, 
202, 207, 221, 240, 247 


BLASIUS, P.R.H. 193 
Blasius equation 200 
BORDA, J.C. 47 
boundary layer, definition 275 
bubbles 

pressure within 15 
BUCKINGHAM., E. 99, 106 
Buckingham II method. 106 
buoyancy 26 


€ 
calibration of flowmeters 124, 132, 134 
capacity coefficient 104, 105, 249, 





252, 258, 260 
capillary action 2 
CARNOT, L.N.M. 47 
cavitation 236, 244, 245, 
247, 259, 260, 261 
definition 275 
centrifugal force 51 
centrifugal pumps 237, 244, 245, 
249.256 
characteristic 252 
definition 
duty point 
efficiency 
in series and parallel 
matching 
scale-up 
channels 
circular 177, 179 
economical depth 175 
‘open 167-191 
rectangular 171-176 
CHEZY, A. 167 


Chézy coefficient 170, 171, 173, 177 
Chézy formula 169, 171, 173, 178 


circular channel flow 177, 179 
coefficient 
capacity 104, 105. 249, 252, 258, 260 
Chézy 170. 171, 173, 177 
contraction 47. 143. 146. 275 








discharge 122, 128, 134, 
143, 146, 276 
head 104, 105, 249, 252. 258, 
Manning 
power 104, 105 
pressure 108 
velocity 141, 143, 145, 276 
COLEBROOK, C.F. 194 
collar bearings 93 
conservation of 
energy 
mass 
momentum 
continuity 


continuity equation 
contraction coefficient 








definition 275 
conversion factors 285 
critical time 4l 
crystal concentration 

measurement of 13 
cylinder pumps 263 
D 
DARCY, H.P.G. 193 
Darcy equation 198, 199 
Darcy-Weisbach equation 198, 199 
da VINCI, L. as 
deceleration of liquid 40 

sity, definition 276 
diaphragm pumps 263 
differential head 

flow measurement by 113-139 
differential manometers 17, 20, 119 
differential pressure 2, 11, 13, 17, 23 

across a venturi meter 126 
dimensional analysis 99-111 

Buckingham II method 99, 106 

for centrifugal pumps 104 

for notches 102 

for orifices 100 

Rayleigh method 99. 100 
dip legs 





INDEX 








density measurement 11-12 

interface measurement 11-12 
discharge 

from reciprocating pumps 265 
discharge coefficient 128, 134, 143, 146 

definition 276 
drainage 

between two reservoirs 153 
drainage from tanks 

cylindrical cross-section 151 

hemispherical cross-section 147 

uniform cross-section 147 
ducts 118 

non-circular 223 
duty point 240, 242, 254 
dynamical similarity, definition 276 
dynamic pumps 235 
dynamic viscosity, definition 276 
E 
economic pipe diameter 196 
economical depth of flow 175 
eddies 44 
elbows, definition 276 
energy 
enlargements 44 
entrance losses, definition 276 
entrance transition length 65 

definition 276 
equations 


Bernoulli 42, 43, 45, 50, 115, 119, 
126, 129, 131, 202, 
207, 221, 240, 247 


Blasius 200 
Darcy 198, 199 
Darcy-Weisbach 198, 199 
Fanning 198 
Hagen-Poiseuille 68, 70, 71, 73 
Prandil 202 
Torricelli 141, 143, 145 
von Kármán 214 
equivalent hydraulic diameter 224, 225 
definition 216 
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equivalent length, definition 276 
equivalent length method 214,217 
EULER, L. 35 
exitheadlosses — 46, 206, 207, 213, 221 
definition 276 
F 
FANNING, J.T. 195 
Fanning equation 198 
Fanning friction factor 195 
film thickness in a channel 79 
flow 
charts for 293 
down an inclined plate 81 
down a vertical wire 82 
in branched pipes 36 
in circular channels 177, 179 


in rectangular channels. 171, 173, 175 


in tapering pipes/tubes 71, 226 
laminar 61-97 
multiphase fluid 103 
open channel 103 
over notches 102 
pulsed 235 
radial 51 
through a gap 84, 88 
through an orifice 100 
through non-circular ducts 223 
turbulent 106, 223 
velocity through a gap 86, 87 
vertical in pipes 78 
flow measurement 
by differential head 113-139 
flowmeters 113 
calibration 124 
calibration curve 132 
nozzle 113 
orifice 113, 114, 122, 128, 278 
Pitot tubes 113, 115, 117, 279 
rotameter calibration 134 
variable area 130 
venturi 113, 119, 122, 124, 
126. 134, 281 
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fluids 
definition 
pressure of 
fluid statics 
flumes 
force 
centrifugal 
inertial 
on pipe reducers 
on U-bends 
radial 
viscous 
within a hydraulic ram 
forced vortices 
definition 
free jets 
definition 
free surfaces, definition 
free vortices 
definition 
frictional head, definition 
friction factor, definition 


G 

gas-liquid interface 
gate valves, definition 
gauge pressure 








277 
52.53 
27 
277 
276 


79 

277 

1, 2, 268 
277 


geometrical similarity, definition 277 


gyroscopic effect 


H 
HAGEN, G.H.L. 


130 


61 


Hagen-Poiseuille equation. 68, 70, 71, 73 


head 
definition 
entrance loss 
inertia 
loss due to friction 
pressure 
static 

head coefficient 
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47 

229 

197 

43 

43 

104, 105, 249, 
252, 258, 260 





HERSHELS. C. 119 
hydraulic rams 7 
hydrometers 27 
1 
impact pressure 115 
impellers 237.256 
inclined leg manometers 2 
inertia head 229 
inertial force 
interface 
79 
liquid-liquid 9.11 
inverted manometers 19 
J 
JOUKOWSKY, N. 4l 
K 
kinematic viscosity 63 
L 
laminar boundary layer 65 
laminar flow 61-97 
definition 277 
down a wire 82 
general equation for pipe friction 199 
local velocity 68, 75, 86 
maximum pipe diameter T 
plate 81 
through tapered tubes 7 
through a gap 84 
LANGHAAR, H. 65 
liquid-liquid interfaces 9,11 
lubrication 61-97, 235, 268 
collar bearings 93 
shafts 91 
M 
MANNING. R. 167 
Manning coefficient 292 
Manning formula 170, 171, 173 
manometers 2, 17,115 





INDEX 
differential 17, 20, 119 
inclined leg 23 
inverted 19 
single leg 21 
U-tube 17 

mass flowrate, definition 277 

mass loading 63 
definition 278 

mean hydraulic depth 168, 172, 173, 

175,177 
definition 278 

menisci 2 

Metric system 1 

momentum 35-59 

Mono pumps 268 

MOODY, L.F. 194 

Moody plot 205, 218, 223, 225, 293 

multiphase fluid flow 103 

N 

net positive suction head (NPSH) 245, 

247,251 
definition 278 

Newtonian fluids, definition 278 

NIKURADSE, J. 193, 194 

non-circular ducts 223 


non-Newtonian fluids, definition 278 


no-slip condition 65, 67, 79, 82 


definition 278 
notches 167-191 
definition 278 
flow over 102 
rectangular 180 
trapezoidal 189 
triangular 186 
v 188 
nozzle flowmeters 113 


NPSH (see net positive suction head) 


o 

open channels 
definition 278 
flow 103 
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orifice flow under constant head 
orifice meters 








143 


113, 114, 122, 128 


definition 278 
orifices 
coefficient of discharge 122, 128 
definition 278 
dimensional analysis 100 
flow through 100 
P 
PANNELL, J.R, 193 
PASCAL, B. 1 
Pascal's law 3 
piezometers 2 
definition 279 
pipe diameter 196 
for laminar flow 70 
pipe enlargement 44 
pipe friction 193-233 
pipelines, definition 279 
pipe reducers 49 
pipes 
acceleration in 228 
in pa 209 
211,214 
rough-walled 204, 214 
smooth-walled. 202 
pistons, reciprocating 264 
PITOT, H. 115 
Pitot traverse 117 
Pitot tubes 113, 115, 117 
definition 279 
PLIMSOLL, S. 30 
POISEUILLE, J.L.M. 61 
positive displacement pumps 235, 268 
power, definition 279 
power coefficient 104, 105 
power for pumping 264 
PRANDTL, L. 193 
Prandtl equation 202 
pressure 
absolute 
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atmospheric 

at a point 

definition 

differential 2, 

gauge 

head 

impact 

static 

vacuum 1 

vapour 206, 244, 245, 251 

within a closed vessel 5 

within a gas bubble 15 
pressure coefficient 108 
pressure differential, definition 279 
pressure drop 

due to enlargements 44 

in a pipeline, scale model 108 
pressure head 43 

definition 279 


pressure measurement 
differential manometers 17, 20, 119 
inclined leg manometers 
inverted manometers 
single leg manometers 
U-tube manometers 

pressure rise by valve closure 

pulsation dampers 

pulsed flow 




















pump characteristic 243, 

pumps 
axial 
centrifugal . 239, 243, 244, 245, 

249, 254, 256, 258, 275 

characteristic curve 243, 253, 255 
cylinder 263 
diaphragm 263 
duty point 2, 254 
dynamic 235 
efficiency 
Mono 
positive displacement 
power 
ram 











reciprocating 235, 246, 262, 
rotary 235, 246, 








265 
268 


scale-up of centrifugal 104 
specific speed 249 
suction specific speed 260 
R 
radial force 51 
ram pumps 263 
RAYLEIGH, Lord 99, 100 
Rayleigh method 107 
reciprocating pistons 264 
reciprocating pumps 235, 246, 262 
discharge from 265 
rectangular channels 171-176 
rectangular notches 180 
rectangular weirs 182. 184. 
relative roughness, definition 279 
REYNOLDS, O. 61 
Reynolds number 63, 70, 76, 77, 
80. 105, 108, 199 
defi 279 
ring mains 219 
rotameters 114, 130, 134 
definition 279 
rotary pumps 235, 246, 268 
rough-walled pipes 204,214 
ROUSE. H. 218 
5 
safety relief valves 268 
scale models 108 
scale-up of centrifugal pumps 104 
separation 244 
definition 279 
separators 9 
definition 279 
sewerage pipes 167 
shear rate, definition 280 
shear stress 93 
definition 280 
for flow through a gap 88 
similarity 


INDEX 
dynamical, definition 276 
geometrical, definition 277 

single leg manometers 21 

siphons 206 
definition 280 

sluices, definition 280 

smooth-walled pipes 202 

sound transmission 4l 

specific gravity 
definition 280 
measurement by hydrometer 27 

specific speed 249, 251, 259, 260 
definition 280 

STANTON, T.E. 193 

static head. 43 
definition 280 

static pressure 115, 245 

steady flow, definition 280 

STEVIN, S. 1 

streamlines 43, 51, 52, 54, 61, 143 
definition 280 

suction specific speed 260 


surface tension 15, 102, 181, 185 
Système International 1 


T 
tank drainage 141-165 
cylindrical cross-section 151 
hemispherical cross-section 147 
instantaneous. 157 
through a pipe 221 
through V-notches 188 
uniform cross-section 147 
with laminar flow 161 
tank inflow and outflow 155, 159 
tapered sections, head loss 226 
lapered tubes, laminar flow 7 
tension 26 
torque 89,91, 04 
TORRICELLI, E. 141, 146 
Torricelli equation 141, 143, 145 
tracer dilution 124 
transition flow, definition 280 
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transmission velocity 4l 
trapezoidal notches 189 
triangular notches 186 
turbulent flow 193-233 
definition 280 
U 
U-bends 38 
uniform flow, definition 280 
U-tube manometers 17 
V 
vacuum 12 
valves 
closure 40 
pressure relief 41. 268 
vapour bubbles 244, 246 
vapour pressure 206, 244, 245, 251 
variable area flowmeters 130 
VASCHY, A. 99 
velocity 
average 275 
of flow through a gap 86, 87 
of fluid in a pipe. 73,75 
velocity coefficient — 141, 143, 145. 276 
velocity gradient 75.88 
definition 281 
velocity head 43 
definition 281 
velocity head method 211,217 
velocity profile in a pipe 66 
vena contracta 47. 141. 143, 145 
definition 281 
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VENTURI, G.B. 119 
venturi meters 
calibration 124, 134 
coefficient of discharge 114 
definition 281 
horizontal 119, 122 
vertical 126 
vibration 238, 244 
viscometers 89 
viscosity, definition 281 
viscous drag 64,91 


viscous flow, definition 281 
viscous force 64 
viscous shear stress 75,88 
V-notches 188 
volute chamber, definition 281 
von KARMAN, T. 194 
von Kármán equation 214 
vortex motion 51 
vortices 

forced 52,53 

free 52,53 
w 
wall shear stress 88, 169, 197, 223 
water 

physical properties 288 
weir crest, definition 281 
weirs 103, 167-191 

definition 281 

rectangular 182, 184 
WEISBACH, J. 193, 198 
whirlwinds 52 
WHITE, C.M. 194 
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